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NOTATION 


A  =  wetted  cross-sectional  area  in  the 

physical  component,  L^ 

A  =  area  served  by  combined  sewers,  acres 

c 

A  =  total  area  of  catchment,  acres 

AR  =  urban  area  runoff,  inches  per  hour 

A(t)         =  fluctuating  cross-sectional  area  due  to  a 
variable  depth,  l2 

B  =  benthal  demand  of  bottom  deposits, 

mg/1-hour 

BOD         =  Biochemical  Oxygen  Demand 

BOD         =  mixed  BOD  concentration  in  the  combined 
^  sewer,  mg/1 

BCTDpM        =  average  field  measured  BOD  concentra- 


'CM 


'f 


m 


tion  in  combined  sewer  overflows 
mg/1 


BOD.        =  BOD  concentration  of  dry-weather  flow 
°  treatment  facility  effluent,  mg/1 

BOD^        =  BOD  concentration  of  municipal  sewage. 


mg/l 


BOD         =  mixed  BOD  concentration  in  receiving 


water,  mg/1 


BOD         =  BOD  concentration  of  urban  stormwater 
^  runoff,  mg/1 
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BOD 

BOD^ 
BOD^ 
BOD 

V 

C 
G 


sc 


mm 
c .  V  . 


mixed  BOD  concentration  of  urban  runoff 
from  the  separate  storm  sewer  and  com- 
bined sewer  flows,  mg/1 

hourly  BOD  concentration  of  total  urban 
runoff,  mg/1 

mixed  BOD  concentration  from  sources  up- 
stream of  urban  area,  mg/1    , 

BOD  concentration  of  treated  WWF  storage/ 
treatment  facility  effluent,  mg/1 

concentration  of  water  quality  parameters 
(pollutant),  M/l3 

concentration  of  D.O.  in  the  stream, 
mg/1 

value  of  input  concentration  from  t. 
to  (t.  +  At),  M/l3  V 

concentration  of  D.O,  at  maximum  deficit,: 
mg/1 

dissolved  oxygen  saturation,  mg/l 


sample  coefficient  of  variation  of  the 
water  quality  parameter 

3 
constant  influent  concentration,  M/L 

factor  to  convert  mg/1  to  Ibs/cf 


^2 
CR. 


C,(I) 


=  concentration  of  pollutant  in  the  inflow, 
dimension! ess 

=  factor  to  convert  FF  from  Ib/hr  to  cfS" 
mg/1 

3 
=  mean  effluent  concentration,  M/L 


composite  runoff  coefficient  dependent 
on  urban  land  use 

concentration  of  wastewater  in  the  storage/ 
treatment  system  at  t  =  0,  for  storm  event 
I,  M/l3 


XV 


C2(I-1) 


Ci(t) 


concentration  of  wastewater  in  the  system 
during  the  last  hour  of  runoff  of  the  pre- 
vious storm  event,  M/L-^ 

concentration  forcing  function,  dimension- 
less  : 


C,(T.) 


discrete  concentration  forcing  function, 
d i men  si  on  1  ess 


c^(t) 


concentration  forcing  function  of  material 
in  inflow,  M/L"^ 


C2(t) 


C2(t.) 


At 


C2(t.  +  ^) 


CgCx,  t  -  t^  )  = 


C^(x-Ut, 
t-^V 


D  =  Cg  -  C 


concentration  of  material  in  the  tank  and 
outflow,  as  a  continuous  function  of  time, 

concentration  value  obtained  by  equation 
(5.12)  for  the  previous  time  step,  M/l3 

time-averaged  concentration  for  each  time 
step,  M/l3 

pollutant  concentration  at  distance  x  along 
the  flow  axis  and  time  t  >^  t,  ,  mass  of 
pollutant/mass  of  fluid,  dim4nsionless 


initial  concentration  at  a  time  x/U 
earlier  and  a  distance  Ut  upstream, 
M/L3 

dissolved  oxygen  deficit,  mg/1 

critical  (maximum)  D.O.  deficit,  mg/1 

initial  D.O.  deficit,  mg/1 


D 
u 

D.O. 
DWFSEP 


D.O.  deficit  in  receiving  waters  upstream 
of  inflow  point,  mg/1 

Dissolved  Oxygen 

DWF  contribution  from  seoarate  sewer  area 
cfs 


DWFCMB 


DWF  contribution  from  combined  sewer  area, 
cfs 
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DWH 

DWT(I) 

E 

f  =  K2/K^ 

FF 

f(t  -  t) 

g 

H 
h 

I(t) 


number  of  dry  weather  hours  preceding  each 
runoff  event 

number  of  DWH  preceding  storm  event  I 

longitudinal  dispersion  coefficient,  L^/T 

self-purification  ratio,  dimensionless 

available  urban  depression  storage,  inches 

first  flush  factor,  pounds/DWH-acre 

system  step  response,  dimensionless 

exponent  in  solution  to  steady-state  dif- 
ferential equation  for  dissolved  oxygen 
deficit,  T-1 

stream  depth,  feet 

sampling  interval  =  t^^,  -  t.  =  1  hour 

3 
system  input  function,  M/L 

exponent  in  solution  to  steady-state  dif- 
ferential equation  for  dissolved  oxygen 
deficit,  T-i 

first-order  decay  rate  of  pollutant  in  the 
fluid  medium,  T~ ' 

number  of  hourly  lags 

oxidation  coefficient  of  nitrogenous  BOD, 
hours-! 


first-order  BOD.  decay  constant ,  day 
and  hours-! 
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atmospheric  reaeration  coefficient,  hours 

length  of  the  storage/treatment  basin, 
feet  : 

remaining  carbonaceous  BOD  concentration, 
mg/1 

ultimate  first-stage  BOD  demand,  mg/1 
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Pe 
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At 
At, 


equalized  variation  of  the  flow,  dimen- 
sion 1  e  s  s 

exponent  in  the  steady-state  differential 
equation  for  dissolved  oxygen  deficit, 
L-l 

total  number  of  increments  desired  per 
time,  step 

remaining  nitrogenous  BOD  concentration, 
mg/1  V 

any  desired  upper  bound,  dimensionless 
number 

total  number  of  data  points  or  observations 

total  number  of  inputs 

oxygen  production  rate  by  algal  photosyn- 
thesis, mg/1 -hour 

hourly  rainfal  1  /sndwmel  tin  inches  over  ~--"- 
the  urban  area 

Peclet  dispersion  number 

volumetric  flow  into  and  out  of  tank, 
l3/T 

mean  influent  fluid  flow  rate,  L'^/T 

combined  sewer  overflow  rate,  cfs 

DWF  treated  effluent,  cfs 

value  of  input  flow  rate  from  t.  to 
(t.  +  At),  l3/T  V 

urban  runoff  carried  by  the  separate 
storm  sewer,  cfs 

total  (storm  plus  combined)  urban  runoff, 
cfs 

upstream  flow,  cfs 
. WWF  storage/treatment  effluent,  cfs 
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QCOM 


q"(T) 


Q,(I) 


Qo(l  -  1) 


total  flow  generated  by  the  combined 
sev;er  area,  including  DWF  contribution, 
for  all  periods  of  urban  runoff,  cf/yr 

variable  time  rate  of  mass  injection  per 
unit  area  at  the  plane  source,  M/L^T 

system  outflow  rate  at  t  =  0,  for  storm  , 
event  I,  L^/T 

influent  fluid  flow  rate  per  unit  width, 
L2/T 

influent  fluid  flow  rate,  L^/T 

outflow  rate  from  the  storage/treatment 
system jduring  the  last  hour  of  runoff  of 
the  previous  storm  event,  L-^/T 

outflow  rate,  L-^/T 


Ql(t) 


Q2(t) 
Q2(t=0) 

Q2(t^.  +  At) 


QgCt^  +  ^) 


fluid  flow  rate  into  tank,  as  a  continuous 
function  of  time,  L^/T 


fluid  flow  rate  out  of  the  tank, 
continuous  function  of  time,  L-^/T 


as  a 


initial  outflow  of  the  system,  l^/l 

outflow  value  at  the  end  of  the  previous 
time  step,  L^/T 

output  flow  rate  at  the  end  of  the  time 
step,  l3/T 

time-averaged  outflow  rate  for  each  time 
step,  L3/T 

total  number  of  hourly  runoff  occurrences 
during  the  year 

fraction  removal  of  BOD  achieved  by  the 
DWF  treatment  facility 

algal  respiration  rate,  mg/1-hour 

sources  or  sinks  of  the  substance  C, 

m/l3t    : 
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s  .  =   standard  deviation  of  observations  of  the 

water  quality  parameter  from  i ts  computed 
mean 

s^  =  unbiased  estimate  of  the  variance  of  ob- 

served magnitudes  of  the  water  quality 
.     parameter 

T  =  water  temperature,  °C 

t  =  Time,  T 

t  =  elapsed  time  at  which  the  critical  deficit 

^  occurs,  hours 

t_  =  detention  time,  T 

t"p.  =  average  detention  time,  T 

t.  =  the  beginning  of  the  time  interval  for 

V  which  the  system  response  is  being  eval- 

uated, T 

t .  =  any  time  t  >^  t, 

T  =  correlation  time  constant  for  influent 

"  fluid  flow  rates,  hours 

T  =  correlation  time  constant  for  mass  rate 


w 


inputs,  hours 


t,   .        =  duration  of  injection,  T 

TDWH        =   total  number  of  DWH  in  the  year,  DWH/yr 

TL  (95%)     =   tolerance  limits  at  a  95%  probability 
level 

U  =  longitudinal  velocity  in  the  storage/ 

treatment  system,  or  in  the  stream,  L/T 

'    3   ■  '   '    ■ 

V  =  volume  of  tank,  L 

3 
V{t)         =   time-varying  volume  of  the  basin,  L 

W.  =  value  of  input  mass  rate  from  t.  to 

^  (t.  +  At),  M/T 


XX 


mean  influent  mass  rate,  M/T 


-  /4.  J.  At\ 
W2{t.  4-  -^) 


X  . 

1 


^1 

h 

At 

V,. 


'CQ 


time-averaged  output  mass  rate  for  each 
time  step,  M/T 

unbiased  estimate  of  the  mean  value  of 
the  water  quality  parameter 

=  discrete  data  series  (observations)  of  a 
hydrologic  process  I,  for  i  =  1,  2,  .  .  . 
n 

=  unbiased  estimate  of  the  flow-weighted 
mean  value  of  the  pollutant  concentration 

=  time  constant,  T~ 

=  regression  coefficient 

regression  CO effi  cient 
-constant  of  integration 
=   regression  coefficient 

regression  coefficient 

=   length  of  time  interval,  say,  1  hour 

=  coefficient  of  variation  of  the  effluent 
concentration,  dimensionl ess 

=  coefficient  of  variation  of  the  effluent 
mass  flow  rate,  dimensionless 

coefficient  of  variation  of  the  influent 
fluid  flow  rate,  dimensionless 

=  coefficient  of  variation  of  the  influent 
mass  rate,  dimensionless 


=  /u  +4KE,  has  the  dimensions  of  velocity, 
L/T 

3 
=  density  of  the  receiving  fluid,  M/L 

standard  deviation  of  the  influent  fluid 
flow  rate,  l3/T 

standard  deviation  of  the  influent  mass 
rate,  M/T 
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variance  of  the  influent  fluid  flow  rate, 
L6/T2 

2   2 
variance  of  the  influent  mass  rate,  M  /T 


-^i 


-^Q 


=  dummy  variable  of  integration 

iAt,  dimensions  of  time,  T 

=   normalized  detention  time  for  fluid  flow, 
dimensionless 

=   normalized  detention  time  for  pollutant 
mass  input,  dimensionless 

data  series  of  observed  magnitudes  of 
the  water  quality  parameter,  for  i  =  1, 
2,  3,  .  .  . ,  n 
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Pollutant  transport  within  each  phase  of  the 
hydrologic  cycle,  and  through  the  various  components 
of  the  physical  system,  is  governed  by  the  principle  of 
conservation  of  mass.   Deterministic  mathematical  models 
are  derived  from  this  unifying  concept  to  represent  the 
movement,  decay,  storage,  and  treatment  of  stormwater 
runoff  pollutants  and  dry  weather  wastewater  flows 
through  the  urban  environment  and  ,the  receiving  body  of 
water.  The  general  one-dimensional ,  transient  conserva- 
tion of  mass  equation  may  be  simplified  for  application 
to  the  various  pathways  by  retention  of  only  the  dominant 
terms  in  each  instance.   The  detention  time  and  decay 
coefficient  are  key  parameters  in  establishing  the  final 
form  of  the  governing  differential  equation. 
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The  precipitation  time  series  and  the  urban  runoff 
series  generated  by  a  continuous  hydrologic  simulation 
model  are  characterized  stochastically  by  autocorrelation 
analysis  to  define  independent  wet  weather  events.   The 
transient  response  of  storage/treatment  systems  to  vari- 
able forcing  functions  of  flow  and  concentration  is  deter- 
mined for  completely  mixed  systems  of  constant  and  vari- 
able volumes  and  for  one-dimensional  advective  systems 
with  and  without  dispersion,  and  results  are  compared. 
Frequency  analyses  are  performed  on  input  and  output 
concentrations  and  mass  rates  for  a  single  event  and  for 
al 1  wet  weather  events  during; the  year  of  record/   .- 

The  unifying  concept  of  continuity  is  extended 
to  determine;  the  receiving  water  response  to  waste  inputs 
from  (1)  wet  weather  urban  sources,  (2)  dry  weather 
urban  sources,  and  (3)  upstream  sources.   The  results 
are  presented  in  terms  of  minimum  dissolved  oxygen  cumu- 
lative frequency  curves.   Interpretations  are  based  on 
established  stream  standards  for  the  study  area:   Des 
Moinesv  Iowa  and  the  Des  Moines  River.   The  wet  weather 
flow  storage/treatment  facility  for  the  receiving  water 
analysis  was  represented  by  a  well-mixed  constant  volume 
model.   The  receiving  water  response  was  found  to  be 
sensitive  to  the  length  of  the  detention  time  in  the 
storage/treatment  device,  during  periods  of  urban  runoff. 
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CHAPTER- I 
INTRODUCTION 

Historically,  scientific  concepts  have  been  modes  of 
reasoning  rather  than  absolute  realities.   In  this  sense, 
models  are  tools  of  thinking  to  formulate  hypotheses  which 
must  be  verified  by  observations  (Vollenweider,  1975).   The 
nature  of  the  models  themselves  depends  largely  on  points  of 
view  and  objectives.   Problems  involving  physical  phenomena 
seem  to  fit  naturally  into  three  classes  (Domenico ,  1 972)  : 

(1)  those  studied  from  a  scientific  point  of 
view  in  order  to  promote  a  better  de- 
scription of  physical  phenomena, 

(2)  those  studied  from  an  engineering  point  of 
view  to  achieve  useful  purposes,  such  as 
prediction  by  using  the  laws  of  science, 
and 

(3)  problems'  studied  from  a  management  and 
planning  pointofview  to  achieve  some 
degree  of  control  over  the  state  of  the 
system. 

Scientific  studies  are  based  on  the  premise  that  all  natural 
phenomena  and  processes  are  interrelated  and  interactions  are 
governed  by  certain  laws,  thus  involving  numerous  subsystems. 
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The  research  problem  addressed  in  this  work  is  the  interaction 
of  urban  stormwater  runoff,  control  measures,  and  receiving 
water  response.   The  water  quality  cycle  is  a  dynamic  system 
interacting  with  each  phase  of  the  hydrologic  cycle.   The 
amount  of  pollution  entering  or  leaving  a  water  body  is 
determined  by  the  quantity  of  flow  and  concentration  of 
pollutants  in  each  of  the  hydrologic  components  of  the  physi- 
cal system.   The  retention  of  pollutants  in  the  water  body 
is  not  a  sole  function  of  the  quantity  and  quality  of 
hydraulic  .flows,since  it  also  depends  upon  the  location  of 
the  pollutants  within  the  water  body.   The  pollutants  exist 
in  the  water,  bottom  sediments,  and  the  aquatic  life. 
The  engineering  problem  is  finding  analytical 
relationships  between  the  variables  characterizing  the  inflow 
and  outflow  processes  and  parameters  defining  the  state  of 
the  system.   In  principle,  these  analytic  relationships  are 
provided  in  all  cases  by  solution  of  the  complete  equations 
of  energy,  mass,  momentum,  and  state  (Eagleson,  1970). 
However,. it  is  seldom  possible  to  formulate  all  of  these 
equations  accurately  (Eagleson,  1970)  because  of 

(1)  inadequate  knowledge  of  physical  behavior, 

(2)  unknown  system  heterogeneities  and 
anisotropics, 

(3)  unknown  temporal  variabilities  of  system 
parameters ,  and 

(4)  numerical  approximations  introduced  for 
computational  economy,  or  to  obtain  a 


solution  where  direct  mathematical  analysis 

is  impossible. 
Whereas  state  variables  (density,  volume,  temperature,  etc.) 
define  the  condition  of  system  components,  decision  vari- 
ables  act  to  modify  the  state.   For  example,  storage  and 
treatment  may  modify  the  concentration  of  pollutants  in  an 
accelerated  manner  to  prevent  damaging  shock  loadings  from 
entering  receiving  bodies  of  water.   The  degree  of  control 
defines  the  management  problem. 

Urban  hydrology,  including  water  quality  processes, 
is  a,  combination  of  concepts  and  parameters  that  pertain 
to  scientific,  engineering,  and  management  points  of  view. 
The  essence  of  a  rational  approach  to  water  quality  control 
in  the  urban  environment  is  the  development  of  a  conceptual 
model,  based  on  scientific  principles,  which  has  the  pre- 
dictive capabilities  necessary  to  the  decision-making  pro- 
cess.  The  high  cost  of  pollution  control  faci 1 i ties ,  in 
terms  of  energy  utilization,  land  requirements,  engineer- 
ing manpower,  and  long-term  financial  burden,  obligates 
the  planning  agency  to  select  an  optimum  strategy  for  area- 
wide  wastewater  management.   Such  a  process  must  focus  on 
a  systematic  procedure  that  identifies  and  defines  (1) 
the  cause/effect  relationships  of  the  physical  environ- 
ment, (2)  the  efficiency  of  control  alternatives,  and  (3) 
the  benefits  to  be  derived  from  implementation  of  these 
controls.   The  mathematical  models  applied  need  not 


incorporate  all  phenomena  but  rather  should  be  relevant 
to  the  problem  under  consideration. 

Two  philosophies  have  emerged  in  recent  years  to 
provide  decision  perspectives  on  urban  stormwater  pollu- 
tion:  (1)  the  micromodel ing  approach,  and  (2)  the  macro- 
model  ing  approach.   The  more  complex  urban  hydrologic 
simulation  models  typify  the  micromodeling  philosophy, 
where  short  time  increment  mathematical  models  are  applied 
to  a  specific  catchment  subjected  to  a  known  or  synthetic 
rainfall  history.   Unquestionably,  as  the  degree  of  detail 
increases,  so  does  the  difficulty  to  generalize  the  results 
obtained.   This  is  parti cul arly  true  of  sewer  routing  models, 
where  hydraulic  computations  depend  on  di verse  sewer  cross- 
sections.   For  example,  13  different  sewer  cross-sections 
are  represented  in  the  U.S.  EPA  Storm  Water  Management 
Model  (Huber  et  al.  ,  1975).   Thus,  varying  shapes  and 
dimensions  from  urban  area  to  urban  area  will  modify  the 
runoff  hydrographs.   These  models  are  invaluable  in  iden- 
tifying problem  areas  in  sewer  networks.   The  Storm  Water 
Management  Model  can  predict  the  elements  where  surcharging 
will  occur  and,  as  an  option,  enlarges  the  downstream 
conduit  by  standard  amounts  until  capacity  exists  to  accept 
the  flow  (Heaney  e_t  al.  ,  1975).   However,  experience  at 
the  University  of  Florida  in  applying  the  detailed  simula- 
tion models  to  large  metropolitan  areas  suggests  that  the 
data  needs  alone  may  be  quite  substantial  and  the 


discretization  for  the  mathematical  abstraction  of  the 
physical  drainage  system  is  time  consuming  (Medina,  1974). 

At  the  other  end  of  the  mathematical  modeling 
spectrum  is  the  macroapproach .   The  methodology  is  typi- 
fied by  a  simplified  version  of  the  physical  process  to 
compute  a  scalar  quantity  which  ranks  pollution  severity 
and  also  allows  study  of  the  scalar's  sensitivity  to  the 
interrelationships  involved  (Young,  1976).   The  approach 
focuses  on  an  aggregated  receiving  water  pol 1 utograph, 
specifically  the  ultimate  BOD  concentration  versus  time. 
This  pollutograph  is  assumed  to  represent  water  quality 
downstream  from  an  urban  area  as  a  flow-weighted  average 
of  concentrations  from  sewage  residuals,  urban  storm 
washoff  and  upstream  sources.   As  computed,  it  is  indepen- 
dent of  hydrologic  rainfall -runoff  features,  and  varies 
with  the  receiving  water  streamflow  and  the  public  treat- 
ment works  of  the  urban  system.   Decision  rules  are  based 
on  calculation  of  scalar  pollutograph  maxima  for  ranking 
the  control  needs  of  the  various  urban  areas  studied. 
This  approach  has  considerable  merit  because  of  its  sim- 
plicity of  application.   However,  there  are  several 
limitations; 

(1)   estimates  for  design  streamflow,  urban 
runoff  quality,  removal  efficiencies, 
and  hydraulic  capacity  of  interceptor 
sewers  are  largely  dependent  on  the 
analyst's  judgment  and  experience; 


(2)  \ independence  of  the  methodology  to  vari- 

ations in  volume  of  runoff  and  associated 
pollutants  excludes  the  possibility  of 
evaluating  realistically  the  relative 
efficiency  of  control  measures  and  the 
receiving  water  response  to  these  mea- 
sures ;  and 

(3)  a  single  scalar  index  precludes  obtaining 
any  information  about  the  frequency  of 
occurrence  of  events  that  result  in  dele- 

:  terious  water  quality  conditions  in  the 
receiving  waterways.  ;;    :   ' 

There  is  a  real  need  for  , an  approach  which  is 
intermediate  between  the  mi cromodel i ng  and  macromodel i ng 
philosophies,  provides  for  continuous  hydrologic  simula- 
tion, and  recognizes  that  pollutant  transport  through  the 
various  components  of  the  urban  environment  is  governed 
by  the  same  principle  that  describes  transport  in  the 
natural  environment--conservation  of  mass.  The  justifica- 
tion for  continuous  hydrologic  simulation  in  dealing  with 
problems  of  urban  storm  runoff  and  urban  storm  runoff 
quality  is  the  probability  of  occurrence  of  events  of 
various  magnitudes  (Linsley  and  Crawford,  1974).   Conse- 
quently, the  research  objectives  are 

(1)   to  present  a  unified  physical  system 
conceptualization; 


(2)   to  develop  a  model  based  on  the  unifying 

concept  which  is  representative  of  the  land 
use,  hydrology,  and  climatology  of  the 
drainage  area  while  providing  an  analytical  , 
framework  to 

-generate  stormwater  runoff  pollutant  loads 
and  dry  weather  sanitary  flow  pollutant 
1 oads , 
-simulate  the  pollutant  removal  capabilities 

of  various  storage/treatment  alternatives, 
-simulate  the  conveyance  system,  including 
mixing  in  combined  sewers  of  wet  and  dry  . 
weather  pollutants, 
-mix  the  various  pollutant  inflow  combina- 
tions with  pollutants  already  in  the  re- 
ceiving water  (from  upstream  sources) ,  and 
-predict  the  oxygen  balance  of  the  polluted 
waters  downstream  from  the  waste  sources, 
subject  to  the  constraints  imposed  by  the 
quantity  and  quality  of  the  data  base;  and 
(3)   to  examine,  in  particular,  the  response  of 
storage/treatment  systems  to  variable  wet 
weather  flows  by  a  detailed  mathematical 
application  of  the  continuity  equation. 
Scientific  principles  derived  from. hydrologic 
theory,  chemical  reactor  engineering  and  process  dynamics 


unit  operations  of  wastewater  engineering,  and  the  mathe- 
matics of  diffusion  are  applied  to  characterize  urban 
runoff  quantity  and  quality,  evaluate  control  measures, 
and  determine  receiving  water  response.   A  major  portion 
of  this  work  is  devoted  to  storage/treatment  systems 
because  of  their  importance  in  reducing  pollutant  loads 
to  the  receiving  waters  from  highly  variable  wet  weather 
flows.   These  systems,  as  stated  in  the  objectives,  are 
represented  by  mass  transport  models  satisfying  continuity. 
Reservoir  management  models  have  been  proposed  using 
inventory  theory  (Sobel,  1975),  but  these  models  are  best 
suited  to  formulate  release  policies  which  are  a  function- ; 
of  the  quantitative  demand  for  water  supply.   It  is  diffi- 
cult to  foresee  their  application  in  a  context  where  water 
quality  and  interaction  between  the  urban  environment  and 
the  receiving  water  are  important. 

Chapter  II  discusses  the  characterization  of  urban 
wastewater  discharges  and  polluted  surface  waters,  and 
presents  a  unified  physical  system  conceptualization.   An 
oxygen  demand  parameter  is  selected  to  characterize  the 
strength  of  all  waste  sources  in  the  system.   The  principle 
of  conservation  of  mass  is  shown  to  be  equally  applicable 
to  pipe  segments  of  a  sewer  system,  storage/treatment 
devices,  and  receiving  bodies  of  water. 

Chapter  III  describes  the  study  area,  including 
the  receiving  stream,  and  presents  the  necessary  modeling 


techniques.   Typical  pollutant  loads  are  summarized  and 
pertinent  water  quality  standards  are  defined.   Data 
sources  are  also  identified. 

Chapter  IV  presents  a  stochastic  characterization 
of  urban  runoff.   Techniques  of  time  series  analysis  are  . 
discussed,  as  well  as  an  approximate  graphic  procedure,  to 
define  a  wet  weather  event. 

Chapter  V  addresses  the  response  of  storage/treat- 
ment systems  to  variable  inputs  of  wet  weather  flows,  and 
their  associated  pollutant  concentrations  and  mass  rates. 
Storage/treatment  units  are  modeled  as  well-mixed  constant 
volume  systems,  well-mixed  variable  volume  systems,  and 
variable  volume  systems  where  nonideal  mixing  is  repre- 
sented by  advection  and  dispersion.   Numericalappl i cations 
are  shown  for  each  type  of  system  for  various  residence 
times,  and  statistical  summaries  of.  the  results  are  pre- 
sented. 

Chapter  VI  introduces  a  statistical  approach  to 
storage/treatment  response  for  the  case  where  the  system 
is  assumed  to  be  represented  by  a  well-mixed  constant 
volume  model.   The  results  are  compared  to  those  obtained 
from  solution  of  the  continuity  equation  in  Chapter  V.   A 
discussion  of  the  concentration  and  mass  flow  rate  equali- 
zation provided  by  the  well -mixed  constant  vol ume  model 
i  s  incl uded. 
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Chapter  VII  discusses  the  methodology  adopted  for 
the  receiving  water  analysis.   The  receiving  water  response 
to  waste  input  combinations  is  presented  in  the  form  of 
minimum  dissolved  oxygen  frequency  curves.   The  interpre- 
tation of  these  results  is  based  on  established  water 
quality  standards. 

Finally,  Chapter  VIII  summarizes  the  conclusions 
drawn  from  this  work  and  the  recommendations  made  for 
further  research. 


CHAPTER  II 
A  UNIFIED  PHYSICAL  SYSTEM  CONCEPTUALIZATION 

2.1   Characterization  of  Urban  Wastewater  Discharges 
and  Polluted  Surface  Waters 

Urban  areas  represent  only  3  percent  of  total  land 
use  in  the  Uni ted  States  (U.  S.  Dept.  of  Agriculture, 
1968),  yet  they  constitute  the  centers  of  most  intense 
human  activity.   Point  di  scharges  resulting  from  com- 
mercial, industrial,  and  residential  wastes  generally 
enter  the  receiving  stream  within  relatively  short  dis- 
tances of  each  other,  and  in  some  cases  all  such  wastes 
are  processed  by  municipal  facilities  and  discharged  to 
the  water  body  at  one  location.   Thus,  these  continuous 
waste  products  are  concentrated  within  a  relatively  small 
volume  of  the  receiving  water.   Intermittent  precipitation 
falling  on  urban  areas  becomes  contaminated  as  it  enters 
and  passes  through  or  within  the  man-made  environment. 
The  first  quality  degradation  occurs  when  the  rainwater 
comes  into  contact  with  pollutants  in  the  air.   The 
Cincinnati  Water  Research  Laboratory  of  the  Robert  A. 
Taft  Sanitary  Engineering  Center  started  rainfall  sampling 
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in  August  1963  (Wei  be!  et  al- ,  1966).   Constituents  found 
in  the  rainfall  from  August  1963  to  December  1964  are 
listed  in  Table  2-i.   The  concentrations  represent  average 
values  from  35  storms. 

Next,  the  surface  runoff  passes  over  ground  and 
building  surfaces,  carrying  suspended  sediment  from  ero- 
sion sites  and  dissolving  other  impurities..  Table  2r2 
summarizes  typical  ranges  of  pollutant  concentrations 
found  in  urban  stormwater  runoff  (Field  and  Taf uri ,  1 973) . 
h'nally,  the  stormwater  runoff  comes  into  contact  with 
(1)  solid  residues  deposited  from  earlier  storms  through- 
out the  conveyance  system  and  appurtenances,  and  (2)  dry 
weather  flow  (DWF)  in  combined  sewer  systems.   This  storm, 
runoff  is  well  mixed  with  sanitary  sewage  under  condi- 
tions of  turbulent  flow  in  a  combined  sewer  system,  and 
it  eventually  discharges  to  the  receiving  stream.   Table 
2-3  lists  ranges  of  pollutant  concentrations  found  to  be 
typical  of  combined  sewer  overflows  (Field  and  Taf uri , 
1973).   From  an  examination  of  Tables  2-1,  2-2,  and  2-3, 
it  can  be  verified  that  the  degradation  undergone  by 
urban  stormwater  passing  through  the  surface  runoff  phase 
of  the  hydrologic  cycle  (including  the  collection  system) 
can  be  several  orders  of  magnitude  greater  than  that 
experienced  by  rainwater  during  the  precipitation 
phase. 
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Table  2-1 

Constituents  in  Rainfall 
(Weibel  et  al-  ,  1966) 


Constituent 


Suspended  solids,  mg/l 

COD,  mg/l 

Total  N,  mg/l  (N) 

Inorg.  N,  mg/l  (N) 

Hydrolyzable  PO^, 
mg/l  (PO^) 

Org.  Chlorine*,  mg/l 


Concentration 


13 

16 
1.27 
0.69 

0.24 
0.28 


*From  pesticides  in  the  air. 
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Table  2-2 


Characteristics  of  Urban  Stormwater 
(Field  and  Tafuri  ,  1973) 


Constituent 

Concentration  Range 

BOD5 

;  1 

- 

700  mg/1 

COD 

5 

-■ 

3,100  mg/1 

TSS 

2 

-  . 

11,300  mg/1 

Total  Sol  ids 

450 

- 

14,600  mg/1 

Volatile  Total  Solids 

12 

- 

1  ,600  mg/l 

Settleable  Solids 

0.5 

- 

5,400  mg/1 

-       Organic  N  . 

0.1 

-- 

.   16  mg/1 

- NH3-N 

0.1 

- 

2.5  mg/1 

Soluble  PO^  (as  PO^) 

0.1 

10  mg/1 

Total  PO^  (as  PO^) 

O.T 

- 

125  mg/1 

Chlorides 

2 

- 

25,000  mg/1* 

Oils 

0 

- 

110  mg/1 

Phenols 

0 

- 

0.2  mg/1 

Lead 

0 

- 

1 .9  mg/1 

Total  Coliform 

200 

- 

146x10^/100  ml 

Fecal  Col iform 

55 

- 

112x10^/100  ml 

Fecal  Streptococci 

200 

- 

1.2x10^/100  ml 

*With  highway  deicing 
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Table  2-3 


Characteristics  of  Combined  Sewer  Overflows 
(Field  and  Taf uri ,  1973) 


Constituent 


Concentration  Range 


BOD5 

30 

TSS 

20 

Total  Solids 

150 

Volatile  Total  Solids 

15 

PH 

4.9 

Settleable  Solids 

2 

Organic  N 

1.5 

NHj-N 

0.1 

Soluble  PO^  (as  PO^) 

0.1 

Total  Coliform 

20 

,000 

Fecal  Coliform 

20 

,000 

Fecal  Streptococci 

20 

,000 

600  mg/1 

1 ,700  mg/1 

2,300  mg/1 

820  mg/1 

8.7 

1  ,550  mg/1 

33.1  mg/1 

12.5  mg/1 

6.2  mg/1 

90x10^/100  ml 

17x10^/100  ml 

2x10^/100  ml 
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The  various  constituents  listed  so  far  represent 
a  large  and  diverse  number  of  chemical,  physical,  and 
bacteriological  indicators  of  water  quality.   The  most 
common  parameters  measured  are  listed  in  Table  2-4  under 
their  appropriate  categories.   The  bacteriological  pro- 
cedures are  designed  to  detect  potential  health  hazards 
from  contamination  of  the  water  with  human  or  animal 
feces.   The  sampling  technique,  frequency  of  sampling,  and 
method  of  preservation  should  be  tailored  to  the  indicators 
chosen  for  measurement.   Field  observations  are  extremely 
valuable  for  verification  of  mathematical  models.   Greater 
precision  may  be  obtained  in  the  laboratory,  but  if  the 
sample  is  unrepresentative,  greater  accuracy  will  be 
achieved  by  j_n  si  tu  procedures. 

The  wastewater  constituents  that  affect  the  dis- 
tribution of  dissolved  oxygen  (P.O.)  in  a  natural  water 
system  are  wel 1  documented.   The  oxygen  demand  of  sewage, 
sewage  treatment  plant  effluent,  polluted  stormwater  run- 
off or  industrial  wastes  is  exerted  by  three  types  of 
materials  (American  Public  Health  Association,  1971): 

1.  carbonaceous  organic  matter  oxidized  by 
heterotrophic  bacteria  for  energy  and  cell 
synthesis; 

2.  organic  nitrogen  compounds  hydrolyzed   into 
ammonia -nitrogen  (NH3-N),  then  oxidized  by 
autotrophic  bacteria  (Ni trosomonas  europaea) 
to  nitrite-nitrogen  (NOJ-NTI   further  oxidized 


17 


Table  2-4 
Pollution  and  Contamination  Indices 


Physi  cal  Parameters : 

Temperature 

Turbidity 

Color 


Chemical  Parameters: 

Oxygen  Demand 

Biochemical  Oxygen  Demand  (BOD' 
Chemical  Oxygen  Demand  (COD)  ^ 
Total  Organic  Carbon  (TOC) 

Nitrogen  Compounds 

Organic 
Nitrite 
Nitrate 

Phosphorus  Compounds 

Ortho  Phosphorus 
Poly  Phosphates 

Total  Sol  ids 

Dissolved 
Suspended 

Volatile  and  Fixed 
Settleable 

Chlorides 

Sulfates 

pH 

Alkalinity 

Hardness 


Table  2-4  continued 


Chemical  Parameters  continued 

Heavy  Metals 

Lead 
Copper 
Zinc 

Chromi  urn 
Mercury 

Biological  Parameters: 

Plankton 

Peri  phyton 

Macrophyton 

Macro  in vertebrates 

Fish  Bioassays 


Bacteriological  Parameters 

Total  Coliform  Count 

Fecal  Coliform 

Fecal  Streptococci 

Total  Plate  Count 
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by  Ni trobacter . wi nogradskyi  to  ni trate-n1 trogen 
(NO3-N};  and 

3.  certain  chemical  reducing  compounds  (ferrous 
iron,  sulfite,  and  sulfide)  which  will  react 
with  molecularly  dissolved  oxygen.    , 

This  oxygen  demand  is  the  response  of  aquatic  biota 
to  an  adequate  food  supply  and  is  commonly  referred  to  as 
the  biochemical  oxygen  demand  (BOD).   The  laboratory  BOD 
technique  is  an  empirical  bioassay-type  procedure:   the 
D.O.  consumed  by  microbial  life  in  an  incubated  bottle  is 
measured  with  respect  to  time  at  a  specified  temperature. 
The  actual  environmental  conditions  of  temperature  changes, 
biological  population,  water  movement,  sunlight,-  and  zones 
of  aerobic  and  anaerobic  processes  cannot  be  faithfully 
reproduced  in  the  laboratory.  .  Thus,  the  "bottled"  system, 
on  a  kinetic  comparison,  is  completely  accurate  in  repre- 
senting itself  but  may  be  relatively  unreliable  as  a 
representation  of  the  source  from  which  the  sample  was 
taken.   The  basic  assumption  that  consumption  of  .D.O.  is  an 
absolute  and  complete  parameter  of  biological  decomposition 
in  the  BOD  bottle  constitutes  a  simplification  of  complex 
interactions  (Liptak,  1974). 

Other  laboratory  methods  have  been  developed  to 
measure  the  oxygen  demand  exerted  by  organic  matter.   The 
chemical  oxygen  demand  (COD)  and  total  organic  carbon 
(TOC)  tests  are  more  precise  chemical  methods,  but  the 
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analytical  results  are  not  accurate  if  the  organic  ma- 
terial measured  is  not  equivalent  to  the  organic  matter 
actually  being  utilized  by  the  microorganisms  in  the 
stream.   Much  can  be  done  to  improve  the  accuracy  of  the 
BOD  test  by  using  dilution  water  from  the  receiving  stream, 
thus  introducing  a  natural  "seed"  of  diverse  organisms 
into  the  bottle  system.   With  all  of  its  limitations,  the 
BOD  procedure  is  still  considered  to  be  the  best  method 
for  evaluating  the  effect  of  waste  inputs  on  the  oxygen 
balance  of  a  stream  (Nemerow,  1974). 

2.2  A  Unifying  Conceptualization  of  the  Physical  System 

The  principle  of  conservation  of  mass  may  be 
applied  universally  throughout  the  urban  and  natural 
environments  to  describe  the  transport  of  pollutants. 
Figure  2-1  represents  a  generalized  component  of  the 
physical  system,  which  may  characterize  (1)  a  pipe 
segment  of  the  sewer  system,  (2)  a  storage/treatment 
device,  or  (3)  a  reach  of  the  receiving  body  of  water. 
Essentially,  all  of  these  subsystems  may  be  approximated 
by  the  one-dtmensional  version  of  the  convecti ve-dispersion 
equation. 


—  =  —  TE—  -  UCl  +  FS 
3t     8x   '-''3X   .  ^-I  -  ^^ 


(2.1) 
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Figure  2-1.   A  Generalized  Component  of  the  Physical 
System. 
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where   C  =  concentration  of  water  quality  parameter 
(pollutant),  M/l3, 

t  =  time,  T, 

-E—  =  mass  flux  due  to  longitudinal  dispersion 
^^   along  the  flow  axis,  the  x  direction, 
M/l2t, 

UC  =  mass  flux  due  to  -advection  by  the  fluid  con 
taining  the  mass  of  pollutant,  M/L^T, 

S  =  sources  or  sinks  of  the  substance  C, 
M/l3t, 


U  =  flow  velocity,  L/T,  and 


E  =  longitudinal  dispersion  coefficient,  L  /I. 


The  source/sink  term  accounts  for  biochemical  processes 
(e.g.,  decay,  photosynthesis,  algal  respi ration) ,  bound- 
ary losses  such  as  stream  benthic  deposits,  and  boundary 
gains  (e.g.,  reaeration,  and  point  or  distributed  waste 
discharges).   Assuming  that  the  longitudinal  dispersion 
coefficient  and  the  advective  velocity  are  constant  along 
the  flow  axis,  equation  (2.1)  may  be  expanded  for  the 
generalized  component  in  Figure  2-1  to 


3C  _  r  i^ 


U  |£ 

3x 


-  KC  +  ^  (C^ 


C) 


(2.2) 


where   K  =  first-order  reaction  rate  coefficient, 
1/T, 
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q^    =  influent  fluid  flow  rate  per  unit  width, 
^    L2/T, 

C-j  =  concentration  of  water  quality  parameter 
in  the  inflow,  M/L-^,  and 


A  =  wetted  cross-sectional  area  in  the  com- 
ponent, l2. 


In  equation  (2.2)  the  influent  fluid  flow  rate  per  unit 
width,  influent  concentration,  and  the  wetted  cross- 
sectional  area  may  all  be  variable  functions  of  time. 

For  a  nondi spersi ve  pipe  segment,  equation  (2.2) 
may  be  reduced  to  the  form  for  simple  advection. 


Il  ^   ^If  ^   ^^ 


(2.3) 


The  solution  to  equation  (2.3)  is  given  by 


C(x,t)  =  C^(x-Ut,  t  -  ^)exp[-Kt] 


(2.4) 


where   C,(x-Ut,  t  -  tj-)  =  initial  concentration  at  a 

time  x/U  earlier  and  a  dis- 
tance Ut  upstream,  M/L^. 


Equation  (2.4)  states  that  the  reactive,  nondi spersi ve 
advective  system  has  an  instantaneous  response  which 
travels  downstream  at  a  velocity  U,  and  the  magnitude 
of  the  response  decreases  exponentially  because  of 
first-order  decay.   Chemical  engineers  refer  to  a  system 
governed  by  equation  (2.3)  as  a  plug  flow  reactor.   If 
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there  is  no  decay,  or  if  it  can  be  assumed  that  it  is 
negligible  because  the  travel  time  in  the  system  is  very 
short  (a  large  U),  then  equation  (2.4)  reduces  to 


C(x,t)  =  C^(x  -  Ut.  t  -  ^) 


(2.5) 


This  is  a  reasonable  assumption  to  make  for  most  sewer 
systems  when  the  study  objectives  are  not  design-oriented. 

Extensive  studies  have  been  conducted  by  chemical 
engineers  concerning  the  kinetic  analysis,  flow  charac- 
teristics, and  dispersion  characteristics  of  chemical 
process  reactors.   Danckwerts  (1953)  studied  the  residence- 
time  distribution  of  continuous-flow  systems.   The  basic, 
differential  equation  for  a  continuous-flow  reactor  is 
gi  ven  as 


2 

8  C 


U  ii 
3x 


KC 


(2.6) 


for  a  first-order  reaction.   Equation  (2.6)  is  obtained 
from  equation  (2.2)  for  steady-state, conditions  (8C/3t  = 
0).   Wehner  and  Wilhelm  (1956)  solved  equation  (2.6) 
analytically.   Their  solution  has  been  recommended  by 
Thirumurthi  (1969)  for  the  design  of  large  waste  stabili 
zation  ponds,  aerated  lagoons,  and  activated  sludge  aera 
tion  tanks . 

For  completely  mixed  storage/treatment  systems 
(3C/3X  =  0),  equation  (2.1)  reduces  to 
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—   =   E  S 


(2.7) 


which  may  be  rewritten,  for  the  sources  and  sinks  of  the 
generalized  component  in  Figure  2-1,  as  the  ordinary  dif- 
ferential equation 


v^.  =    ^1^1  -  V 


KC  V 


(2.8) 


where   Qi  =  inf 1 uent  f 1 uid  flow  rate,  L  /T, 


V  = 


outflow  rate,  L  /T,  and 

volume  of  the  fluid  mass  in  the  component, 
L3. 


Solutions  to  equation  (2.8)  are  presented  in  Chapter  V 
for  the  constant  volume  and  variable  vol ume  cases ,  for 
time-varying  influent  flow  rates  and  concentrations.   In 
chemical  engineering  texts  (e.g..  Cooper  and  Jeffreys, 
1971),  the  design  equation  for  a  continuous-flow  stirred 
tank  reactor  is  given  as 


1  +  Kt, 


(2.9) 


where   C  =  constant  influent  concentration,  M/L  , 
°   and 


t„  =  detention  time  in  the  reactor,  T, 
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This  equation  is  used  by  Eckenfelder  (1966)  in  designing 
aerated  lagoons  for  industrial  water  poll ution  control . 
It  should  be  noted  that  equation  (2.9)  predicts  the 
steady-state'  value  of  the  general  solution  to  equation 
(2.8)  for  a  constant  volume  system.   A  methodology  is 
also  presented  in  Chapter  V  to  obtain  solutions  to  equa- 
tion (2.1)  for  the  more  difficult  case  of. a  dispersive 
variable  volume  system. 

A  receiving  body  of  water,  such  as  a  river  reach 
or  an  estuary,  is  goverened  by  equation  (2.1).   Whereas 
C  represented  the  BOD  concentration  of  wastewaters  in 
the  urban  environment,  equation  (2.1)  is  usually  solved 
to  obtain  the  dissolved  oxygen  concentration  in  the 
natural  system.   The  BOD  inputs  from  the  urban  area  be- 
come oxygen  sinks,  along  with  benthal  demand,  respira- 
tion by  aquatic  plants,  and  others.   An  oxygen  concen- 
tration gradient  is  developed  between  the  atmosphere 
and  the  waterway,  leading  to  reaeration  as  a  dissolved 
oxygen  source.   This  process  is  supplemented  by  algal 
photosynthesis.   A  steady-state  analysis  of  equation 
(2.1)  is  presented  in  Chapter  VII. 

As  noted  in  Figure  2-1 ,  all  of  the  systems 
described  by  the  generalized  component  are  essentially 
characterized  by 

(1 )  a  flow  velocity  U, 

(2)  a  reaction  rate  coefficient  K, 
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(3)  a  residence  time  tp.,  and 

(4)  a  longitudinal  dispersion  coefficient  E, 
all  of  which  are  a  function  of  the  system  geometry  (vol- 
ume V,  cross-sectional  area  A,  length  x).   Velocities  in 
the  sewer  systems  tend  to  be  much  higher  than  in  either 
the  storage/treatment  devices  or  the  receiving  waters. 

In  a  receiving  stream  reach,  lake  or  bay  the  detention 
time  is  much  longer  than  in  the  other  systems  due  to  the 
large  volume  of  the  fluid  mass  relative  to  the  wastewater 
inflows.   From  performance-time  curves  of  activated  sludge 
units,  an  overall  BOD  first-order  decay  coefficient  of 
0.2  per  hour  (4. 8/day)  may  be  obtained  for  the  biologi-_ 
cal  waste  treatment  process  (Fair,  Geyer  and  Okun,  1968). 
For  tidal  rivers  and  estuaries  a  first-order  decay  rate 
for  BOD  not  much  higher  than  0.02  per  hour  (0.5/day) 
can  be  expected  (Hydroscience ,  1971).   In  lakes,  K  may 
vary  from  0.01/day  for  the  dry  season  to  0.02/day  for 
the  wet  season  (Huber  e^t  aj_.  ,  1976),  and  the  detention 
time  is  on  the  order  of  months.   By  contrast,  a  deten- 
tion time  of  10  to  11  hours  is  common  design  practice 
for  conventional  activated  sludge  plants  (American 
Society  of  Civil  Engineers,  1959).   Estimates  of  K,  t^, 
and  the  product  Ktj,  are  summarized  for  various  systems 
in  Table  2-5.   The  reaction  rate  group  Kt„  provides  a 
rough  idea  of  treatment  potential. 
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Table  2-5 


Treatment  Potential  of  Various  Components 
of  the  Physical  System 


Component 


Selected  Decay        Selected  Detention 
Coefficient  Time,  t^  Kt^ 

K  at  20° 


Pipe  Segment^ 
of  Sewer  System 

0.005/hour 

0.56  hours 

0.003       ; 

Activated-Sluge 
Plant  for  Dry 
Weather  Flow 

0.2/hr 

11  hours 

2.20 

Waste  Stabili- 

ziation  Pond° 
for  Raw  Sewage 

0.061 /day 

39.8  days 

2.40 

Wet  Weather 

Flow  Storage/ 

Treatment 

System 

0.005/hour 

24  hours 

0.12 

Kissimmee 
River^ 

0.3/day 

9  days 

2.70 

Lake  Toho 

0.015/day 

5  months 
(150  days) 

2.25 

Tidal  River 
or  Estuary 

0.5/day 

5  to  50  days 

2.5  to  25.0 

^Segment    length    10,000   feet    (3    km)    with   mean    velocity   of 
5   ft/sec    (1 .5  m/sec) . 

''Average   of   five   ponds    in   Fayette,   Missouri    (Thirumurthi , 
1969). 

*^Kissimmee    River,    Florida    for   a    length    of   90    (144    km) 
from     Huber   ei  al- »    1976;    and   Hydroscience,    1971. 

^Lake   Toho,    Florida,    decay   rate    for    total    phosphorus, 
from   Huber   et   al . ,    1976. 
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Comparisons  based  on  the  reaction  rate  group 
Ktp  in  a  dispersive  system  must  be  qualified.   For  ex- 
ample, the  longitudinal  dispersion  coefficient  in  the 

aeration  tank  of  an  activated  sludge  process  may  be  ex- 

2 
pected  to  be  around  1.6  ft  /sec  (Murphy  and  Timpany, 

Z 
1967).   For  a  tidal  river  or  estuary,  about  300  ft, /sec 

represents  the  lower  limit  of  the  longitudinal  dispersion 
coefficient.   The  differences  in  mass  flux  due  to  advec- 
tion  (U)  and  dispersion  (E)  will  vary  the  system  response 
considerably.   Additional  discussion  of  the  importance 
of  these  parameters  is  presented  in  Chapter  V.   Further- 
more, bays  and  lakes  may  be  treated  in  a  detailed  manner 
by  a  multi-dimensional  form  of  equation  (2. 1 ),  namely , 
a  2-dimensional  or  3-dimensional  expression  of  the  con- 
vective  diffusion  equation.   Crude  approximations  may 
be  obtained  by  simply  applying  equation  (2.8),  which 
assumes  completely  mixed  systems. 


CHAPTER  III 

DESCRIPTION  OF  THE  STUDY  AREA  AND 
NECESSARY  MODELING  TECHNIQUES 


3.1   The  Need  for  a  Study  Area 

Pitt  and  Field  (1974)  applied  available  informa- 
tion on  urban  runoff  pollution  to  a  hypothetical  test  area, 
subjected  to  a  worst-case  storm  event,  to  demonstrate 
potential  problems  and  their  solutions.   However,  it  has 
al ready  been  establ ished  that  there  is  a  need  for  contin- 
uous hydrologic  simulation  to  assess  the  frequency  with 
which  runoff  events  cause  adverse  effects  in  the  receiving 
waters.   It  would  be  quite  difficult  to  generate  a  realis- 
tic pollutant  distribution  for  a  long  sequence  of  synthetic 
flows.   There  are  other  important  reasons  that  justify 
selection  of  a  real  study  area.   The  only  way  to  estab- 
lish the  necessary  validity  which  renders  a  mathematical 
model  of  water  quality  useful  for  planning  purposes  is 
to  conduct  verification  procedures  and  calibrate  against 
field-measured  data  (Hydroscience ,  1971).   Furthermore, 
when  evaluating  the  effectiveness  of  proposed  control  mea- 
sures it  is  important  to  base  comparisons  against  existing 
conditions  in  the  study  area. 
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3.2   Des  Moines,  Iowa 


The  city  of  Des  Moines,  Iowa  is  located  at  the 
confluence  of  the  Des  Moines  River  and  the  Raccoon  River 
as  shown  in  Figure  3-1.   Within  its  limits  are  approximately 
200,000  people  out  of  a  total  of  288,000  for  the  metro- 
politan area.   The  mean  annual  precipitation  is  about  31 
inches  (787  mm),  approximately  equal  to  the  State  of  Iowa 
average  and  the  United  States  average  (Waite,  1974). 
The  urban  area  covers  49,000  acres  (19,830  ha)  of  land 
which  has  gently  rolling  terrain.   Most  of  the  area,   . 
45,000  acres,  is  served  by  separate  sewers,  while  4,000 
acres  are  served  by  combined  sewers. 

Selection  of  the  study  area  was  based  primarily 
on  data  availability.   Davis  and  Borchardt,  of  Henningson, 
Durham  &  Richardson,  Inc.  ,  Omaha,  Nebraska,  conducted  an 
extensive  sampling  program  of  combined  sewer  overflows, 
stormwater  discharges,  and  surface  waters  in  the  Des  Moines, 
Iowa  Metropolitan  Area  for  the  U.  S.  Environmental  Pro- 
tection Agency.   The  objective  was  a  combined  sewer 
overflow  abatement  plan  for  the  metropolitan  area.   The 
sampling  program  was  conducted  from  March  1968  to  October 
1969.   Other  considerations  revolved  around  the  fact  that 
Des  Moines,  Iowa  is  somewhat  typical  of  many  urban  centers 
throughout  the  country. 

(1)   it  has  a  medium-sized  population; 


VICIMITY       MAP 


Figure  3-1.   Map  of  Des  Moines  Area  (Davis  and  Borchardt 
1974). 
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(2)  its  domestic  and  industrial  dry  weather 
flows  receive  secondary  treatment; 

(3)  its  wastewaters  are  discharged  into  a  non- 
tidal  receiving  stream;  and 

(4)  the  urban  area  receives  a  mean  annual  pre- 
cipitation equal  to  the  national  average. 

The  mean  annual  temperature  is  50°C  (10°C),  vary- 
ing from  an  average  of  21°F  (-6°C)  in  January  to  an  average 
of  75°F  (24°C)  in  July. 


3.3 


Pollutant  Loads 


Annual  pollutant  unit  loads  upstream  from  the  city 
are  summarized  in  Table  3-1 ,  as  determi ned  by  Davis  and 
Bprchardt  (1974).   From  data  gathered  by  wastewater  treat- 
ment plant  personnel  during  the  period  August  to  October, 
1969,  the  average  daily  flow  was  estimated  at  35.3  mgd 
(133,600  cu  m/day)  with  an  organic  load  from  domestic 
and  industrial  sources  of  85,800  pounds  of  BOD  (43,450  kg 
BOD  per  day).   The  high-rate  trickling  filter  facility 
achieved  84  to  85  percent  BOD. removal,  producing  an 
effluent  BOD  concentration  of  49  mg/1  (Office  of  Water 
Planning  and  Standards,  1974).   A  simplified  profile  of 
the  DWF  treatment  facility  is  shown  in  Figure  3-2.   The 
estimated  annual  loading  from  the  urban  area's  separate 
and  combined  sewer  systems  is  presented  in  Table  3-2. 
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Table  3-1 

Pollutant  Unit  Loads  for  Drainage  Area* 
Above  Des  Moines,  Iowa 
(Davis  and  Borchardt,  1974) 


Des  Moines 
Ri  ver 


Raccoon 
River 


Total 


Drainage  Area,  acres 
(ha) 


3,738,000     2,202,000   5,940,000 
(1,512,769)    (891,149)    (2,403,918) 


Unit  Average  Annual 

Runoff,  acre-ft/acre  0.42  0.40  0.41 

(ha-m/ha)  (0.13)  (0.12)  (0.12) 

Unit  BOD,  lbs/acre  13.40  6.93  11.01 

(kg/ha)  (15.02)  (7.77)  (12.34) 

Unit  NOo,  lbs/acre  3.75  3.74  3.75 

(kg/ha)  (4.20)  (4.19)  (4.20) 

Unit  PO/i,  lbs/acre  0.54  0.42  0.50 

(kg/ha)^  (0.61)  (0.47)  (0.56) 


*0n  an  annual  basis 
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Taking  the  total  upstream  drainage  area  for  the 
Raccoon  and  Des  Moines  Rivers,  the  pollutant  contributions 
are  65,225,000  pounds  of  BOD  (29,586,000  kg);   22,222,000 
pounds  of  NO3  (10,080,000  kg);  and  2,940,000  pounds  of 
PO.  (1,334,000  kg).   The  urban  area  loadings  (when  added 
to  upstream  values)  represent,  respectively,  15  percent, 
3  percent  and  51  percent  of  the  total  BOD,  NO3,  and  PO^ 
mass  loadings  to  the  Des  Moines  River  below  the  metropoli- 
tan area.   The  Davis  and  Borchardt  estimates,  made  from 
river  sampling  data  taken  below  Des  Moines,  indicate 
the  following  average  annual  river  loadings:   70,000,000 
pounds  of  BOD  (31,851,466  kg);  25,400,000  pounds  of  NO3 
(11,521,250  kg),  and  7,950,000  pounds  of  PO^  (3,606,059 
kg).   These  figures  reveal  that  (1)  6,610,000  pounds  of 
BOD  (2,998,246  kg)  are  "lost"  in  transit  through  the  urban 
section  of  the  stream,  and  (2)  by  contrast  -2,474,360 
pounds  of  NO3  (1,122,351  kg)  and  1,917,400  pounds  of  PO^ 
(869,718  kg)  are  gained  in  addition  to  the  measured  urban 
sources. 

The  authors  of  the  report  offer  some  explanations. 

The  'sometimes'  decreases  in  organic 
load  through  the  metro  area  may  be  attribu- 
table to  treatment  realized  in  the  low 
head  impoundments  at  Scott  and  Center 
Streets  on  the  Des  Moines  River  and  just 
below  Fleur  Drive  on  the  Raccoon.   To  some 
extent  these  impoundments  may  be  serving _ 
as  intermittent  sedimentation  and  stabili- 
zation units. 

All  BOD  data,  including  those  used 
from  the  two  other  studies,  were  obtained 
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from  unfiltered  samples.   However,  since 
the  analytical  technique  was  the  same  for 
all  samples,  the  relative  magnitude  of  the 
data  should  not  be  affected. 

There  has  been  some  speculation  that 
treated  wastewater  effluents  may  exert  an 
antagonistic  or  retardant  effect  on  the 
BOD  exertion  rate  of  the  receiving  stream. 
If  true,  this  may  be  due  to  surfactants  or 
to  the  expected  lower  exertion  rate  of  the 
effluent.   In  this  regard,  the  decreased 
BOD  in  4  or  5  measurements  between  R-5 
and  R-6  is  of  interest.   Increased  loads 
between  the  summation  of  R-4  and  R-9  ver- 
sus R-5  are  likely  due  to  raw  and  combined 
sewage  bypassing  the  intervening  area. 

Another,  and  probably  the  most  prac- 
tical, possibility  for  the  discrepancies 
is  the  fact  that  the  data  are  biological 
and  biochemical  in  nature  and  such  data 
do  not  always  provide  predictable  compara- 
tive summations,  (p.  108) 

The  river  sampling  stations  {R-4,  R-5,  R-6  and 
R-9)  are  shown  in  Figure  3-3.   It  is  possible  that  inter- 
vening creeks,  such  as  Beaver  Creek  which  carries  nutrient 
loads  of  2,860,000  pounds  of  NO3  per  year  (1,297,274  kg 
per  year)  and  390,000  pounds  of  PO^  per  year  (176,900  kg 
per  year),  may  answer  observed  differences  in  the  nitrate 
loads.   However,  the  phosphate  totals  remain  unbalanced 
and  the  cause  unresolved. 

Combined  sewer  overflows  were  monitored  at  five 
locations  and  stormwater  discharges  at  three  locations. 
The  average  BOD  concentration  of  stormwater  was  reported 
as  53  mg/1 ,  while  the  average  BOD  concentration  in  com- 
bined sewer  overflows  was  72  mg/1. 
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Figure    3-3.      Location   Map:      River    Sampling    Points    (Davis 
and   Borchardt,    1974)  . 
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3.4  The  Receiving  Stream 


From  Figure  3-3,  it  is  clear  that  although  the 
Raccoon  River  receives  urban  stormwater  discharges  from 
parts  of  the  urban  area,  the  ultimate  recipient  is  the 
Des  Moines  River.   All  stormwater  discharges,  combined 
sewer  overflows,  and  wastewater  treatment  plant  effluent 
are  assumed  to  enter  the  Des  Moines  River  at  its  con- 
fluence with  the  Raccoon  River.   The  Des  Moines  River 
stretches  for  200  miles  (322  km)  from  the  City  of  Des 
Moines  to  its  junction  with  the  Mississippi  River  andj  in 
general,  the  river  is  wide  and  swift  with  occasional  deep 
holes  and  a  broad  flood  plaint'  According  to  thef/State  ;  ;  ,■ 
Hygienic  Laboratory  (1974)  bottom  material  is  composed  of 
silt  deposits,  sand,  gravel  and  rubble  providing  numerous  ^ 
habitats  for  fish  and  other  aquatic  life.   Recreational 
activities  such  as  fishing  and  boating  are  quite  heavy. 
The  entire  reach  is  classified  as  warm  water  "B"  stream    . 
by  the  Iowa  Water  Quality  Standards,  such  that  the  absolute 
minimum  dissolved  oxygen  level  specified  is  4.0  mg/T.  ■  The 
Iowa  Standards  also  require  a  minimum  of  5.0  mg/1  during 
at  least  16  hours  per  day  (State  Hygienic  Laboratory, 
1970). 

For  the  heavy  precipitation  months  of  June,  July 
and  August,  1968^  Davis  and  Borchardt  reported  the  follow- 
ing nutrient  concentrations  at  a  point  approximately  5.5 
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miles  (9.0  km)  downstream  from  the  confluence  of  the  Raccoon 
and  Des  Moines  Rivers: 

(1)  total  organic  nitrogen  ranged  from  1.5 
to  3.7  mg/1 , 

(2)  nitrate  nitrogen  ranged.from  0.2  to  7.8 
mg/1,  averaging  3  to  4  mg/1 ,  and 

(3)  orthophosphate  (0-PO^)  ranged  from  0.6  to 
1.8  mg/1,  averaging  si ightly  over  1 . 0  mg/1. 

Davis  and  Borchardt  also  observed  high  algal  densities  in 
both  the  Des  Moines  and  Raccoon  Rivers,  and  state  that 
nutrient  concentrations  are  almost  always  present  at  levels  re- 
ported by  Sawyer(T966)  to  be  suf  f  i  cient  for  nui. sance  al  gal    ' 
growths:   0.3  mg/1  for  inorganic  nitrogen  (NH^,  NO2,  NO3) 
and  0.015  mg/1  of  inorganic  phosphorus. 

3.5  An  Abstraction  of  the  Physical  System 

To  study  the  interaction  of  urban  runoff, .control 
measures,  and  receiving  water  response,  an  analytic  frame- 
work is  sought  that  describes  the  key  elements  of  the  physi- 
cal system.   Conceptual  models  are  normally  expressed  in 
mathematical  terms  (Vollenweider,  1975).   A  conceptual 
model  is  derived  from  an  abstraction  of  the  physical  system, 
and  such  an  image  is  presented  in  Figure  3-4. 

The  urban  community  served  by  a  separate  sewer 
system  will  convey  stormwater  runoff  and  municipal  sewage 
through  conduits  which  are  not  connected  together.   The 
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Figure  3-4.   An  Abstraction  of  the  Physical  System, 


44 


BOD  concentration  of  the  storm  sewer  runoff  is  mixed  with 
the  dry  weather  flow  (DWF)  and  accumulated  sewer  solids. 
An  interceptor  carries  the  sanitary  design  f 1 ow  to  the 
municipal  sewage  treatment  plant.   Since  complete  mixing 
is  assumed,  the  BOD  concentrations  of  the  combined  sewer 
overflow  and  the  flow  intercepted  for  treatment  by  the  DWF 
plant  are  identical.   The  combined  sewer  overflow  and  the 
separate  storm  sewer  flow  are  both  processed  through  a 
central  wet  weather  flow  (WWF)  storage/treatment  facility 
before  discharge  to  the  receiving  body  of  water. 

An  examination  of  Figure  3-4  from  top  to  bottom 
reveals  three  distinct  phases  of  the  water  quality  cycle  ,. 
within  the  urban  environment.   At  the  top  level  is  the 
pollutant  build-up  phase,  represented  by  surface  runoff 
generated  within  the  separate  and  combined  sewer  areas  as 
well  as  industrial  and  domestic  wastewater.   Estimation 
of  DWF  sanitary  loads  is  a  relatively  simple  task.   The 
design  flow  for  a  sanitary  wastewater  treatment  plant  is 
based  on  water  consumption  within  the  urban  area,  and  the 
organic  waste  load  is  a  function  of  population  density  and 
data  obtained  from  industrial  and  commercial  establishments 
However,  the  estimation  of  surface  runoff  quantity  and 
qual i  ty  is  more  compl ex. 

Urban  runoff  is  a  direct  result  of  precipitation 
over  an  urban  area.  However,  due  to  temporal  and  spatial 
variabilities  associated  with  atmospheric  motions,  all 
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hydrologic  processes  are  more  or  less  stochastic  (Chow, 
1964).   If  a  deterministic  system  receives  a  random  input, 
output  from  the  system  will  also  be  random.   A  single  time 
history  representing  a  random  phenomenon  is  referred  to  as 
a  sample  record  when  observed  over  a  finite  time  interval 
(Bendat  and  Piersol  ,  1971).-  By  selecting  a  sample  record 
of  point  precipitation,  a  random  element  in  the  hydrologic 
process  is  eliminated:,  chance  of  occurrence.  The  precipi- 
tation time  series  may  then  be  fed  into  an  urban  hydrologic 
simulation  model  to  generate  the  urban  runoff  series.  The 
spacing  and  sizing  of  individual  events  in  the  sequence, 
however,  is  probabi 1 i stic  (Eagleson,  1970) .   Therefore, 
a  combination  of  probabilistic  and  deterministic  methods 
must  be  used.   At  this  stage,  analysis  of  time  series,  a 
classical  statistical  technique,  is  applied  to  define  a 
minimum  interevent  time  and  group  runoff  occurrences  into 
Independent  storm  events.   A  detailed  presentation  is  made 
in  Chapter  IV. 

The  second  phase  in  the  urban  water  quality  cycle 
is  at  the  middle  level  In  Figure  3-4,  represented  by  the 
man-made  conveyance  and  pollutant  control  systems.   Since 
travel  times  i n  the  sewer  system  are  usually  relatively 
short,  it  has  been  assumed  that  no  pollutant  decay  occurs 
in  the  conduits.   The  BOD  concentration  of  the  DWF  treat- 
ment plant  effluent  is  given  by 
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BOD 


[BOD^  •  DWFSEP  +  BOD^  •  DWFCMB]  (l-R^,) 


DWFSEP  +  DWFCMB 


(3.1) 


where   BOD^  =  BOD  concentration  of  municipal  sewage, 
^   mg/1  . 

BOD  =  mixed  BOD  concentration  in  the  combined 
^   sewer,  mg/1 , 

DWFSEP  =  DWF  contribution  from  separate  sewer 
area,  cf s , 

DWFCMB  =  DWF  contribution  from  combined  sewer 
area ,  cf s ,  and 

R.  =  fraction  removal  of  BOD  achieved  by  the 
.   DWF  treatment  facility. 

Treatment  of  the  DWF  facility  is  accomplished 
in  a  single  time  step.   That  is,  a  mass  of  poll utantj. 
entering  the  DWF  plant  at  a  particular  time  step  is  treated 
and  released  during  that  time  step,  with  no  effects  extend- 
ing beyond  the  discrete  time  interval .   The  effects  of 
storage/treatment  in  the  WWF  facility  are  dynamic  and 
extend  beyond  the  time  step  length,  depending  on  the  de- 
tention time  primarily.   The  BOD  concentration  of  the 
input  to  the  WWF  unit  is  given  by 


BOD. 


BOD 


Qe  +  BOD^ 
^s      c 


sc 


Q  +  Q 
^s    c 


(3.2) 


where   BOD   =  mixed  BOD  concentration  of  urban  runoff 
^-   from  the  separate  storm  sewer  and  com- 
bined sewer  fl ows ,  mg/1 , 
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BOD  =  BOD  concentration  of  urban  stormwater 
^   runoff ,  mg/1 , 

Q  =  urban  runoff  carried  by  the  separate 
^   storm  sewer,  cfs,  and 

Q  =  combined  sewer  flow  below  the  inter- 
^   ceptor,  and  to  the  WWF  storage/treatment 
uni ts  ,  cf s . 

To  characterize  the  response  of  the  WWF  storage/ 
treatment  system,  mass  transport  models  derived  from  the 
concepts  of  continuity  presented  in  Chapter  II  are  applied 
to  the  wet  weather  quantity  and  quality  inputs.   The 
complete  mathematical  analysis  is  presented  in  Chapter  V. 
Frequency  analyses  are  performed  on  concentration  and 
mass  rate  inputs,  and  outputs  from  selected  models  of  the 
storage/treatment  system. 

At  the  bottom  level  of  Figure  3-4  is  the  third 
phase  of  the  water  quality  cycle  in  the  urban  environment: 
the  natural  pollutant  control  system  provided  by  the  waste 
assimilation  capacity  of  the  receiving  body  of  water.   An 
analogy  has  been  drawn  in  Chapter  II  between  the  man-made 
storage/treatment  system  discussed  above  and  the  receiving 
water.   Both  systems  are  governed  by  the  same  principles 
of  conservation  of  mass.   The  mass  of  pollutant  entering 
either  body  is  subjected  to 

(1)  mixing  with  the  impurities  already  present 
in  the  system  within  a  specified  volume; 

(2)  first-order  decay  of  nonconservati ve  sub- 
stances ; 
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(3)  mass  flux  due  to  longitudinal  dispersion 
along  the  flow  axis;  and 

(4)  mass  flux  due  to  advection  along  the  flow 
axis  by  the  fluid  containing  the  pollutant 

The  concentration  of  the  combined  BOD  inputs  in  the  re- 
ceiving water  is  given  by 


BOD, 


BOD^Q^  .  BOD^Q^  +  BOD^Q^ 
Qu  -^  ^d  ^  ^w 


(3.3) 


where   BOD  =  mixed  BOD  concentration  in  receiving  water, 
^        mgAl, 

BOD  =  mixed  BOD  concentration  from  sources  up- 
"   streamofurbanarea,mg/l, 

BOD.  =  BOD  concentration  of  dry  weather  flow 
treatment  plant  effluent,  mg/1 , 

BOD  =  BOD  concentration  of  WWF  storage/treatment 
^   facility  effluent,  mg/T, 

Q  —upstream  flow,  cfs, 

Q  .   =    DWF   treated   effluent,    cfs,    and 

Q  =  WWF  storage/treatment  effluent,  cfs. 
^w  ^^  ' 

The  initial  conditions  of  BOD  in  the  river  are 
defined  by  equation  (3.3),  and  the  hypothetical  impact  on 
the  oxygen  balance  of  the  receiving  stream  is  estimated 
using  simplified  mathematical  modeling  approaches.   Criti- 
cal D.O.  deficits  are  computed  as  a  function  of  several 
stream  parameters:   temperature,  flow,  oxygen  concentration. 
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deoxygenation  and  reaeration  rates,  BOD  concentrations, 
velocity,  and  the  longitudinal  dispersion  coefficient. 
Minimum  DO's  are  calculated  subsequently  and  frequency 
analyses  are  performed. 


3.6  Technique  for  Calculation  of  Urban  Runoff  Quantity 
and  Quality 


This  section  briefly  describes  the  methods  used 
to  generate  storm  runoff  and  pollutant  concentrations. 
The  Hydrologic  Engineeri ng  Center  model  ,  STORM,  is 
utilized  to  obtain  hydrographs  and  pollutographs  for 
Des  Moines  for  the  year  1958  on  an  hourly  time  s.tep. 

Urban  Runoff  Quantity 

STORM  computes  urban  runoff  as  a  function  of  land 
use  and  rai nfal 1/snowmel t  losses  (Hydrologic  Engineering 
Center,  1975) .  " 


.*^  =  "u(^-■ru)  ,:„:.:  ^'■'^ 


where   AR  =  urban  area  runoff,  in/hr, 

CR  =  composite  runoff  coefficient  dependent 

"  on  urban  land  use, 

P  =  hourly  rainfal 1/snowmel t  in  inches  over 

"  the  urban  area,  and 
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f  =  available  urban  depression  storage, 
u    - 
in . 


A  maximum  depression  storage  of  a  hundredth  of 
an  inch  (0.25  mm)  is  assumed  for  Des  Moines,  Iowa.   The 
hourly  urban  runoff  values,  expressed  in  cfs,  are  saved 
in  a  file  for  later  recall  by  the  simplified  mathematical 
model . 

Urban  Runoff  Quality 

The  basic  water  qual ity  parameters  modeled  by 
STORM  are  suspended  and  settleable  solids,  BOD,  total 
nitrogen  (N ),  and  total  phosphate  (PO^).   It  is  important 
to  emphasize  that  the  BOD  values  are  expressed  in  terms  of 


the  standard  BOD^  test:   incubation  at  a  temperature  of 
20°C  for  five  days.   These  values  represent  most  of  the 
carbonaceous  oxygen  demand  exerted  by  organic  matter 
present  in  the  urban  runoff,  and  include  the  BOD  contri- 
bution from  suspended  and  settleable  solids.   The  BOD 
loading  rates  generated  by  STORM  are  based  on  land  use 
and  other  factors  such  as  number  of  dry  days  without 
runoff  since  the  last  storm  and  the  street  sweeping  intervals 

3.7   Separate-,  Combined  and  Dry  Weather  Loading  . 

All  of  the  following  methodology  can  be  used, 
regardless  of  the  technique  employed  to  generate  storm 
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runoff  and  quality,  as  long  as  these  values  pertain  to  the 
entire  area  being  modeled. 

Separate  Flows  and  Loadings 

Apportionment  of  the  total  flow  and  BOD  loading  is 
made  on  the  basis  of  the  relative  area  served  by  separate 
and  combined  sewers.   Runoff  from  separate  sewered  areas 
is  thus  (refer  to  Figure  3-4) 


A 


where   Qg  =  stormwater  flows  from  separate  sewered  areas, 
cfs , 


A  =  area  served  by  separate  sewers,  acres. 


Q^  =  total  (storm  plus  combined)  urban  runoff, 
cfs ,  and 


A^  =  total  area  of  catchment,  acres. 

The  concentration  of  BOD  in  separate  storm  sewers, 

BOD  ,  is  simply  the.  hourly  value  computed  by  STORM,  BOD, 

t 

(mg/1),  for  the  total  urban  runoff. 

Dry  Weather  Flow  and  Loadings 

,  Dry  weather  flow  and  BOD  loadings  are  assumed  known 
from  data  on  point  sources  in  the  area.   Thus,  Q.  represents 
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the  flow  (cfs)  into  receiving  waters  of  treated  wastewater, 
and  BODj  represents  the  BOD  concentration  at  68°F  (20°C) 
for  five  days,  mg/1.   The  amount  of  treatment  can  be 
varied  in  the  analysis. 

Combined  Flows  and  Loadings 

Dry  weather  flow  is  assumed  to  cause  only  a  negli- 
gible increase  in  flow  in  a  combined  sewer  during  a  storm 
event.   However,  two  factors  related  to  DWF  may  increase 
significantly  the  BOD  concentration  of  the  combined  sewer 
storm  water. 

(1)  The  BOD  strength  of  the  municipal  sewage  with 
which  it  mixes,  and 

(2)  the  BOD  exerted  by  sediment  accumulation  in 
each  section  of  the  sewer  under  DWF  conditions 
which  is  subject  to  the  "first  flush"  effect 
induced  by  the  initial  runoff. 

Data  collected  at  various  combined  sewer  overflow 
sampling  stations  in  Des  Moines,  Iowa,  support  the  first 
flush  theory  (Davis  and  Borchardt,  1974).   BOD  and  total 
suspended  solids  (TSS)  concentrations  decreased  with  time 
with  little  or  no  relation  to  the  flow  pattern.   Further- 
more, poll utographs  (BOD  vs.  time  and  TSS  vs.  time)  for 
these  stations  seem  to  indicate  that  the  flushing  occurs 


mostly  during  the  f i  rst  hour  of  runoff  generated  by  the 
storm  event.   A  model  based  on  shear  stress  considerations 
was  developed  to  predict  dry  weather  flow  deposition  for 
30  combined  sewer  collection  systems  in  Dorchester, 
Massachusetts  (Pisano,  1976).   The  overall  percentage  of 
DWF  suspended  sol  ids. loadings  depositing  daily  in  the  sys- 
tems was  estimated  at  10.3  percent- 

To  incorporate  the  first  flush  effect,  it  is  assumed 
that  the  hourly  in-sewer  sediment,  build-up  is  constant  over 
consecutive  dry  weather  hours  (DWH).   This  assumption  is 
reasonable  although  particle  size  and  specific  gravity, 
depth  of  flow  and  the  slope  of  the  conduit  are  important 
factors  affecting  deposition.   The  computation  of  combined 
sewer  flows,  as  with  separate  sewer  flows,  is  accomplished 
by  apportionment  of  the  total  hourly  flow  predicted  by 
the  urban  hydrologic  model  on  the  basis  of  relative  area 
served. 


Qc  =^  Qt 


(3.6) 


where   A  =  area  served  by  combined  sewers,  acres,  and 


Q  =  combined  sewer  overflow  rate,  cfs 
c 


The  BOD  load  contributed  from  accumulation  of 
sewer  solids  during  consecutive  DWH  is  lumped  into 
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the  first  hour  of  a  runoff  sequence,   and  is  estimated 
by 


FF  =  FFLBS  •  DWH 


(3.7) 


where 


FF  =  first  flush  BOD  load,  lbs. 


FFLBS  =  first  flush  factor,  Ibs/DWH,  and 


DWH  =  number  of  dry  weather  hours  preceding 
each  runoff  event. 


The  first  flush  factor  is  site-dependent  and  must 
be  determined  from   (1)  the  total  flow  generated  by  the 
combined  sewer  area,  including  the  DWF  contribution,  for 
all  periods  of  urban  runoff  during  the  year;  (2)  the  total 
number  of  DWH  during  the  year;  and  (3)  the  average  field 
measured  BOD  concentration  in  the  combined  sewer  overflows. 
For  Des  Moines,  Iowa,  it  was  estimated  from  digital  computer 
simulation  that  about  14  percent  of  the  average  BOD  load 
in  combined  sewer  overflows  may  be  attributed  to  first 
flush  effects.   Thus,  the  first  flush  factor,  in  pounds  of 
BOD  per  DWH,  may  be  obtained  from 


0.14  BOD 


FFLBS  = 


CM 


QCOM 


TDWH 


(3.8) 


where   BODp^  =  average  field  measured  BOD  concentration 


'CM 


in  combined  sewer  overflows,  mg/1. 
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QCOM  =  total  flow  generated  by  the  combined  sewer 
area,  including  DWF  contribution,  for  all 
periods  of  urban  runoff,  cf/yr. 


TDWH  =  total  number  of  DWH  in  the  year,  DWH/yr, 
and 


c,   =  factor  to  convert  mg/1  to  Ibs/cf 
=  0.000062435. 

The  total  flow  generated  by  the  combined  sewer  area 
is  obtained  from 

.   .  :r 
QCOM  =   Z    (Q^  +  DWFCMB).  (3.9) 

where   i  =  1 ,  2 ,  3,  . . . ,  r ,  and 


r  =  total  number  of  hourly  runoff  occurrences  during 
the  year. 


Finally,  the  mixed  BOD  concentration  in  the  combined 
sewer,  BOD  (mg/1),  is  computed  by  the  following  expression: 


BOD   .  ^°t:jc  "  '°°f  V^d  •  (V^t)^  f^^-  ^2 


Q,  ^   Qd  •  (A^/A^) 


(3.10) 


where   Cp  =  factor  to  convert  FF  from  Ib/hr  to  cfs  • 
■    mg/1  =  4.45. 
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The  first  flush  factor,   FFLBS,  is  computed  in 
the  next  section  from  available  data  sources.   Further 
research  is  necessary  to  develop  a  mathematical  relation- 
ship that  predicts  the  first  flush  BOD  load,  independent 
of  calibration  to  measured  data,  with  a  reasonable  degree 
of  accuracy. 

3.8  Data  Sources 

Numerous  agencies  provided  the  information  required 
to  meet  the  substantial  data  needs  for  mathematical  model- 
ing of  urban  runoff,  storage/treatment  systems,  and  receiv- 
ing water  effects.   All  land  use  classifications,  population 
density  figures,  area,  curb  lengths,  etc.,  were  obtained 
from  data  prepared  by  the  American  Public  Works  Associa- 
tion for  specific  use  in  STORM  simulations.   Precipitation 
records  for  1968,  collected  at  the  Des  Moines,  Iowa 
Municipal  Airport,  were  obtained  from  the  National  Weather 
Records  Center  at  Asheville,  North  Carolina.   Receiving 
water  upstream  flows,  temperatures,  BOD  and  DO  levels 
were  taken  from  the  Davis  and  Borchardt  report;  however, 
supportive  material  which  provided  considerable  insight 
on  the  characteristics  of  the  Des  Moines  River  was  obtained 
from  the  State  of  Iowa  Department  of  Environmental  Quality. 

The  total  hourly  urban  runoff  (Q^)  and  its  BOD 
concentration  (BOD.)  are  obtained  from  the  STORM  simulation 
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of  Des  Moines,  Iowa.   BOD^;,  Q^,  BOD^,  Q^,  BOD^,  and 

Q  are  computed  from  input  data  and  the  appropriate  mix- 
w 

ing  equations  presented  earlier.   BOD^  is  obtained  from 
the  mathematical  models  of  the  storage/treatment  process 
discussed  in  Chapter  V.  The  first  flush  factor  is  ob- 
tained from  calibration  of  the  digital  computer  model  to 
measured  data. 

(1)   From  rainfall  records. 


TDWH  =  6,993  DWH/yr 

(2)  From  equation  (2.10) ,  :.:/._/ ^^Z-   :;  -i:..  .^ :..-., 
QCOM  =  156,530,672  cf/yr  , 

(3)  From  Davis  and  Borchardt  { 1 974)  ,   - 
MD^^   =    72  mg/1 

(4)  From  digital  computer  simulation,  disre- 
garding in-sewer  solids  build-up,  the 
average  BOD  concentration  in  the  combined 
sewers  =  61 .92  mg/1 

From  the  observed  difference  between  items  (3) 
and  (4)  above, 

BOD  difference  =10.08  mg/1 

=  0.0006293  Ibs/cf 
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BOD  load 


(0.0006293  lbs/Gf).(156,530,672) 


98,511.76  Ibs/yr  (44,684  kg/yr) 
that  can  be  attributed  to  first 
flush  effects 


FFLBS 


98,511.76    Ibs/yr 
6,993    DWH/yr 


=    14.09    lbs    BOD/DWH 


=    6.39    kg    BOD/DWH 


This  factor  is  then  used  in  equation  (3.7)  to 
estimate  the  first  flush  BOD  load,  FF,  during  the  first 
hour  of  runoff  generated  by  each. storm  event.   First 
flush  effects  account  for  14  percent  of  the  average  BOD 
load  in  the  combined  sewer  overflows. 


CHAPTER  IV 
STOCHASTIC  CHARACTERIZATION  OF  URBAN  RUNOFF 

4.1   Hydrologic  Time  Series 

Before  the  unifying  concept  can  be  applied  to 
wet  weather  and  dry  weather  flows,  it  is  important  to 
examine  and  analyze  the  hydrologic  time  series.   A  time 
series  is  a  sequence  of  values,  arranged  in  order  of"  their 
occurrence,  which  may  be  characterized  by  statistical 
properties.   Data  representing  a  purely  random  physical 
phenomenon,  or  even  a  physical  process  that  contains  a 
random  element,  cannot  be  described  by  an  explicit  mathe- 
matical relationship.   A  hydro! ogic  time  series  may  be 
considered  to  be  the  sum  of  two  components:   a  random 
element  and  a  nonrandom  element.   For  a  hydrologic  time 
series,  the  occurrence  of  an  event  is  not  necessarily 
independent  of  all  previous  events  (Dawdy  and  Matalas, 
1964).   In  fact,  when  known  deterministic  components  are 
removed  (e.  g.  ,  annual  cycles)  the  dependence  between  hydro- 
logic  observations  decreases  with  an  increase  in  the  time 
base.   For  example,  the  dependence  between  monthly  obser- 
vations is  less  than  between  daily  observations  (Dawdy 
and  Matalas,  1964) . 

59 


60 


As  stated  previously,  rainfall  input  to  STORM  is 
prepared  as  a  sequence  of  consecutive  hourly  values  (in- 
cluding zeros  for  no  measurable  precipitation).   These 
inputs  are  used  by  STORM  to  generate  the  corresponding 
series  of  hourly  urban  runoff.   The  basic  approach  to  de- 
fine a  wet  weather  event  is  to  analyze  the  hydrologic  time 
series  and  establish  the  minimum  number  of  consecutive  dry 
weather  hours  (DWH)  that  separates  independent  storm  events. 
The  independence  of  these  storm  events  is  not  strictly 
climatological  and  is  discussed  later.   The  dry  weather 
hours  refer  to  periods  during  which  no  runoff  was  produced. 
Thus,  depression  storage  and  evaporation  rates ^must  be 
satisfied  before  any  runoff  is  generated  by  STORM. 

The  precipitation  time  series  for  Des  Moines,  Iowa 
for  the  year  of  record  1968  is  presented  in  Figure  4-1. 
The  abscissa  represents  the  10-month  period,  in  hours,  from 
March  1  to  December  30.   An  examination  of  the  rainfall 
record  provides  considerable  insight  as  to  storm  groupings, 
their  intensity  and  duration,  and  frequency  of  occurrence. 
The  broken  line  on  the  abscissa  indicates  dry  weather 
periods  at  least  nine  hours  in  length.   The  annual  pre- 
cipitation total  is  27.59  inches  (701  mm),  producing 
10.28  inches  (261  mm)  of  urban  runoff  as  predicted  by  the 
urban  hydrologic  model. 


4.2  Analysis  of  Time  Series 

The  analysis  of  the  hydrologic  time  series  is 
performed  at  two  levels:   (1)  an  analytic  approach. 
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autocorrelation,  and  (2)  an  approximate,  graphic  procedure. 
The  results  of  the  analytic  approach  are  subjected  to  a 
parametric  test  of  significance  for  serial  dependence. 

Two  approaches  may  be  used  to  compute  sample 
autocorrelation  coefficients:   a  circular  series  approach 
or  an  open  series  approach.   In  the  circular  series  approach 
the  end  of  a  sample  series  is  assumed  to  continue  with  its 
beginning  (Yevjevich,  1972).   For  both  continuous  and 
discrete  series  in  hydrology,  the  use  of  the  circular 
series  approach  is  not  recommended  because  the  first  and 
the  last  part  of  the  series  may  be  independent,  but  cor- 
responding adjacent  parts  may  be  hi  ghly  dependent.   The 
circular  series  approach  may  decrease  the  degree  of  de- 
pendence, thus  introducing  bias  into  the  computation 
(Yevjevich,  1972).   Thus,  it  is  preferable  to  estimate  the 
autocorrelation  coefficients  by  an  open  series  method. 


rj(k)  = 


n-k 


"-^  li=i  U  b=m  '] 


B -Mm"  lifi- MX '']T 


(4.1) 


where   r^(k)  =  sample  estimate  of  lag-k  autocorrelation 
coefficient  for  hydrologic  process  I, 


I 
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discrete  data  series  (observations)  of 
hydro! ogic  process  I,  for  i  =  1,  2,  ...,  n 


total  number  of  data  points  or  observations 
and      . 


=  number  of  hourly  lags 


In  the  open  series  approach  decribed  by  equation  (4.1),  the 
first  part. of  the  series  (n-k)  long  is  used  for  x.  and  the 
last  part  of  the  series  (n-k)  long  is  used  for  x.^.. 
Equation  (4. 1 )  provi des  an  unbiased  and  efficient  estimate 
of  the  1 ag-k  autocorrelation  coefficient,  which  is  also 
referred  to  as  the  sample  estimate  of  the  1 aq-k  seri  al 
correlation  coefficient  (Fiering  and  Jackson,  1971). 

A  plot  of  the  serial  correlation  coefficients, 
r(k),  against  the  number,  of  lags,  k,  is  called  a  corre- 
logram.   The  technique  of  autocorrelation  analysis  is  es- 
sentially a  study  of  the  behavior  of  the  correlogram  of  the 
process  under  investi  gation  (Quimpo ,  1968).   The  correlogram 
shape,  or  curve  joining  each  point  to  the  next,  is  hence- 
forth referred  to  as  the  autocorrelation  function.   An 
analysis  of  the  precipitation  time  series  of  Figure  4-1 
results  in  the  curve  shown  in  Figure  4-2.   At  a  lag  of 
of  zero  hours,  the  correlation  of  the  discrete  open  series 
is  unity  because  this  point  on  the  curve  represents  the 
linear  dependence  of  the  data  series  on  itself.   The 
number  of  observations  (including  zero  values)  totals 
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Figure  4-2.   Lag  k  Autocorrelation  Function  of  Des  Moines, 
Iowa,  Hourly  Rainfall,  1968. 
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Figure  4-2  continued 
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7,372  consecutive  values,  and  lags  up  to  720  hours  were 
investigated.   The  first  minimum  of  the  autocorrelation 
function  occurs  at  a  lag  of  10  hours,  and  the  value  of  the 
function  is  also  zero  at  this  point.   The  physical  in- 
terpretation is  that  periods  without  rainfall  for  at  least 
10  hours  separate  uncorrel ated,  and  therefore  independent, 
storm  events. 

Actually  any  point  of  the  autocorrelation  function 
which  lies  outside  of  the  95  percent  tolerance  limits 
indicated  in  Figure  4-2  suggests  a  significantly  non-zero 
correlation  between  storm  events  at  that  particular  time 
lag.   The  Des  Moines  rainfall  record  obviously  exhibits  : 
nonrandom  behavior  at  lags  of  377  hours  ('^  16  days)  and 
421  hours  ('x-  18  days)  in  particular.   The  tolerance 
limits  for  a  normal  random  time  series  of  N  values,  and 
an  open  series  approach  at  a  95  percent  probability  level, 
are  given  by  Yevjevich,  1972. 


TL  (95%) 


-1  ±  1  .645  V  N-k^l 
N-k 


(4.2) 


where   TL  (95%)  =  tol erance  .1 imi ts  at  a  95%  probability 
1 evel . 


As  the  number  of  lags,  k,  increases,  the  tolerance 
limits  diverge.  However,  the  divergence  is  not  noticeable 
for  large  N.   Values  of  the  autocorrelation  function  between 
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lags  of  100  to  300  hours  and  500  to  720  hours  fell  between 
the  95  percent  tolerance  limits  and  are  not  shown  in 
Figure  4-2. 

Similarly,  autocorrelation  analysis  was  performed 
on  the  sequence  of  hourly  runoff  values  generated  by  STORM 
from  the  rainfall  input.   The  lag-k  serial  correlation 
coefficients,  rQ(k),  are  plotted  against  the  number  of  lags 
in  Figure  4-3.   The  analytic  technique  establishes  that 
the  minimum  interevent  time  of  consecutive  DWH  that  separates 
independent  runoff  events  is  nine  hours.    Examination  of 
Figure  4-3  reveals  that  the  runoff  time  series  is  not 
purely  random  either.   Li  near  dependence  is  observed  at  :.  - 
time  lags  of  377  hours  ('v  16  days)  and  436  hours  ('v^  18 
days),  as  expected,  because  of  the  high  correlation  between 
rainfall  and  runoff  processes. 

Slight  differences  are  observed  between  the  correlo- 
grams  of  Figure  4-2  and  4-3.   These  are  due  to  the  fact  that 
depression  storage  and  evaporation  rates  must  be  satisfied 
before  runoff  is  generated  by  STORM.   Thus,  the  digital 
simulation  of  the  runoff  process  by  STORM  acts  as  a  filter 
and  has  a  slight  smoothing  effect. 

The  graphic  procedure  requires  the  number  of  dry 
weather  hours  immediately  preceding  each  hourly  runoff 
occurrence.   These  values  are  determined  directly  by  the 
chronological  record  provided  by  STORM  of  all  the  runoff 
events  it  generates  from  the  i nput  rai nfal 1 .   If  a  hydrologic 
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Figure   4-3.      Autocorrelation    Function   of  Hourly   Urban    Runoff 
for   Des   Moines,    Iowa,    1968. 
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Figure   4-3   continued 
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model  such  as'STdRM  i  s  not,  avai  1  aBl.e ,  a  close  approximation 
may  be  obtained  by,  assuming  that  the  same;  numbers  of  DWH 
precede  the  rainfall  and  runoff  events.   Thus,  the  infor- 
mation provided  by  the  precipitation  records  or  the  rain- 
fall time  series  (such  as  Figure  4-1)  is  sufficient.   A 
plot  of  the  number  of  wet  weather  events  obtained  by  vary- 
ing the  minimum  interevent  time  is  shown  in  Figure  4-4.  . 
It  is  evident  that  a  time  value  exists  after  which  an 
increase  in  the  minimum  intereivent  time  does  not  result  . 
in  a  correspondingly  significant  reduction  in  the  number 
of  storm  events.   By  the  graphic  procedure,  a  period  of 
eight  consecutive  DWH  i  s  obtained  as  the  minimum;-  i  nter- 
event  time.  This  result  may  be  checked  by  constructing 
a  divided  difference  table.   From  Table  4-1,  the  first 
minimum  second  difference  occurs  after  the  point  of  maxi- 
mum curvature,  at  about  8  hours.   These  results  indicate 
that  the  graphic  procedure  yields  an  approximate  solution. 
Whenever  possible,  the  analytic  approach  that  defines  a 
correlation  should  be  applied  to  investigate  the  sequen- 
tial properties  of. the  hydrologic  series. 

4.3  Definition  of  an  Event 


Based  on  the  above  analyses,  a  wet  weather  event 
and  its  duration  are  defined  in  the  mathematical  model 
as  follows: 
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Table    4-1 

Divided   Difference    Table    for   Determination 
of   Point   of  Maximum   Curvature 


Minimum  Interevent  No.  of  Divided  Differences* 

Time,  hours  Events  A  A^ 


2  82 

6.00 

3  76  3.50 

13.00 

4  63  3.83 

1.50 
6  60  0.25 

0.50 
8  59  0.04 

0.75 
12  56  0.06 

0.33 
15  55  0.03 


20  52 


0.60 


2 
♦Maximum   second   difference,    A    ,    defines   the   point   of 


maximum  curvature 
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1.  Any  runoff  occurrence  having  nine  or  more 
DWH  preceding  it  denotes  the  beginning  of 
the  event  (see  below). 

2.  The  event  continues  as  long  as  all  of  the 
subsequent  runoff  occurrences  have  a  DWH 
value  immediately  preceding  them  equal  to 
or  less  than  eight  hours. 

3.  The  event  runoff  duration  (in  hours)  is 
equal  to  the  sum  of  all  the  runoff  occur- 
rences in  (2) . 

4.  The  actual  duration  (in  hours)  must  be 
determined  by  examining  the  date  and  hour 
of  the  first  runoff  value  and  the  date  and 
hour  of  occurrence  of  the  last  runoff 
value  within  the  event. 

When  applied  to  the  hydrologic  time  series  of 

1968  for  Des  Moines,  Iowa,  the  above  definition  results 

in  a  total  of  58  storm  events  encompassing  65  days  and   ; 

ranging  from  an  event  runoff  duration  of  one  hour  to  an 

actual  event  duration  of  30.  hours.   The  longest  sequence 

of  consecutive  runoff  hours  is  a  wet  weather  event  24 

hours  in  length,  denoted  hereafter  as  Wet  Weather  Event 

No.  52.   The  precipitation  series  generating  this  event 

begins  at  a  time  of  6137  hours  on  the  abscissa  of  Figure 

4-1,  corresponding  to  a  late  December  storm.   The  response 

of  various  storage/treatment  systems  to  BOD  inputs  from 

this  runoff  event,  in  particular,  and  all  wet  weather 

events  in  general  is  presented  in  Chapter  V. 


4. -4  Statistical  Summaries  of  Deterministic  Mechanisms 

The  mathematical  relationships  that  explicitly 
describe  system  behavior,  or  response  to  time-varying 
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inputs,  produce  results  which  are  specific  to  those  in- 
puts;  It  is  therefore  helpful  to  complement  these  tech- 
niques with  statistical  measures  that  allow  the  generali- 
zation of  results  (Benjamin  and  Cornell ,  1970;  Freund, 
1971).   Thus,  a  better  perspective  may  be  obtained.   To 
illustrate,  it  is  of  interest  to  investigate  the  sta- 
tistics of  a  water  quality  parameter  entering  the  WWF 
storage/treatment  facility  shown  in  Figure  2-1 ,  and  also 
the  statistics  of  that  same  water  quality  parameter  after 
leaving  the  system.   Regardless  of  the  deterministic 
mechanisms  used  to  describe  the  input  to  and  output  from 
that  system,  the  relationships  develped- be  low. will  apply. 
In  such  a  manner,  the  effectiveness  of  various  system 
models  may  be  studied. 

For  a  sequence  of  observed  (or  computed)  values, 
the  sample  mean  is  simply  given  by 


»=i 


n 
1=1 


(4.3) 


where  x  =  unbiased  estimate  of  the  mean  value  of  the 
water  quality  parameter, 

X.  =  data  series  of  observed  magnitudes  of  the 
^   water  quality  parameter,  for  i=l,  2,  3,  ..., 
n,  and 

n  =  total  number  of  observations. 

The  sample  estimate  of  the  variance  is  computed 
conveniently  from 
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(4.4) 


where   s  =  unbiased  estimate  of  the  variance  of  ob- 
served magnitudes  of  the  water  quality 
parameter. 


The  standard  deviation  is  simply  the  square  root 
of  equation  (4.4).   The  coefficient  of  variation  is  ob- 
tained from  the  standard  deviation  and  the  mean 


c.  v, 


(4.5) 


where   c.v.  =  sample  coefficient  of  variation  of  the 
water  quality  parameter,  and 

s  =  standard  deviation  of  observations  of  . 
the  water  quality  parameter  from  its  com- 
puted mean. 


Statistical  measures  of  pollutant  concentrations, 
unlike  mass  rates,  require  flow-weighting  to  obtain  mean- 
ingful comparisons.   A  flow-weighted  average  is  obtained 
as  fol lows: 


Z  V  .  Q  . 
i=l  ^'      ' 
n 

1^  '^ 


(4.6) 


where   X  =  unbiased  estimate  of  the  flow-weighted 

^   mean  value  of  the  pollutant  concentration, 
and 

,Q.  =  flows  associated  with  observed  magnitudes 
^   of  the  pollutant  concentration,  L^/T. 
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Flow-weighted  variances  may  be  obtained  similarly;  how- 
ever, a  statistical  approach  based  on  recent  work  by 
DiToro  (1975)  is  presented  in  Chapter  VI  which  provides 
a  technique  for  calculation  of  the  coefficient  of  varia- 
tion of  concentration  for  the  well-mixed  constant  volume 
equalization  basin.   Thus,  variances  are  computed  in 
Chapter  V  only  for  flow  rates  and  mass  rates,  which  pro- 
vide sufficient  information  on  the  performance  of  the 
storage/treatment  devices. 


CHAPTER  V 

MATHEMATICAL  APPROACH  TO  STORAGE/TREATMENT  RESPONSE 
TO  VARIABLE  FORCING  FUNCTIONS 


5. 1   System  Models  and  the  Unifying  Concept 

The  applicability  of  the  one-dimensional  version 
of  the  conservation  of  mass  equation  to  various  components 
of  the  physical  system  has  been  discussed  in  Chapter  II- 
Three  mathematical  models  of  the  storage/treatment  unit 
are  investigated  throughout  this  chapter: 

-the  well-mixed  constant  volume  model, 
-the  well-mixed  variable  volume  model,  and 
-the  dispersive  variable  volume  model. 
These  mass  transport  models  are  derived  from  simplifica- 
tion of  the  same  differential  equation,  equation  (2.1). 
The  methodology  presented  for  completely  mixed  systems 
may  also  be  used  to  approximate,  for  example,  pollutant 
transport  through  lakes  and  bays. 

The  storage/treatment  unit  has  become  an  essen- 
tial  element  of  decision  perspectives  on  urban  stormwater 
control.   As  shown  in  Figure  3-4,  it  has  been  assumed  for 
modeling  purposes  that  the  urban  runoff  generated  by  pre- 
cipitation over  areas  served  by  separate  or  combined 

81 


82 


sewers  is  processed  through  a  single,  central  storage/ 
treatment  facil i ty  pri or  to  release  into  receiving  waters. 
As  a  subsystem,  it  is  acted  upon  by  external  forces  and 
as  a  result  produces  certain  effects.   Its  inputs  are 
referred  to  as  forcing  functions  and  output  as  responses. 
The  response  of  the  system  to  a  given  forcing  function 
depends  on  the  type  of  mathematical  model  chosen  to  de- 
scribe its  properties  and  behavior. 

Forcing  functions  may  be  continuous  functions  of 
time,  or  defined  only  at  discrete  intervals  of  time.   A 
discrete  function  may  be  treated  as  piecewise  continuous 
if  the  value  of  the  function  is  constant  for  each  time 
interval  where  it  is  defined.   Thus,  the  value  of  a 
discrete  forcing  function  may  change  abruptly  from  one 
time  interval  to  the  next.   In  that  sense,  a  discrete 
forcing  function  may  be  regarded  as  a  step  function  input. 
Regardless  of  the  form  of  the  input  function,  the  mathe- 
matical models  define  the  response  of  the  system  as  a 
continuous  function  of  time.   To  compare  effluent  quality 
and  quantity  among  the  three  mathematical  models  of  the 
storage/treatment  system,  a  common  parameter  was  selected: 
the  residence  time  of  inflows.   It  was  desired  to  select 
a  broad  range  of  residence  times  so  that  generalizations 
on  system  behavior  could  be  made.   In  sewage  treatment 
works  practice,  the  size  of  sedimentation  tanks  is  calcu- 
lated by  the  design  overflow  rate.   It  represents  the  unit 
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volume  of  flow  per  unit  of  time  divided  by  the  unit  of 
tank  surface  area  (American  Society  of  Civil  Engineers, 
1959).   The  selection  of  the  overflow  rate  is  a  matter 
of  engineering  judgment  and  experience,  but  usually  the 
recommendations  on  the  Ten-State  Standards  are  followed 
(Great  Lakes-Upper  Mississippi  River  Board  of  State  Sani- 
tary Engineers,  1971).   These  surface-loading  or  surface- 
settling  rates  range  from  600  gpd/ft   (211  1/day/m  )  to 
1000  gpd/ft^  (352  1/day/m^).   In  effect,  this  overflow 
rate  is  the  average  upflow  velocity  as  confirmed  by  ex- 
amining the  units  in  which  it  is  expressed.   Based  on 
average  design  flow  and  tank  depths,  overflow  rates  in 
current  use  result  in  nominal  detention  periods  of  2  to 
2.5  hours  (Metcalf  &  Eddy,  Inc.,  1972). 

In  industrial  water  pollution  control  practice, 
stabil ization  basins  are  the  most  common  method  of, organic 
waste  treatment  where  sufficient  land  area  is  available 
(Eckenfelder,  1966).   In  particular,  well-mixed  aerated 
lagoons  with  depths  of  6  to  15  feet  (1.8  to  4.6  m)  are 
designed  to  provide  2  to  10  days  of  detention  (Eckenfelder, 
1966).   However,  industrial  flows  in  an  urban  area  are 
much  smaller  than  stormwater  flows.   To  illustrate,  the 
combined  industrial  and  municipal  wastewater  flow  for 
Des  Moines,  Iowa,  is  approximately  55  cfs  ( 1 . 6  cu  m/sec). 
The  average  urban  stormwater  flow  (including  combined 
sewer  overflow)  is  approximately  1600  cfs  (45  cu  m/sec) 
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for  a  precipitation  year  close  to  the  national  average. 
This  represents  a  1:29  ratio.   Of  course,  this  ratio  will 
be  highly  variable  depending  on  population  densities, 
regional  climatology,  and  regional  industrialization. 
The  illustration  underscores  the  fact  that  the  surface 
area  required  to  provide  detention  times  on  the  order  of 
days  is  not  likely  to  be  available  in  Standard  Metropoli- 
tan Statistical  Areas  (SMSA's)  in  terms  of  both  space  and 
cost  of  providing  it.   Thus,  four  residence  times  were 
selected  for  investigation  for  each  system  model. 

-a  2-hour  detention  time, 
■    -a  6-hour  detention  time,         ::■::: 

-a  12-hour  detention  time,  and 

-a  24-hour  detention  time. 
In  addition,  a  48-hour  detention  time  was  studied  using 
the  variable  volume  model  to  represent  the  storage/treat- 
ment system,  and  results  are  presented  for  a  single  storm 
event  to  confirm  that  increased  flow  equalization  occurs. 

Although  applications  of  the  mathematical  models 
are  made  to  Des  Moines,  lov/a,  the  results  are  not  site 
limited.   As  mentioned  in  Chapter  III,  Des  Moines  and  the 
Des  Moines  River  represent  an  urban  area-receiving  water 
configuration  common  throughout  the  country.   The  data 
base  is  typical  of  that  available  to  engineers  involved 
in  studies  of  areawide  waste  management.   Generalizations 
are  discussed  in  the  summary  at  the  end  of  this  chapter. 
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5.2   Well-Mixed  Constant  Volume  Model 

Continuous  Forcing  Functions 

If  the  storage/treatment  facility  for  urban  run- 
off is  modeled  as  a  well-mixed  constant  volume  tank,  the 
system  may  be  represented  as  shown  in  Figure  5-1.   The 
volumetric  rates  of  flow  into  and  out  of  the  tank  are 
equal  and  constant.   However,  the  concentration  forcing 
function  and  response  are  continuous  functions  of  time. 
The  pollutant  mass  in  the  tank  is  assumed  to  undergo  a 
first-order  reaction.   From  a  mass  balance  of  the  dis- 
solved material  across  the  tank,  the  governing  equation 
is  given  by 


dc. 


Q  c^(t)  -  Q  C2(t)  -  K  C2(t)  V 


(5.1) 


where   V  =  volume  of  tank,  L  . 

Cj,(t)  =  concentration  of  material  in  the  tank  and 
outflow,  as  a  continuous  function  of  time, 
.M/L3, 

Q  =  volumetric  flow  into  and  out  of  tank, 
l3/T, 

c,(t)  =  concentration  forcing  function  of  material 
in  infl ow,  M/L^, 

K  =  first-order  decay  rate  constant,  1/T,  and 

t  =  time,  T. 


Equation  (5.1)  is  the  mathematical  expression  that 
describes  the  response  of  the  system,  C2(t),  to  a 
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Figure  5-1.   Well-Mixed  Constant  Volume  Model 


time -varying  forcing  function  of  concentration ,  c-j  (t) . 
Since  this  model  is  a  linear,  first-order  ordinary  differ- 
ential equation  with  constant  coefficients,  analytic  solu- 
tions are  not  difficult  to  obtain.   Rearranging  equation 
(5.1)  to  a  more  convenient  form. 


dc 
d 


V-'      [v  ^  kJc^U)  =J-c^(t) 


(5.2) 


and  letting 

equation  (5.2)  is  further  simplified  to 


dc 


TF  +  a  C2(t)  =^c^(t) 


(5.3) 


The  inverse  of  a  is  a  measure  of  the  time  required  by  the 
system  to  respond  to  a  change  in  the  input  concentration, 
c,(t).   The  general  solution  to  equation  (5.3)  is  ob- 
tained by  the  method  of  integrating  factors  (Ross,  1964) 
and  is 


Tadt  r^  /adt  „         "1 
C2(t)  =  e""^     Te^     ^c^(t)dt  +  bJ 


(5.4) 


where   3  =  constant  of  integration. 


By  examining  the  initial  condition  of  the  system,  it  is 


observed  that  at  t  =  0,  C2{t=0)  =  g,  and  equation  (5.4) 
is  reduced  to  the  form 


C2(t)  =  f e 


=  Ve"^*/  c^(t)e^*dt  +  C2(t=0)e-^*       (5.5) 


At  this  stage,  complete  solution  of  equation  (5.5)  is 
dependent  on  the  form  of  the  concentration  i'orcing  func- 
tion, c,(t).   This  continuous  function  of  time  may  be 
constant,  linear,  exponential,  or  of  other  types.   Solu- 
tions to  various  common  forcing  functions,  for  constant 
volume  and  constant  flow  rate,  are  presented  in  Table 
.5-1,. as  adapted  from  Rich  (1973).   It  is  important  to 
emphasize  that  these  solutions  pertain  to  continuous- 
flow  systems.   Thus,  the  forcing  function  and  the  re- 
sponse have  been  expressed  as  continuous  functions  of 
time:-  c,(t)  and  C2(t)i  respectively.   The  reference  to 
function  of  time,  t,  has  been  deleted  in  the  derivative 
terms  for  notational  simplicity. 

The  above  analysis  may  be  extended  to  a  storage 
facility  maintaining  constant  volume,  but  subjected  to 
time-varying  forcing  functions  of  both  concentration  and 
flow.   Referring  to  Figure  5-1,  the  flow  rate  Q  is  now 
a  continuous  function  of  time,  Q(t),  and  the  governing 
equation  is  given  by 


dCr 

dt" 


Q(t)  [ci(t)  -  C2(t)l  -  K  C2(t)V  (5.6) 
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Table  5-1 


Solutions  to  Equation  (5.5)  for  Various 
Forcing  Functions  of  Concentration* 
(Rich,  1973) 


Type  of  Input   Form  of  c-jCt)   Output  Concentration  C2(t) 


Constant 


il  _5L  (l-e"^^)+C9(t=0)e 
a   V  ^ 


•at 


Linear 


c,  +  bt 


^  _Q_  d-e"^*)  + 
a   V 


b Q_  (l-e"^*-at) 

.2   V 


+  C5  (t=0)e 


•at 


r  4.-   1         ^    ^  +  bt 

Exponential     c,  e  — 


lL__Q_  (ei^*  -  e-"M 
a+b  V 


+  C2  (t=0)e 


at 


*  b  =  arbitrary  constant 


a  =  Q  +  K,  where  K  =  first-order  decay  rate  constant 
V 
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Again,  by  rearranging  into  a  more  convenient  form. 


^.[aM.K]c3,t,=Mti 


V   ^i(t) 


(5.7) 


and  letting 


(t)  =  Sf^  +  K,  then 


dc. 


+  a(t)  C2(t)  =  — ^ 


Q(t)  c^(t) 


(5.8) 


The  general  solution  to  equation  (5.8) ,  a  1 inear,  first- 
order,  ordinary  differential  equation  with  variable  co- 
efficients  is  also  obtained  by  the  method  of  integrating 
factors  (Ross,  1964). 


C2(t)  =  e 


-^    a(t)dt  r^     ^  a(t)dt  Q(t)  c^(t)  dt 


-/*  a(t)dt 


+  C2(t=0)  e 


(5.9) 


where    f^   ^^^^^^ 
-  0 


=  integrating  factor, 


The  complete  integration  of  equation  (5.9)  is  dependent 
on  the  form  of  the  two  forcing  functions,  Q(t)  and  c.j(t) 
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Solutions  to  various  types  of  both  forcing  functions  are 
presented  in  Table  5-2.   Integration  has  been  accomplished 
as  far  as  analytically  possible.   An  alternative  method 
is  expansion  into  series  and  integration  of  a  number  of 
terms  until  the  approximation  is  considered  adequate. 
However,  it  is  simpler  to  discretize  the  input  to  obtain 
a  numerical  solution  on  the  digital  computer. 

Discrete  Inputs 


Observed  data  provide  information  about  physical 
phenomena  and  may  be  broadly  classified  as  being  either 
deterministic  or  nondetermi ni stic  (random).   In  practice, 
no  physical  data  are  truly  deterministic  nor  truly  ran- 
dom (Bendat  and  Piersol,  1971),  as  discussed  previously. 
The  forcing  functions  of  flow  rate  and  concentration 
presented  in  Tables  5-1  and  5-2  are  explicit  mathematical 
relationships.   However,  information  systems  about  hydro- 
logic  phenomena  and  water  quality  processes  are  usually 
designed  to  deliver  parameter  measurements  at  discrete 
intervals  for  a  finite  period  of  time.   Typical  configura- 
tions of  these  inputs  are  shown  in  Figure  5-2. 

The  response  of  the  storage/treatment  system  to 
discrete  inputs  may  be  derived  from  the  analysis  performed 
for  continuous  forcing  functions.   From  equation  (5.9), 
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Table    5-2 


Analytic   Solutions    to   Equation    (5.9) 
for   Various    Forcing    Functions 
of   Concentration   and    Flow    Rate* 


Form  of 

Form  of 

c^{t) 

Q^lt) 

Linear 

Exponential 

C]  +  bt 

Qe^^* 

where 

b,  y.  may  be 

either  positive 

or  negative 

Exponential        Exponential 


c-,e  +  bt 


Qe 


+Yt 


Output   Concentration 


C2(t) 


Q  .+Tt 


exp   [-^(e^"'-l)-Kt] 
[exp  (-:^[{c^  +bt)  exp:{:^e 

+  Kt)  -  c^  exp  (:^)  -      [v 

/t 
exp  (:^  e""^*  +  Kt)  dt 

-/^exp  (^e""^*  +  kt)  tdt 
0  ^ 

=0)] 


+  C2  (t 


exp[-^(e^^^-l)-Kt]- 
[ci  exp  (- ;^)  [exp  (;^e-^^^  + 


Kt  +  bt)  -  exp  (;^) 


(K 


+  b)/*exp  [;^e^^*] 


exp  [{K  +  b)  t]  dt 


;]+  C2   (t=0)j 
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Table  5-2  continued 


Form  of      Form  of 
c^{t)        Q^(t) 


Exponential   Linear 
c,  e 


I.  I  iicai  p  —t 

Q  =  at  exp     -  ^  t  +  27  t^  +  Kt  I • 

where  a  may  be  [  (i  t  +  ^  t^  +  Kt  +  bt) 

positive  or  negative  I 

-  c^  -  c^   {K  +  h)f    exp  (^  + 
2I  t^  +  Kt  +  bt)  dt**  + 
C2  (t=0)| 


♦Constant  Volume  Storage/Treatment  System. 

**Note  that  when  a  is  negative  the  integral  is  in  the  form 
of  an  error  function,  allowing  simplification. 
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Figure  5-2.   Discrete  Inputs  of  Flow  Rate  and  Concentration 
at  Equal  Time  Intervals. 
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CpCt.  +  At)  -  exp 


r   r-t.+At  /  0.    \  1 

-X:      (r^44 


+  c^Ct.) 


(5.10) 


where   Q.  =  value  of  input  flow  rate  from  t.  to 
^    (t.+At),  l3/T, 

c.  =  value  of  input  concentration  from  t. 
■•    to  (ti+At),  M/l3, 

At  =  length  of  time  interval ,  say ,  1  hour, 
and 

t.  =  the  beginning  of  the  time  interval  for 
^   which  the  system  response  is  being 
eval uated ,  T. 


For  computational  convenience,  the  forcing  function  may 
be  treated  as  a  step  function  input.   Thus,  initial  con- 
ditions are  assumed  to  be  defined  at  time  t.  by  C2(t.), 
the  solution  at  the  end  of  the  previous  time  step.   This 
allows  simplification  of  equation  (5.10)  to 


j(t.  +  At)  =  exp  T-y^  |y-+K)dtl 


j^''exp|-(?i.K)t]dt 


(continued) 
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c^Ct.)  exp   -/  (/+  k]  dtj         (5.11) 


Equation  (5. 11  )i  integrates  to 


C2{t.  +  At)  = 


)  =  exp  [-(v-+  KJAtj 


Q.c.   /     rQ.At     "]    X 

irriv  H  [-f-+KAtJ  -ij 


+  C2(t.) 


(5.12) 


The  average  concentration  for  each  time  step  may 
be  obtained  by  integrating  equation  (5.12)  over  the  length 
of  the  time  step,  as  follows: 


VS-  ^  ^)  =  A^ /''^2^^)  ^^ 


(5.13) 


where   t  =  dummy  variable  of  integration. 

Substituting  equation  (5.12)  for  C2(t),  and  noting  that 
Cj,(t.)isaninitialvalue, 

^2<*i  ^  "'  =  HvC"  ^TW   [  1  -  exp  (  -  5jl  -  K.)]  d. 


and  simplifying  further. 
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V 

V 


C  (t   +  — )  =  ^  ^    -   5-l_   .  -J-  /     pxD  f  -  -1 

^2^^i    2^   0^.  +  KV    Q^.  +  KV    At  V^   ^^P  \^   V 

(5.15) 

Final ly,  integrating  equation  (5. 1 5) ,  the  solution  is  ob- 
tained. 

e(t  .^)  =  -^iiL    r_Vi-.V(t)l. 

.     2^^i        2^       Q^.  +  KV        [  Q^.  +  KV      '^2^^i^J 


_1 

At 


1   -  expr-^yi+  K^Atl 


Q.  +  KV 


(5.16) 


where       ^^(t.) 


C2(t,.At, 


concentration  value  obtained  by  equa- 
tion (5.12)  for  the  previous  time  step, 
M/l3,  and 


time-averaged  concentration  for  each 
time  step,  M/L^. 


When  values  of  equation  (5.16)  are  plotted  with  respect 
to  time,  the  average  value  for  each  interval  of  the  step 
function  is  assumed  to  occur  at  the  mid-point  of  the  time 
step.   Thus,  the  notation  (t.  +^)  has  been  used. 

Since  the  input  and  output  mass  rates  are  also 
of  interest,  these  are  computed  from  the  known  concentra- 
tion and  flow  rates.   The  discrete  mass,  rate  input  is 
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given  by 


W.  =  Q. 
1   ^1 


(5.17) 


where   W.  =  value  of  input  mass  rate  from  t.  to 
^    (t.  +  At),  M/T. 

Because  inflow  equals  outflow  during  each  time  step,  even 
though  the  flow-through  rate  may  change  from  one  time  in- 
terval to  another,  the  response  of  the  storage/treatment 
system  to  the  mass  rate  input  is  given  by 


W2{t.  + 


At 


2)  =  Qi 


Co(t,  +  % 


(5.18) 


At. 


where   W2(t.  +  -y)  =  time-averaged  mass  rate  output 

for  each  time  step,  M/T. 


Again,  when  values  of  equation  (5.18)  are  plotted  with 
respect  to  time,  the  average  value  for  each  interval  of 
the  step  function  is  assumed  to  occur  at  the  mid-point 
of  the  time  step.   If  the  flow  rate  and  concentration  in 
equations  (5.17)  and  (5.18)  are  not  expressed  in  compatible 
units,  an. appropriate  conversion  factor  is  used. 


Organic  Decay  and  Sedimentation  Rates 

Colston  (1974)  investigated  the  characteristics 
of  urban  land  runoff  from  a  1 . 67-square-mi 1 e  (4.33  sq. 
kms)  watershed  in  Durham,  North  Carolina.   Thirty-six 
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separate  storms  were  sampled;  peak  discharges  varied  from 
2.25  to  1740  cfs,  and  the  number  of  dry  days  preceding 
each  storm  varied  from  0.5  to  34.   The  standard  BOD  test 
was  performed  at  several  dilutions  (from  0.5  percent  to 
5.0  percent  of  the  stormwater  concentration).   The  GOD 
test  was  performed  at  100  percent  strength  of  the  storm- 
water  concentration.   Both  the  BOD  and  COD  organic  decay  ■ 
rates,  Kj,  were  determined  by  the  method  of  moments.   This 
procedure  is  one  of  four  generally  accepted  methods  for 
computing  the  reaction  rate  of  the  BOD  curve  (Nemerow, 
1974).   There  are  numerous  problems  associated  with  plan- 
ning BOD  experiments  (Berthouex  and  Hunter ,  1 971  )  and  / 
obtaining  precise  parameter  estimates.   Colston  judged 
the  standard  BOD  test  as  inappropriate  for  determining 
the  ultimate  oxygen  demand  of  urban  runoff. 

However,  it  is  significant  that  the  COD  K^  rates 
varied  from  0.06  to  0.20  per  day,  while  those  determined 
by  the  conventional  BOD  test  variied  from  0.08  to  0;27  per 
day.   From  these  studies,  Colston  concluded  that  the  COD 
K^  exertion  rate  is  equivalent  to  that  obtained  by  a 
standard.  BOD  test.   The  average.  COD  decay  rate ,  at  1 00 
percent  sample  strength,  was  reported  to  be  0.13  per  day 
(base  e) . 

Based  on  results  reported  by  Colston,  the  organic 
decay  rate  for  BOD  prior  to  mixing  with  the  receiving 
waters  is  assumed  to  be 
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K^  =  0.005  hour 


-1 


(5.19) 


Thus,  while  stormwater  is  held  in  storage,  treat- 
ment will  occur  naturally  at  the  above  deoxygenation  rate. 
In  addition,  even  in  relatively  well -mixed  storage/treatment 
units  some  degree  of  sedimentation  will  occur.   Figure 
5-3  shows  the  BOD  removal  efficiency  of  primary  sedimen- 
tation tanks  treating  municipal  sewage  (modified  after 
Fair,  Geyer,  and  Okun,  1968).   The  exponential  relation- 
ship shown  in  Figure  5-3  is  accurate  for  BOD  influent 
concentrations  ranging. from  100  to  300  mg/1  ,  and  for 
relatively  long  rectangular  tanks  with  a  flow  regime 
approaching  a  pi ug  flow  pattern. 

In  the  calculation  of  the  BOD  washoff  rate  from 
an  urban  watershed,  STORM  includes  contributions  of 
10  percent  of  the  suspended  solids  washoff  rate  and  2 
percent  of  the  settleable  solids  washoff  rate  on  an  hourly 
basis  (Hydrologic  Engineering  Center,  1975).   The  washoff 
rate  of  suspended  solids  from  an  urban  watershed  may  be, 
depending  on  land  use  and  soil  cover,  quite  large.   Thus, 
it  is  possible  that  a  large  percentage  of  the  predicted 
BOD  load  in  urban  stormwater  runoff  is  derived  from  sedi- 
ment.  In  order  to  simulate  the  treatment  achieved  by 
settling  in  a  relatively  well-mixed  storage  facility, 
a  maximum  BOD  removal  of. 10  percent  is  imposed  on  the 
exponential  relationship  shown  in  Figure  5-3.   Then,  an 
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Figure  5-3.   BOD  Removal  from  Sewage  in  Primary  Sedimentation 
Tanks.   (Adapted  from  Fair,  Geyer,  and  Okun, 
1968.) 
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equivalent  sedimentation  rate  constant,  K  ,  is  derived  as 
f ol lows : 


1  -  exp  {-  K^  At)  =  0.10  [1  -  exp(-0.70  t^)] 


(5.20) 


where   At  =  length  of  time  step,  T,  and 
tp  =  detention  time,  T. 


Solving  for  K  ,  the  following  expression  is  obtained 


K  =  -  4-  £n[0.90  +  0.10  exp(-0.70  t^,)] 


At 


(5.21) 


where,  for  constant  volume 


'D 


Qi 


(5.22) 


The  organic  decay  rate  and  the  sedimentation  rate 
are  combined  into  an  overall  decay  rate  constant 


d    s 


(5.23) 


which  is  subsequently  used  in  equations  (5.12)  and  (5.16) 
to  obtain  the  output  concentration.   When  the  time  step 
length  and  the  detention  time  are  expressed  in  hours,  the 
units  of  K  are  hours"  . 

A  storage/treatment  unit  for  wet  weather  flows 
is  essentially  a  large  equalization  basin.   Its  function 
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is  to  dampen  extreme  fluctuations  of  concentration  and 
mass  flow  rate  by  providing  a  well-mixed  volume  with  a 
sufficiently  long  detention  time.   Equalization  is  pro- 
vided by  both  the  volumetric  mixing  of  the  influent  with 
the  wastewater  already  present  in  the  facility  and  the 
decay  action  in  the  system.   It  is  clear  from  equation 
(5.22)  that  the  nominal  detention  time  of  the  basin  will 
vary  according  to  the  influent  fluid  flow  rate,  Q..   An 
average  of  the  time  the  fluid  is  held  in  the  system  may 
be  obtained  by 


■  -     i  =  l 


(5.24) 


where   tr^  =  average  retention  time,  T, 

V  =  constant  volume  of  the  storage/treatment 
system,  L^, 

Q.  =  variable  influent  fluid  flow  rate, 
y      L3/T, 

1  =  1  ,  2,  3 ,  . . . ,  n ,  and 

n  =  total  number  of  time  intervals  of  equal 
duration. 


It  is  important  to  note  that  an  erroneous  average  deten- 
tion time  would  be  obtained  by  dividing  the  basin  volume 
by  the  computed  average  flow  rate. 
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Interevent  Decay  Process 

It  is  important  not  to  overlook  the  fact  that  the 
wastewater  in  the  storage/treatment  system,  between  storm 
events,  is  also  undergoing  a  first-order  decay  process. 
From  knowledge  of  the  number  of  DWH  between  storm  events, 
and  an  accountability  of  the  mixed  basin  concentration 
during  the  last  hour  of  runoff  of  each  storm  event,  the 
following  relationship  may  be  developed: 


c^d)  =  c^d-l)  exp  [-  K^  DWT(I)] 


(5.25) 


where  'C„(I) 

0 


c^d-i) 


concentration  of  wastewater  in  the 

storage/treatment  system  at  t  =  0, 
for  storm  event  I,  M/L^, 

concentration  of  wastewater  in  the 

system  during  the  last  hour  of  runoff 
of  the  previous  storm  event,  M/l3, 
and  . 


DWT{I)  =  number  of  DWH'  preceding  storm  event  I 


It  has  been  assumed  that  the  first-order  decay  process 
during  interevent  times  is. represented  adequately  by  the 
organic  decay  rate.   Thus,  any  sedimentation  occurring 
in  the  basin  during  dry  weather  periods  between  storms 
is  assumed  to  have  a  negligible  effect  on  the  BOD  concen- 
tration of  the  quiescent  wastewater. 
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Model  Appi ication 


A  description  of  the  study  area  has  been  pre- 
sented in  Chapter  III.   An  analysis  of  the  hydrologic 
time  series  and  a  statistical  definition  of  a  wet  weather 
event  are  discussed  in  Chapter  IV.   In  brief  review,  an 
urban  area  consisting  of  45,000  acres  (18,212  ha)  of 
separate  sewers  and  4000  acres  (1619  ha)  of  combined 
sewers  must  provide  storage/treatment  for  urban  runoff 
generated  by  a  total  annual  precipitation  of  27.59  inches 
(701  mm).   The  total  urban  runoff  predicted  by  the  hydro- 
logic  simulation  model  is  10.28  inches  (261  mm),  but 
hourly  runoff  occurrences  are  grouped  into  58  storm 
events  by  defining  a  minimum  interevent  time  of  9  dry 
weather  hours.   Input  and  output  hydrographs  and  polluto- 
graphs  are  presented  for  one  storm  event,  but  frequency 
analyses  include  all  runoff  occurrences  during  the 
year. 

Discrete  forcing  functions  of  flow,  BOD  concen- 
tration and  BOD  mass  rate  for  Wet  Weather  Event  No.  52 
are  shown,  respectively,  in  Figures  5-4,  5-5,  and  5-6. 
From  digital  computer  simulation  using  equation  (5.24) , 
the  volumes  (for  all  events)  required  to  maintain  average 
detention  times  of  2,  6,  12  and  24  hours  in  the  storage/ 
treatment  facility  are  presented  in  Table  5-3.   As  stated 
previously,  it  is  assumed  that  the  entire  urban  runoff 
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Figure  5-5   Input  BOD  Concentration  for  Wet  Weather  Event 
No.  52. 
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Figure    5-6.      Input   BOD   Mass    Rate    for  Wet   Weather   Event   No.    52, 
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Table  5-3 


Required  Basin  Volume*  for  Desired  Average 
Detention  Time** 


Average  Detention  Time  :^ 


Hours 


Basin  Volume  V, 
Cubic  Feet  (cu  m) 


2 

6 

12 

24 


4,300,000 
(121  ,776) 

12,900,000 
(366,328) 

25,800,000 
(730,556) 

51  ,600,000 
(1  ,461 ,312) 


*One  storage/treatment  facility  for  entire  urban  area. 
**As  defined  by  equation  (5.24)  for  all  wet  weather  events 


no 


generated  by  the  study  area  is  processed  at  a  single  cen- 
tral facility.   In  a  well-mixed  constant  volume  unit,  the 
time-varying  output  flow  rate  is  identical  to  the  time- 
varying  input  flow  rate,  as  shown  by  comparing  Figure 
5-4  to  Figure  5-7.   Figure  5-7  is  the  output  runoff 
hydrograph  regardless  of  the  basin  residence  time.   The 
response  of' the  well-mixed  constant  volume  model  to  BOD 
concentration  and  BOD  mass  rate  inputs  is  shown,  respec- 
tively, in  Figures  5-8  and  5-9  for  varied  average  deten- 
tion time. 

Figures  5-8  and  5-9  indicate  that  significant 
reductions  in  peak  BOD  concentrations  and  mass  rates  are 
obtained  by  increasing  the  residence  time  of  inflows  to 
the  storage/treatment  system.   The  residence  time,  of 
course,  is  increased  by  simply  enlarging  the  capacity  of 
the  holding  basin.   It  is  useful  to  examine  input  and  out- 
put statistics  for  the  storm  event,  presented  in  Table 
5-4.   Percentage  reductions  in  peak  BOD  concentrations 
and  mass  rates  and  mean  BOD  concentrations  and  mass  rates 
are  summarized  in  Table  5-5.   From  an  examination  of  the 
standard  deviations  of  the  BOD  mass  rates,  it  appears 
that  equalization  of  the  effluent  is  being  accomplished 
more  efficiently  as  the  average  detention  time  increases. 
An  alternative  approach  to  computation  of  concentration 
variances,  and  a  detailed  analysis  of  equalization  in  the 
well-mixed  constant  volume  model,  are  presented  in 
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Figure  5-7,   Output  Flow  Rate  for  Wet  Weather  Event  No.  52 
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Figure  5-8.   BOD  Concentration  Response  for  Wet  Weather 
Event  No.  52. 
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Figure  5-9.    BOD  Mass  Rate  Response  for  Wet  Weather 
Event  No.  52. 
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Table  5-5 

Pollutant  Removal  Effectiveness  of 
Well-Mixed  Constant  Volume  Model  for  Wet 
Weather  Event  No.  52 


%  Reduction  for 


Item 


tp=2  hrs   tp=6  hrs   ^^  =  12  hrs   tj^=24  hrs 


Peak  BOD 
Concentration 

28 

52 

68 

Mean  BOD 
Concentration** 

2 

14 

32 

Peak  BOD 
Mass  Rate 

-1* 

6 

27 

Mean  BOD 
Mass  Rate 

1 

14 

32 

79 
52 
51 
52 


♦Slight  increase  in  peak  BOD  mass  rate  is  observed  because 
the  time  sequence  of  inflows  and  concentrations  result  in 
pollutant  accumulation  and  the  detention  time  is  insuffi- 
cient for  decay. 


**From  flow-weighted  average 
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Chapter  VI.  From  Table  5-5,  it  appears  that  an  average 
detention  time  of  at  least  6  hours  would  be  required  to 
effect  adequate  pollutant  removal  for  Wet  Weather  Event 
No.  52. 

Displaying  pol  1  utographs  for  e^^ery   wet  weather 
event  during  the  year  would  provide  little  additional 
information  on  the  effectiveness  of  the  storage/treatment 
reactor.   Thus,  it  is  convenient  to  resort  to  statistical 
techniques  again  to  present  results  in  a  unified  format. 
Frequency  analyses  are  performed  on  BOD  concentrations 
and  mass  rates  for  every  hour  of  all  wet  weather  events. 
Figure  5-10  represents  the  frequency  distribution  of  all 
input  BOD  concentrations  to  the  WWF  storage/treatment 
facility.   These  concentrations  include  zero  values  for 
those  periods  of  dry  weather  within  the  storm  events 
(refer  to  Chapter  IV).   The  frequency  distribution  of 
all  BOD  concentrations  in  the  well-mixed  basin,  and  in 
the  discharge,  is  presented  in  Figure  5-11.   A  good  basis 
for  comparison  is  provided  when  cumulative  frequencies 
are  plotted  together,  as  in  Figure  5-12. 

The  effects  of  mixing  and  first-order  decay  in 
the  system  may  be  observed  by  comparing  the  cumulative 
frequency  curve  for  input  BOD  concentrations  to  the 
curves  representing  the  response  for  different  residence 
times.   As  the  detention  time  increases,  a  greater  number 
of  high  concentrations  mix  with  the  low  concentrations. 
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resulting  in  a  redistribution  of  magnitudes.   Similar 
analyses  are  performed  for  the  BOD  mass  rates.   The 
frequency  distribution  of  input  BOD  mass  rates  for  all 
wet  weather  events  is  shown  in  Figure  5-13.   Likewise, 
the  frequency  distribution  of  the  output  BOD  mass  rates 
is  presented  in  Figure  5-14  for  various  detention  times. 
Large  variations  in  maximum  BOD  mass  rates  and  the  higher 
BOD  mass  rates  precluded  the  use  of  a  common  scale  for 
all  graphs.   The  cumulative  frequencies  are  plotted  \ 
together  in  Figure  5-15.   A  relatively  uniform  reduction 
in  the  magnitude  of  all  BOD  mass  rates  in  the  output 
stream  is  observed  as  the  holding, time  increases .,  To  ;: 
supplement  these  figures,  input  and  output  statistics 
for  all  wet  weather  events  are  provided  in  Table  5-6. 
The  pollutant  removal  effectiveness  of  the  well-mixed 
constant  volume  mpdel  is  summarized  in  Table  5-7.   It 
seems,  again,  that  adequate  attenuation  of  extreme  mass 
flow  rates  requires  at  least  6  hours  of  detention  time. 
Reductions  in  mean  BOD  concentration  and  mean  mass  rate  ; 
are  also  quite  adequate  for  tp.  >_  6  hours. 

5.3  Well -Mixed  Variable  Volume  Model 


Continuous  Forcing  Functions 


The  storage/treatment  facility  for  wet  weather 
flows  may  be  modeled  as  a  well-mixed  variable  volume  basin 
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Figure  5-14.   Frequency  Distribution  for  Output  BOD  Mass  Rate 
for  All  Wet  Weather  Events. 
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Figure    5-14   continued 
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Table  5^7 


Pollutant  Removal  Effectiveness  of 
Well -Mixed  Constant  Volume  Model  for 
■All  Wet  Weather  Events 


Item 


%   Reduction  for  Given  Detention 

tp=2  hrs   tp=6  hrs   1^=12  hrs   tp=24  hrs 


Peak  BOD 
Concentration 

Mean  BOD 
Concentration  * 

Peak  BOD 
Mass  Rate 

Mean  BOD 
Mass  Rate 


31 

47 

23 

42 

_4** 

19 

12 

30 

59 
56 
30 

45 


72 
70 
49 
62 


*Flow-weighted  average. 

**Slight  increase  in  peak  BOD  mass  rate  is  due  to  pollutant 

accumulation  while  insufficient  detention  time  is  allowed 

for  decay,  because  of  the  time  sequence  in  flows  and  con- 
centrations. 
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in  which  case  the  system  may  be  represented  as  shown  in 
Figure  5-16.   The  volumetric  rates  of  flow  into  and  out 
of  the  system  are  unequal  functions  of  time.   The  pollu- 
tant mass  in  the  tank  is  assumed  to  undergo  first-order 
decay.   The  concentration  forcing  function,  c,(t),  and 
the  response,  c^Ct),  are  both  continuous  functions  of 
time.   From  a  flow  balance  across  the  basin,  the  time 
rate  of  change  in  volume  is  expressed  in  differential 
form  (Henderson,  1966). 


^     =     Q  (t) 
dt     ^r^^ 


Q2(t) 


(5.26) 


where   Q-i(t) 


Q2(t) 


fluid  flow  rate  into  tank,  as  a 
continuous  function  of  time,  L^/T, 
and 

fluid  flow  rate  out  of  the  tank,  as 
a  continuous  function  of  time,  L3/T. 


It  is  further  assumed  that  the  outflow  is  a  linear  func- 
tion of  the  volume  in  the  basin  to  allow  a  solution  to 
the  continuity  equation,  thus 


Q2(t)  =  aV(t) 


(5.27) 


where   a  =  time  constant,  1/T,  and 

V(t)  =  time-varying  volume  of  the  basin,  L" 


From  equations  (5.26)  and  (5.27),  an  expression  for  the 
outflow  as  a  function  of  the  inflow,  in  differential  form, 
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SYSTEM     SCHEMATIC 


Q,(t)C,(t) 


V(t),C2(t) 


Q|t)C^t) 


FLOW      BALANCE: 


MASS     BALANCE 


INPUT 


Q,(t)C,(t) 


Q^t) 


^^^    =  Q.(t)C.{t)-Q,(t)C^{t)-KC2(t)V(t) 
at 


V*-'    r    '         2""'   2' 


SYSTEM 


? 
OUTPUT 


Figure  5-16.   Well-Mixed  Variable  Volume  Model 


129 


may  be  obtained 


dQ, 


-^  +   aQ^Ct)  =  aQ^(t) 


(5.28) 


where  Q-iCt)  is  a  known  function  of  time.   Therefore,  equa- 
tion (5.28)  is  directly  integrable  by  applying  an  inte- 
.grating  factor  (Ross,  1954). 

QJt)=ae~''^r     Q,(t)e°'^  dt  +  Q„(t=0)e"''^        (5.29) 
0 

where   Q2(t=0)  =  initial  outflow  of  the  system,  L/T. 

From  a  mass  balance  of  the  dissolved  material 
across  the  holding  facility,  the  governing  equation  is 
found  to  be 


d(Vc2) 


^=   Q^(t)c^(t)  -  Q2(t)c2(t)  -  Kc2(t)V(t)       (5.30) 


From  equation  (5.27),  the  volume  in  equation  (5.30)  may 
be  expressed  as  a  function  of  the  outflow  rate.   Thus, 


lQ2(t)#/C2(t)[l  ^]  = 

QT(t)c^(t)  -  [^^]  Q2(t)c2(t 


)       (5.31 
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Rearranging  equation  (5. 31 )-  into  a  more  convenient  form: 
the  following  relationship  is  obtained: 


dc. 
It 


Q2( 


tj\dr)  ^^ 


+  K 


aQl(t)c,(t) 
^2(t)=     Q^U)  r      (5.32) 


Equation  (5.32)  is  a  linear,  first-order,  ordinary  differ- 
ential equation  with  variable  coefficients.   Letting 


P(t) 


Q,i 


tT  Vdt  / 


+  a  +  K 


(5.33) 


the  solution  to  equation  (5.32)  is 


/  P(t) 


dt 


C2(t)  =  e 


X 


t  /P^t^^^ai 


Ql(t)c^(t) 
Q2(t) 


dt 


+  C2{t=0) 


(5.34) 


Integration  of  equation  (5.34)  is  dependent  on  the  form 
of,  essentially,  three  forcing  functions:   c,(t),  Qi(t), 
and  Q2(t)  which  may  be  obtained  from  Q,(t).   For  non- 
trivial  forcing  functions,  integration  of  equation  (5.34) 
by  direct  analysis  is  not  possible.   Two  examples  are  pre- 
sented in  Table  5-8,   For  either  case,  a  numerical 
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Table  5-i 


Integration  of  Equation  (5.34)  for  Various 
Forcing  Functions  of  Concentration 
and  Flow  Rates* 


Form  of 
Forcing  Functions 


Integration  of  C2(t) 


Exponenti  al 

/.^  +bt 

c^(t)    =   c^e 


t        2.0        2^-at   ,   ^-at   ^   2. 
J       /  g    +ida-a    e         -dae         -2a    t 


0    ^       ~cit.«„-at      2.2,     2.  ,^    .  _ 

ae         +de         -a    t    +a    t  +  <:at-a-(i 


Triangul ar 
Q^(t)    =   4Qt(l-t) 


exp 

+a+    K)dt 

■   .                    t      2.„         2^-at    ,    „-at    ^    2. 
t  Pjt/g    +2a-a    e         -^ae         -^a    t 

Tq   exp    L    0^    „-at.,^-at      2.2.    2.  ,^    .        „ 
"  ae         +2e         -a    t    +a    t+2at-a-2 


] 


Combination 
Q2(t)  =  4Q(u|)(i) 

[at-l+e'^"^]  -  4Qt^ 


+a+  K)dt]. 

3,  /^   ^^^    ^   +  bt 
a  Ct  ( t  -  t  )  e 

[___! : 

(a+2)e"°'*-a^t^  +  a^t  +  2at-a-2 


•]  dt+C2(t  =  0) 


Exponential 
c^(t)  =  c^e-" 


t        2^-at,    |/„-at   ,    ^-at   ,|/„-a1 
r  r  ••    fO.   e        +aKe        -aae        -aKe 

XP     -l^n     I TT TT 


Exponential 


exp   -[^0    ( 
dt] 


„  at   ,„-at 
ae       -ae 


ft  rrt    .«^e-^^aKe-^^-aae-^^-aKe-"S 

[_  0  ae       -ae 


Exponential 

Q  (t)  =  5L_a  [e-3t. 
^w^  '   a-a  ■- 


dt] 
[ 


,    X    -{a+b)t 
a  (a-a)c,e  ^    ' 

1 ]  dt  + 

-at   -at     ■' 


ae   -ae 


C2  (t  =  ( 


♦Unequal  inflow  and  outflow  rates,  variable  volume  detention 
facility,  where 

Q2(t)  =  a  V(t). 
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approximation  of  the  original  differential  equation,  equa- 
tion (5.32),  or  a  numerical  integration  of  equation  (5.34) 
is  necessary. 

Discrete  Inputs 

The  response  of  the  well-mixed  variable  volume 
model  to  discrete  inputs  is  derived  from  the  basic  equa- 
tions developed  for  continuous  forcing  functions.   Forcing 
functions  are  assumed  to  be  step  function  inputs,  as  be- 
fore.  Thus,  from  equation  (5.29)  and  after  integration 
over  the  time  step  length.  At, 


Q2(t.  +  At)  =  Q.  [1  -  exp(-a  At)] 


+  Q2(t^.)exp(-a  At) 


(5.35) 


where   Q2(^i  ■*"  ^*^  "" 


"i 


Q2(t,) 


output  flow  rate  at  the  end  of 
the  time  step,  L3/T, 

value  of  input  flow  rate  from  t. 
to  (ti  +  At),  l3/T,  and 

outflow  value  at  the  end  of  the 
previous  time  step,  L3/t, 


The  average  outflow  rate  for  each  time  step  may  be  ob- 
tained by  integrating  equation  (5.35)  over  the  Tength  of 
the  time  step.   That  is. 


\^H'^^-   ae/"  ''2 <-)''- 


(6.36) 
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and  after  appropriate  substitution  and  subsequent  inte- 
gration, 

Q2(t.  +^)  =  Qi  +-^ -1_L  [expC-a  At)  -  1]       (5.37 


where   Q2(t.  +  ^) 


=  time-averaged  outflow  rate  for  each 
time  step,  L^/T. 


To  obtain  the  discrete  form  of  equation  (5.34) 
some  preliminary  analysis  is  required.   For  convenience, 
let  At  be  replaced  by  t,  Q.  by  a,  Q2(t^)  by  b,  and 
Q2(t^  +  At)  by  Q2(t)  in  equation  (5.35).   Thus,  rearranging 
terms , 


Q2(t)  =  a  +  b  exp(-aT) 


(5.38) 


Then, 


dQ.  ^ 

-j—^  =  -   b  exp(-aT) 
d  T 


(5.39) 


and  by  substituting  equations  (5.38)  and  (5.39)  into  (5.33), 


P(t)  =     -"  ^  ^'._  +  a  +  K 


a  +  b  e 


-0.x 


(5.40) 


It  is  now  possible  to  begin  the  derivation  of  equation 
(5.34)  by  solving  the  nonconstant  terms. 
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J      P(t)  dx  =  Jin 
0 


a  +  b  e 


-at 


Jin  (a  +  b)  +  at  +  Kt 


(5.41) 


Thus , 


exp 


-/     P(t)  dx]  =  ^  ^  ^.     exp  (-at  -  Kt) 

.0  J      a  +  b  e  "^ 


and 


[|' 


exp  \j       P(t)  dx 


a  +  b  e 


■at 


a  +  b 


therefore. 


dx  = 


/    exp[/     P(T)dx  I    Q-i^ 
0  Lo  -1  ^2^  -* 

1  r  p(a  +  K)t      ,1 

(a  +  b)(a  +  K)    L  ^  J 


(5.42) 


exp  (at  +  Kt)  (5.43) 


(5.44) 


Substituting  equations    (5.42)    through    (5.44)    into   equation 
(5.34),    replacing   a   and   b   by   their   original    variables,    and 
assuming   step   function   inputs,    the   discrete   form  of  equa- 
tion   (5.34)    is   given   by 


^2^ 


*i  "^  ^*^  ^  ^TT-K  [Q^.-[Q-Q2(V]exp(-a  At)  J 


[1  -  exp(-a  At-K  At)]     + 

Qo(tJc2(t.) 

7, TTri — n  il  v\ — I ITT  exp(-a  At-K  At) 

^i  "  l.(Qi-Q2^^i^-'®^P^~^  At)       "^^ 


(5.45) 
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The  average  concentration  for  each  time  step  is 
theoretically  obtained  by  integrating  equation  (5.45)  over 
the  length  of  the  ti me  step. 


- .,.  .At 
^2^^i    2'   At 


^-rJ,   /2(^^ 


dx 


(5.46) 


However,  an  exact  integration  of  equation  (5.45)  is  impos- 
sible.  Therefore,  the  average  mass  rate  output  is  ob- 
tained from 


t.  .  %  -  ^  /' 


•2-i  ■  2--  At  -    C2(T)Q2(T)dT 


(5.47) 


where  zAx)    is  given  by  equation  (5.45)  and  Q2(t)  is  given 

by  equation. (5. 35) .   This  expression  may  be  integrated 

because  some  complicated  terms  are  canceled  but.   The 
integration  of  equation  (5.47)  yields 


texp(-ctAt-KAt)-l]Q.c.  -  Q2(t.)c2(t^.)  • 

[exp(-aAt-KAt)  -  1] 

Wo{t.  +  ^)  =  time-averaged  output  mass  rate 
-         ^         for  each  time  step,  M/T. 


(5.48: 


where 


The  discrete  mass  rate  input  is  given  by  equation  (5.17), 
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as  for  the  constant  volume  storage/treatment  model ,   The 
average  concentration  for  each  time  step  is  then  obtained 
from  equation  (5.48)  as  follows® 


C2(t.  +4) 


Q2(ti  +  -i) 


(5.49) 


Unlike  the  constant  volume  model,  the  well-mixed 
variable  volume  facility  has  a  constant  detention  time. 


V(t.) 


'D  ~  Q^ITTy 


since  by  equation  (5.27) 


(5.50) 


V(^i)_   i 
'D   a  V(t.)    a 


(5.51) 


Thus,  the  volume  in  the  facility  varies  per  time  step, 
but  the  outflow  is  a  linear  function  of  the  basin  volume 
such  that  every  parcel  of  water, enteri ng  the  storage/ 
treatment  system  will  remain  in  the  system  for  a  time 
length  of  1/a. 

Interevent  Processes 


First^order  decay    between  storm  events  is 
described  by  equation  (5.25).   However,  when  the  storage/ 
treatment  system  is  represented  by  a  variabl  e  vol ume  model , 
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it  is  necessary  to  consider  the  outflow  rate  from  the 

system  several  hours  beyond  the  time  of  the  last  inflow. 

Thus,  initial  conditions  of  system  outflow  are  defined 
by 


Qgd)  =  Q2(I-l)exp[-a  DWT(I)] 


(5.52) 


where   Qq^^^ 


system  outflow  rate  at  t=0,  for  storm 
event  I ,  L3/T,  and 


Q  (I_l)  =  outflow  rate  from  the  storage/treatment 
system  during  the  last  hour  of  runoff 
of  the  previous  storm  event,  l3/T. 


Model  Appl ication 

As  stated  previously,  the  detention  time  of  inflows 
to  the  variable  volume  storage/treatment  system  is  con- 
stant. The  average  storage  volumes  used  by  the  holding 
facility,  for  fi xed.  detention  times,  were  obtained  from 
digital  computer  simulation  of  all  wet  weather  events 
and  are  presented  in  Table  5-9.   Comparing  Tables  5-3 
and  5-9,  it  appears  that  the  variable  volume  model  is 
considerably  more  efficient  in  storage  utilization  as  far 
as  providing  detention  time  to  inflows  is  concerned. 

Discrete  inputs  of  flow,  BOD  concentration  and 
BOD  mass  rate  for  Wet  Weather  Event  No.  52  have  been 
shown  in  Figures  5-4,  5-5  and  5-6.   The  response  of  the 
well-mixed  variable  volume  model  to  these  forcing  functions 
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Table  5-9 


Average  Storage  Volumes*  Used 
for  Given  Dentention  Times** 


Fixed  Detention  Time 
tp.  Hours 


Average  Basin  Volume  V, 
Cubic  Feet  (Cubic  Meters' 


2,287,368 
(64,778) 

2,746,323 
(77,776) 


12 


5,588,860 
(158,277) 


24 


11  ,154,820 
(315,905) 


*One  variable  volume  storage/treatment  facility  for  the 
entire  urban  area's  wet  weather  flows. 


**Defined  by  equation  (5.51) 
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is  shown,  respectively,  in  Figures  5-17,  5-18,  and  5-19 
for  varied  detention  time.   Figure  5-17  indicates,  as  ex- 
pected, that  significant  reductions  in  peak  flow  rates  out 
of  the  holding  facility  are  obtained  by  increasing  resi- 
dence times  of  inflows  to  the  system.   The  variable  volume 
basin,  however,  does  not  provide  significant  reduction  in 
BOD  concentrations,  as  noted  in  Figure  5-18.   On  the 
other  hand,  because  peak  flow  rates  are  reduced,  BOD  mass 
rate  peaks  are  correspondingly  decreased  as  observed  in 
Figure  5-19.   It  is  again  useful  to  examine  input  and  out- 
put statistics  computed  for  the  storm  event,  presented  in 
Table  5-10.   Pollutant  removal  efficiencies  for  peak,  BOD  . 
concentrations  and  mass  rates  and  mean  BOD  concentrations 
and  mass  rates  are  summarized  in  Table  5-11. 

Although  Figure  5-18  indicates  that  the  system 
response  is  practically  identical  for  varied  detention 
time,  a  flow-weighted  average  of  output  BOD  concentra- 
tions reflects  more  accurately  the  impact  to  the  receiv- 
ing body  of  water.   That  is,  the  flows  in  the  high  con- 
centration range  are  decreased  while  the  flows  in  the  low 
concentration  range  are  increased,  with  the  effect  becoming 
more  pronounced  as  the  detention  time  increases.  This 
■phenomenon  can  be  visualized  by  simply  comparing  the  hydro- 
graphs  in  Figure  5-17  to  the  pol 1 utographs  of  Figure 
5-18.   Thus,  it  is  not  surprising  that  when  the  flow- 
weighted  average  input  BOD  concentration  is  compared  to 
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VARIABLE    VOLUME 
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Figure  5-17.   Output  Flow  Rate  for  Wet  Weather  Event  No.  52 
for  Varied  Detention  Time. 
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Table  5-11 


Pollutant  Removal  Efficiency  of 

Well-Mixed  Variable  Volume  Model  for 

Wet  Weather  Event  No.  52 


%    Reduct 

ion  for 

Item 

tjj=  2-  hrs 

tp=6  hrs 

to=12  h 

rs 

tQ=24  hrs 

Peak  BOD 
Concentration 

3.3 

3.9 

4.4 

10.0 

Mean  BOD 
Concentration  * 

14 

38 

54 

65 

Peak  BOD 
Mass  Rate 

32 

,58 

73 

84 

Mean  BOD 
Mass  Rate 

48 

79 

84 

;  89   / 

*From  flow-weighted  average, 
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flow-weighted  averages  of  the  output  BOD  concentrations 
large  percentage  reductions  are  noted.   An  examination 
of  Table  5-11  suggests  that  excellent  reductions  of  peak 
and  mean  BOD  mass  rates  in  the  discharge  are  obtained  for 
tp.  >^  6  hours.   Percentage  reductions  in  peak  and  mean 
flow  rates  are  summarized  in  Table  5-12. 

Although  peak  and  mean  flow  rates  are  decreased 
substantially  for  detention  times  up  to  12  hours,  complete 
flow  equalization  is  not  achieved  in  the  strict  sense  that 
outflow  variability  is  not  reduced.   The  coefficient  of 
variation  of  the  outflow  rate  is  not  improved  over  the 
coefficient  of  variation  of  the  inflow  rate  except  for  the 
24-hour  detention  time.   Thus,  the  variable  volume  system 
response  for  a  detention  time  of  48  hours  was  also  inves- 
tigated.  The  results  are  presented  in  Table  5-13  for  Wet 
Weather  Event  No.  52  and,  as  expected,  both  flow  and  mass 
rate  equalization  are  achieved.   In  addition,  significant 
reductions  in  peak  and  mean  BOD  concentrations  are  ob- 
tained. 

Frequency  distributions  of  input  BOD  concentra- 
tions and  mass  rates  for  all  wet  weather  events  have  been 
presented,  respectively,  in  Figures  5-10  and  5-18  and 
apply  to  all  system  response  model s .   Frequency  analyses 
of  output  BOD  concentrations  and  mass  rates  are  presented 
in  Figures  5-20  to  5-23.   Figure  20  reflects  the  smoothing 
of  high-concentration  values  and  the  increase  in-  frequency 
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Table  5-12 


Reduction  in  Peak  and  Mean  Flow 
Rates  v/ith  Vai^iable  Volume  Model  for 
Wet  Weather  Event  No.  52 


Storage/Treatment 
System  Detention  Time 


%  Reduction  in 


Peak  Flow  Rate     Mean  Flow  Rate* 


2  hrs 

6  hrs 

12  hrs 

24  hrs 


5 
28 
50 
69 


38 
66 
66 
68 


*Averaged  over  period  of  discharge,  up  to  70  hours.  If  each 
were  averaged  over  the  same  overal 1-  time  period,  each  would 
have  the  same  mean  di  scharge. 
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Table  5-13 

Statistics  for  Variable  Volume  Model 
with  Detention  Time  of  48  Hours* 
for  Wet  Weather  Event  No.  52 


Item 


Peak  BOD 
Concentration 


Input  to  System 


Response   Reduction 


180  mg/1 


134 


26 


Mean  BOD 
Concentration** 

Peak  Flow 
Rate 

Mean  Flow 
Rate 

Standard  Deviation 
of  Flow  Rate 

Coeff ici  ent  of 
Variation  of  Flow 
Rate 

Peak  BOD  Mass 
Rate 

Mean  BOD 
Mass  Rate 

Standard  Deviation 
of  BOD  Mass  Rate 


63  mg/1 

18 

71 

2,267  cfs 
(64  cu  m/sec) 

424 
(12). 

81 

'  1 ,070  cfs 
(30  cu  m/sec) 

26  7 
(8) 

75 

771  cfs 
(22  cu  m/sec) 

,104   . 
(3) 

87 

0.72 

,,    0.39 

46 

38,499  Ibs/hr 
(17,462  kg/hr) 

3,356 
(1,522) 

91 

15,213  Ibs/hr 
(6,900  kg/hr) 

1,095 
(497) 

93 

14,262  Ibs/hr 
(6,469  kg/hr) 

926 
(420) 

94 

Coefficient  of 
Variation  of  BOD 
Mass  Rate 


0.94 


0.85 


10 


*From  simulation  of  all  wet  weather  events,  the  average 
storage  volume  used  equals  20,729,423  cubic  feet  (587,057 
cu  m)  . 

**Flow-weighted  average. 
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VARIABLE    VOLUME 
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Figure  5-21.  Frequency  Distribution  of  Output  BOD  Mass 
Rates  for  All  Weather  Events  for  Variable 
Volume  Storage/Treatment. 
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VARIABLE   VOLUME 
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Figure    5-21    continued 
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of  Tow-range  concentrations,  as  the  detention  time  increases, 
The  effect  is  si gni f icantly  more  pronounced  for  BOD  mass 
rates,  as  shown  i n  Fi gure  5-21  .   The  cumulative  frequency 
curves  for  output  BOD  concentration  in  Figure  5-22  also 
reflect  the  redistribution  when  compared  to  the  input  BOD 
concentration  cumulative  frequency  curve.   The  efficiency 
of  the  well-mixed  variable  volume  storage/treatment  system 
in  reducing  BOD  mass  rates  to  the  receiving  waterway  is 
reflected  in  Figure  5-23.   Input  and  output  statistics 
for  all  wet  weather  events  are  presented  in  Table  5-14, 
and  the  pollutant  removal  efficiencies  are  summarized  in 
Table  5-15  and  compared  to  reductions  in  peak  and  mean  flow 
rates.   Results  deviate  from  those  obtained  for  Viet  Weather 
Event  No.  52  in  that  reductions  in  peak  and  mean  flow 
rates  appear  too  high,  which  in  turn  explains  why  such 
high  reductions  in  BOD  mass  rates  are  obtained.   First  of 
all,  a  large  number  of  wet  weather  events  are  much  shorter 
in  duration  than  Event  No.  52,  but  flow  rates  into  the  stor- 
age/treatment system  may  be  much  higher.   The  methodology 
assumes  that  all  inflows  are  given  the  same  detention. 
Thus,  if  a  short  event  with  very  large  hourly  runoff 
values  is  processed  through  the  system,  it  will  occupy 
a  basin  volume  much  larger  than  the  averages  presented  in 
Table  5-9.   Table  5-16  lists  the  maximum  storage  used  by 
the  largest  hourly  inflows  for  each  detention  time  investi- 
gated.  Thus,  if  basins  are  sized  according  to  the  average 


154 


•r— 

1- 

■s- 

a> 

fO 

s: 

> 

4J 

« 

■o 

CU 

QJ  3 

X 

•r- 

-M 

^ 

<U 

1 

3 

,_ 

1 

cu 

r— 

3  < 

s_ 

S- 

o 

o 

**-  M- 

Ki-  ir>    00  o    o  vo 
r^  en    00  <id-    oo  co 


O  1— 

en  Lf) 

V£>  CM 

, — 

.—  00 

«^  in 

(£)<— 

o  r^ 

CM  <£) 

Lf>  «:i- 

CM 

r--  00 

CD  LO 

CO  IT) 

IT) — - 

1 —     . 

1 —       • 

Ln  UT) 

o  «d- 

«*  lO 

CM  ir> 

r^  o 

CO  iXJ 

O 

U3  cr> 

1—  CM 

CM  VC 

CM  «;i- 

00 

in  o 

«;*■>- 

1 — 

CM      . 

• 

o  in 

•s 

CO 

M3 

1 — 

n       r, 

^ 

.  ^-^ 

■    >^  , 

in  en 

VO  CM 

r^  r^  CM  in 
I—  o  <»  CO 
I—  in     CM  O 


o 


I—  CM      00^-- 

o  in    I—  cri 

CTl  CM     CO  -^ — 


CO  "^J-  00  CO  AO  CO 

r«.  CM  i~~»  r-~  CO  CO 

un  1 —  «*  in  o  r^ 

CO  in  CO  1 —  in  CM 

^  vo         —  — 


%."-- 

s-.^-^ 

s-^-^ 

^^ 

, ^ 

, . 

sz  s~ 

^  s- 

J=  s- 

o 

a 

o 

^V-C 

"«^^ 

c 

</)    Oi 

to  0) 

oi   OJ 

CO  -»-. 

V)  ~--. 

CO 

If-  (/) 

4-    (/) 

4-   CO 

X)    O) 

X»    CJl 

J3    C7) 

o  ■»», 

o  -~^ 

o  -^^ 

, — 

, — ii: 

1—  .ii! 

1—  j!!<: 

';l- 

E 

E 

E 

00 

00 

o 

<d- 

O 

• 

CO  en 

1^  r— 

O  00 

■    • 

«a-    3 

CO    3 

cn  3 

o 

VO  CM 

1 —  in 

in  r^ 

o 

cn  o 

in  O 

CM    O 

in  m 

r^  CO 

CM  CM 

M 

ft 

c 

A    •^ 

ff«    *« 

«    f% 

«=J-  ID 

r—  in 

1—  r^ 

in  en 

1—  en 

00  00 

CM  O 

<d- 

CO 

00  CM 

CM--' 

1 

O  •!->  O  ■!-> 

CO    C  CQ    C 

<U  <U 

JS^    O  CO 

<u  o  cu  o 

CL.  <_>  s  CJ 


LU 

^ii  cu 

fO  4-> 

Q-  Ci 


S-   3 


c:  cu      c:  u- 


cu   fO 
s:  o:; 


4J  *4- 


■  3 
o 

o 

•M    O 


•r-    O  O 

0  -r-  O 
•!--(->  CQ 

«+-    n3 

It-  -I-  CO    .i<;  CU 

01  S-  +->       (13  ,+J 

o  fo  rtj     cu  (O 

c_3  >•  Di    a.  cj: 


Q  -tJ 
O  res 

03  q: 


cu    CO 


4J  4J 
•r-   Qi 


-o  o 

C  CQ 
4->  4- 

oo  o 


o 

M-  CQ 

o 

4->    O 

cu  c  cu 

•r-     O   4J 

o  •>-  res 

•r-   4->  Qi 
M-    (O 

M-    •!-      lO 


cu 


155 


Table   5-15 


Pollutant  Removal  Efficiency  and  Outflow 
Rate  Attenuation  of  Well-Mixed  Variable 
Volume   Mode.l    for  All    Wet   Weather   Events 


%    Reduction  for  Given  Detention 

Item 

tp=2  hrs   tp=6  hrs   tp=12 

hrs   tQ=24  hrs 

Peak  BOD 

Concentration 

3          4          4 

4 

Mean  BOD 
Concentration* 

9         22         28 

36 

Peak  BOD 
Mass  Rate 

50         77         87 

■  93 

Mean  BOD 
Mass  Rate 

83         94       \'   95 

■   96 

Peak  Flow 
Rate 

64         82         90 

\-  94 

Mean  Flow 
Rate 

79         92         92 

92 

*From   flow-weighted   average 
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Table  5-16 


Maximum  Storage  Volume*  Required 
for  Desired  Detention  Time** 


Constant  Detention  Time 
tp),  Hours 


Maximum  Basin  Volume 
Used,  Vp^gx  cubic  feet 
(cu  m) 


2 

5 

12 

24 


64,091 ,360 
(1  ,815,057) 

95,283,376 
(2,698,425) 

108,058,528 
(3,060,218) 

118,376,288 
(3,352,416) 


*One  storage/treatment  variable  volume  facility  for  entire 
urban  area. 


**As  defined  by  equation  (5.51). 


157 


volumes,  it  can  be  expected  that  average  removals  would 
be  lower  for  the  entire  record  series.   In  essence,  it  is 
likely  that  these  short  events,  many  1  hour  in  duration 
only,  are  affecting  the,  statistical  averages  dispropor- 
tionately. 

5.4  Dispersive  Variable  Volume  Model 

The  previous  system  models  studied  represent  ideal 
mixing  conditions.   However,  nonideal  mixing  is  a  more 
likely  phenomenon  where  the  transport  of  a  pollutant  mass 
in  a  fluid  medium  is  described  by  advection  and  dispersion. 
Since  the  detention  time  in  a  constant  vol ume  system  varies 
according  to  the  influent  fluid  flow  rate,  the  flow-through 
velocity  is  al so  vari able .   Most  analytical  solutions  to 
the  convective  dispersion  equation  are  for  a  constant 
advective  vel  oci  ty  (Harl-eman  ,  1  971 ).   Thus,  the  model 
of  nonideal  mixing  is  developed  for  a  storage/treatment 
system  represented  by  a  variable  volume  facility  (i.e., 
constant  velocity).   The  dispersive  model  may  also  be 
applied  to  a  river  or  an  estuary. 

Continuously  Discharging  Plane  Source 

In  Figure  5-24,  a  plane  source  perpendicular  to 
the  flow  axis  continuously  injects  a  mass  of  pollutant 
into  a  variable  volume  storage/treatment  system.   A 
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uniform  longitudinal  velocity  advects  the  poll utant  mass 
in  the  direction  of  flow.   The  governing  equation  may  be 
derived  from  the  general  expression  for  the  one-dimen- 
sional convective  dispersion  equation  presented  in  Chapter 
II  as  equation  (2.1).   For  a  constant  longitudinal  disper- 
sion coefficient  and  first-order  decay  (sink  term),  equa- 
tion (2.1)  may  be  rewritten  as 


2 

oX 


(5.53) 


where   C  =  concentration  of  pollutant  in  the  fluid 
medium,  M/l3, 

-^  -    temporal  variation  in  the  pollutant  con- 
,  ^   centration,  M/L^T, 

E  =  longitudinal  dispersion  coefficient,  con- 
stant along  the  flow  axis,  L^/T, 

U  =  longitudinal  velocity  in  the  storage/ 
treatment  system,  L/T, 

f^  =  spatial  variation  in  the  pollutant  concen- 
^^        tration,  M/L4,  and 

K  =  first-order  decay  rate  of  pollutant  in  the 
fluid  medium,  1/T. 


Equation  (5.53)  is  a  linear  partial  differential  equa- 
tion, parabolic,  of  the  second  order,  and  is  formally 
equivalent  to  the  one-dimensional  heat  equation  (Berg 
and  McGregor,  1966).   The  fundamental  mathematical  tech- 
niques to  solve  the  Fickian  equations  of  heat  conduction 


160 


are  given  rigorous  treatment  by  Carslaw  and  Jaeger  ( 1 959) . 
The  most  powerful  technique  applies  the  Laplace  transform 
to  the  convective  diffusion  equation  to  remove  the  time 
variable,  leaving  an  ordinary  differential  equation  of 
the  second  order  the  solution  of  which  yields  the  trans- 
form of  the  concentration  as  a  function  of  the  space 
variable  x  (Crank,  1956).   Then,  by  inverse  transforma- 
tion the  solution  is  interpreted  to  give  an  expression 
for  the  concentration  in  terms  of  both  space  and  time 
variables,  x  and  t,  respectively. 

Such  a  solution  is  presented  by  Harl eman  (1 971 ) 
for  the  continuously  discharging  plane  source  by  a:  time^^  . 
integration  of  the  instantaneous  plane  source  solution 
for  a  constant  time  rate  of  injection.   Thus, 


c^{ 


X,t-T)  =J 


q"(T) 


0  .   p/  4  E-{t-T) 
K-(t-T)|  dx 


exp 


r  [x-U-(t-T 
["   4E-(t-T 


yy 


(5.54) 


where   q"(T)  =  variable  time  rate  of  mass  injection  per 
unit  area  at  the  plane  source,  M/l2t, 
and 

3 

p  =  density  of  the  receiving  fluid,  M/L  . 


If  the  time  rate  of  injection  per  unit  area  is  constant, 
then  the  solution  to  equation  (5.54)  for  any  time  t 


161 


greater  than  or  equal  to  the  duration  of  injection  t, 
is 


it"-'-*!'  =  2^'^'^PL^J 


I   V/l^Et/    L/4E.(t-t.  JJ   V^e; 


/  4E.(t-t^). 


-H7i)-"Liis]} -'(-§) 


./4E-(t-t^) 


(5.55) 


where   Cp(x,t-ti)  =  pollutant  concentration  at  distance 

X  along  the  flow  axis  and  time  t  >^  t-,  , 
mass  of  pollutant/mass  of  fluid, 
dimensionl ess , 


n  =  /u  +  4KE  ,  has  the  dimensions  of 
veloci  ty ,   L/T,  and 

t-j  =  duration  of  injection,  T. 

When  the  time  t  is  equal  to  the  duration  of  injection  t, 
equation  (5.55)  reduces  to 


{-(7i)-H(i) 


(5.56) 
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for  X  >  0. 

Equations  (5.r55)  and  (5.56)  are  the  response  of 
the  system  to  one  continuous  injection  for  a  duration 
equal  to  t.    It  is  of  interest  to  adapt  these  solutions 
to  the  case  where  a  succession  of  injections,  represent- 
ing different  mass  loadings  lagged  in  time,  occurs. 
First,  it  should  be  noted  that 


pQ^C^ 


(5.57) 


where   Q,  =  influent  fluid  flow  rate,  L  /T, 

C,  =  concentration  of  pollutant  in  the.  inflow, 
dimensionless,  and 

A  =  cross-sectional  area  of  the  plane  source, 
equal  to  the  width  of  the  storage/treat- 
ment system  times  the  depth  of  flow  in 
the  tank,  L^. 


Also,  since  the  system  has  a  variable  volume  and  a  con- 
stant velocity. 


U  = 


"ATtT 


constant 


(5.58) 


where   U  =  longitudinal  velocity,  L/T, 

Qi(t)  =  fluctuating  influent  fluid  flow  rate, 
'      l3/T, 

A(t)  =  fluctuating  cross-sectional  area  due  to  a 
variable  depth,  L^, 

L  =  length  of  the  storage/treatment  basin, 
units  of  L,  and 
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tp,  =  dispersive  variable  volume  system  deten- 


tion time,  T 


Therefore,  it  is  convenient  to  define  an  input  function 
to  the  system  as 


I(x)  =-^^ 


(5.59) 


where   I(t)  =  system  input  function,  M/L  ,  and 

Ci{t)  =  concentration  forcing  function,  dimen- 
sion! ess  . 


From  equation  (5.55)  a  system  response  function  may  also 
be  defined .  ;- 


f(t-T)     =     exp 


xU 
2E 


iHm- 


-i^MM^ 


{ 


V/4Et/ 


err 


x-n-(t-T) 
L    4E.(t-T) 


|exp(_-f) 


(5.60) 


where   f(t-T)  =  system  step  response,  dimensionless. 


Equation  (5.60)  may. be  described  as  a  weighting  function, 
which  is  defined  as  the  output  of  the  system  at  any  time 
t  to  a  unit  step  input  ending  a  time  t-x  before. 
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Discrete  Convolution 

The  essence  of  linearity  is  the  principle  of 
superposition  (e. g . ,  Dooge,  1973).   If  f(x)  represents 
the  output  to  an  input  x,  then  the  principle  is  satis- 
fied if  for  any  two  inputs  x, ,  x^  and  constantc, 

f{x^  +  x^)  =  f(x^)  +  f(x2)  (5.61a) 

f(cx)  =  cf(x)  (5.61b) 

Equation  (5.61)  represents,  respectively,  the  additive 
and  homogeneous  properties  of  a  linear  system  (Bendat  and 
Piersol,  1971).   The  principle  of  superposition  may  be 
applied  to  any  number  of  inputs.   Figure  5-25  illustrates 
a  summation  form  of  the  nature  of  discrete  convolution. 
The  dispersive  variable  volume  storage/treatment  system 
is  subjected  to  a  step  function  input  of  pollutant  concen- 
tration, I(t-).   Individual  step  responses  to  each  injec- 
tion are  characterized  by  f(t-T.).   It  should  be  noted 
that  the  effects  of  the  various  injections  are  appropri- 
ately lagged  in  time.   Thus,  the  general  form  of  the 
discrete  convolution  equation  is  given  by 


C(t)  =  H   I(T.)f(t  -  T.)  (5.62) 

1=1    ^       ^ 


165 


DISPERSIVE   VARIABLE   VOLUM] 


im 


cu) 


C(t)  =    21      I(T;)f(t-TJ 


Fiqure  5-25.   Discrete  Convolution, 
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where   C(t)  =  system  concentration  distribution  re- 
sponse to  step  function  inputs,  dimen- 
sion! ess  , 

i  =  1 ,  2,  3,  . . . ,  n, 

n  =  total  number  of  inputs,  and 

t.  =  iAt,  dimensions  of  time,  T. 

For  computational  convenience  the  step  function  input  is 
again  def'ined  in  terms  of  the  discrete  forcing  function 
and  constant  terms,  as  follows: 


I(t^.)  = 


U  C^(t.) 


(5.63) 


where   C-jCx.)  =  discrete  concentration  forcing  funC' 
^    tion,  dimension! ess . 


Since  the  system  response  is  of  interest  only  at  the 
storage/treatment  uni t's  outfal  1,  then  equation  (5.60) 
is  evaluated  for  x  =  L.   To  obtain  adequate  computational 
accuracy,  it  is  necessary  to  solve  equation  (5.60)  over 
a  wide  range  of  times  prior  to  its  invocation  in  the 
convolution  equation.   Any  given  storm  may  be  approxi- 
mated as  a  step  function  consisting  of  n  segments  of 
different  pollutant  injections  I(t^),  each  of  the  same 
duration  tn  (which  is  equal  to  the  time  step  length.  At). 
The  system  step  response,  equation  (5.60),  is  evaluated 
at  short  time  increments  Atp,  and  an  average  is  obtained 
over  the  time  step.  At. 
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f(t 


iJ 


'i' 


i  [ff  (*j-'i'* 


f(t.  +  i  At^  -  t^)   ...   +    ]-  f(t.  +  m  At2  -  t^)] 


(5.64). 


where        t.    =   any  time   t  ^  t,  , 

i    =    1 ,    2,    . . . ,    m-1  , 

m  =  ^rj-  =  total  number  of  increments  desired 
2   per  time  step, 

j  =  1 ,  2,  . . . ,  n,  , . . ,  N,  and 

N  =  any  desired  upper  bound. 

Thus,  a  large  array  of  values  of  the  function  f(t..  -  t, ) 
is  stored  in  the  digital  computer's  memory  for  subsequent 
recall  by  equation  (5.62).   Values  of  t-  up  to  100  hours 
are  evaluated,  using  a  computational  increment  Atp  =  0.2 
hours  in  equation  (5.64).   It  should  be  noted  that  equa- 
tion (5.62)  provides  the  system  response,  at  each  time 
step,  due  to  several  inputs  lagged  accordingly. 

Solutions  are  extended  well  beyond  the  time  of 
the  last  injection  of  urban  runoff  pollutants  for  two 
Important  reasons.   First ,  the  transient  effects  may 
still  be  substantial  several  hours  after  the  storm  event 
has  ended.   Second,  solutions  to  the  convective  disper- 
sion equation  are  limited  in  their  inability  to  handle 
initial  conditions  of  concentration  in  the  storage/ 
treatment  system  (prior  to  injection).   The  fluid  flow 
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rate  computations  presented  for  the  well-mixed  variable 
volume  model  are  applicable  to  the  dispersive  variable 
volume  model,  including  interevent  outflow.   The  disper- 
sive variable  volume  system  mass  rate  response  is  simply 
obtained,  for  each  time  step,  by 


W(t)  =  C(t)  •  Q2(t) 


(5.65) 


where  C(t)  is  defined  by  equatio-n  (5.62)  and  Q2(t)  by 
equation  (5.37) .  , 

Model  Application  . 

The  mathematical  models  representing  the  well- 
mixed  storage/treatment  systems  did  not  require  knowledge 
of  specific  physical  dimensions,  such  as  length,  width 
and  depth.   In  their  application  to  a  real  urban  site, 
the  models  were  able  to  predict  storage  volumes  from 
known  residence  times  and  inflows.   The  system  length 
is  required  to  apply  the  dispersive  variable  volume  model. 
From  an  establ ished  length  and  a  fixed  detention  time,  a 
constant  longitudinal  velocity  is  obtained. 

In  sewage  treatment  plant  design  practices,  the 
length  of  sedimentation  tanks  is  limited  to  a  maximum 
of  300  feet  (91.5  m)  for  rectangular  shapes  (ASCE,  1959). 
A  survey  of  42  plants  indicated  the  depth  of  tanks  ranged 
from  7.5  to  16.5  feet.  (2.3  to  5.0  m),  and  the 
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length-to-width  ratio  averaged  about  4.0/1  (ASCE,  1959). 
The  City  of  Milwaukee  evaluated  the  merits  of  detention 
tanks  as  a  practical  method  for  abatement  of  combined 
sewer  overflow  from  a  570-acre  (231-ha)  drainage  area 
(City  of  Milwaukee  et   aj[.  >  1975).   A  project  detention 
tank  was  constructed  which  measured  420  feet  (128  m) 
in  length,  75  feet  (23  m)  in  width,  and  had  a  depth  of 
16  feet  (4.9  m).   The  pollutant  removal  effectiveness  of 
this  facility  is  discussed  in  the  summary  at  the  end  of 
this  chapter. 

From  these  considerations,  the  length  of  the 
dispersive  variable  volume  system  was  fixed  at  300,  feet   ,: 
(91,5  m).   From  storage  volumes  reported  in  Tables  5-9 
and  5-16,  width  and  depth  requirements  may  be  determined 
according  to  the  dictates  of  good  design  practice,  engi- 
neering judgment,  and  experience.   It. is  recognized  that 
in  an  actual  implementation  several  detention  tanks  would 
be  required  for  practical  reasons.   However,  the  length 
assumed  in  this  analysis  is  realistic  and,  consequently, 
the  results  obtained  from  the  digital  computer  simulations 
should  prove  to  be  quite  helpful  in  selecting  adequate 
levels  of  detention  for  many  applications. 

The  selection  of  a  longitudinal  dispersion  co- 
efficient for  the  storage/treatment  system  proved  to  be 
more  difficult.   Measurements  of  this  important  parameter 
were  not  a  part  of  the  data  collection  effort  in  the 
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Milwaukee  study.   Murphy  and  Timpany  (1967)  reported  an 

2 
average  axial  dispersion  coefficient  of  5,730  ft  /hr 

(532  sq.  m/hr)  in  an  activated  sludge  aeration  tank 
66  feet  (20  m)  long,  30  feet  (9m)  wide,  and  15  feet 
(4.6  m)  deep.   The  tank  was  subjected  to  constant  air- 
flow from  sparger-type  air  diffusers  inducing  spiral 
mixing.   Thirumurthi  (1969)  suggested  that  a  dimension- 
less  dispersion  number  (t^E/L  )  of  0.25  was  reasonable 
for  a  waste  stabilization  pond  65  feet  (19.8  m)  long 

and  with  a  3-day  detention  time.   This  would  yield  an 

2 
approximate  dispersion  coefficient  of  15  ft  /hr  (1.4  sq. 

m/hr).   Recently,  Novotny  and  Stein  (1976)  reported  a 

2 
longitudinal  dispersion  coefficient  of  68  ft  /hr  (6.3 

sq.  m/hr)  in  a  relatively  well -functioning  equalization 

basin  at  a  batch  organic  chemical  wastewater  treatment 

plant.   The  basin  had  a  volume  of  133,648  cubic  feet 

(3,785  m"^)  and  an  average  inflow  of  T.22  cfs  (0.03 

3  3 

m  /sec  or  2,980  m  /day).   Thus,  the  average  detention 

time  provided  was  30  hours,  and  since  no  mixing  devices 
were  present,  equalization  by  dispersion  prevailed. 

Intuitively,  a  forced  aeration  tank  should  have 
a  much  higher  dispersion  coefficient  than  a  storage/ 
treatment  facility  with  no  mixing  devices.   In  fact, 
this  is  demonstrated  by  the  experimental  studies  of 
Murphy  and  Timpany  (1967),  where  the  di spersi on  coeffi- 
cient Was  found  to  increase  as  a  logarithmic  function  of 
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the  air  flow  rate.   Thus,  values  near  the  Novotny  and 
Stein  (1976)  measurement  should  be  expected  for  storage/ 
treatment  systems  characterized  by  nonideal  mixing.   A 
series  of  trials  was  conducted  by  solving  the  system 
response  function,  equation  (5.59),  for  various  disper- 
sion coefficients  and  the  length  selected  above.   Solu- 
tions were  obtained  with  a  dispersion  coefficient  of 

72  ft^/hr  (6.7  sq.  m/hr)  for  detention  times  of  6,  12, 

2 
and  24  hours.   However,  a  minimum  value  of  192  ft  /hr 

(17.8  sq.  m/hr)  was  required  for  solution  at  tp  =  2 
hours  because  of  numerical  problems  resulting  from  actual 
evaluation  of  equation  (5.60)  on  the  digital  computer. 
Discrete  inputs  of  flow,  BOD  concentration  and 
BOD  mass  rate  were  presented  earlier  in  Figures  5-4, 
5-5,  and  5-6  for  a  single  storm  event  having  a  24-hour 
runoff  duration  (Wet  Weather  Event  No.  52).   The  response 
of  the  dispersive  variable  volume  model  to  these  inputs 
is  shown  in  Figures  5-26  to  5-28.   The  output  hydro- 
graphs  presented  in  Figure  5-26  are  the  same  as  those 
predicted  by  the  wel 1 -mixed  variabl  e  volume  model.   As 
stated  previously,  the  same  flow  balance  equations  were 
applied  to  both  the  well-mixed  and  the  dispersive  cases. 
Figures  5-27  and  5-28  illustrate  the  transient  response 
to  BOD  inputs.   In  general,  substantial  reductions  in 
the  peak  BOD  concentrations  and  mass  rates  are  visible 
as  the  residence  time  of  inflows  is  increased.   An 
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exception  occurs  at  a  detention  time  of  6  hours.   The 
corresponding  outflow  hydrograph  segments  (rising  limb, 
peak,  and  receding  limb)  in  Figure  5-26  almost  coincide 
temporally  with  those  of  the  BOD  concentration  polluto- 
graph.  Figure  5-27.   Basin  dynamics  are  such  that  the 
peak  BOD  concentration  occurs  when  outflow  from  the  basin 
is  near  its  peak  discharge.   Thus,  the  peak  BOD  mass  rate 
for  tpj  =  6  hours  is  slightly  higher  than  that  at  tp.  = 
2  hour. 

Input  and  output  statistics  for  BOD  concentra- 
tions and  mass  rates  are  presented  in  Table  5-17  for 
the  single  storm  event.   Pollutant  removal  efficiencies 
for  peak  and  mean  BOD  concentrations  and  mass  rates  are 
summarized  in  Table  5-18.   From  these  tables  it  appears 
that  adequate  pollutant  removal  is  achieved  only  at 
detention  times  of  12  hours  or  higher. 

Frequency  distributions  of  BOD  inputs  for  all 
events  have  been  presented  in  Figures  5-10  and  5-13. 
Frequency  distributions  of  the  output  BOD  concentrations 
and  BOD  mass  rates  are  shown,  respectively,  in  Figures 
5-29  and  5-30.   It  is  clear  from  Figure  5-29  that  there 
is  a  shift  in  the  frequency  distribution  towards  the 
lower  range  BOD  concentrations  as  detention  time  increases 
The  same  effect  is  observed  in  Figure  5-30  for  output 
BOD  mass  rates,  although  apparently  more  pronounced. 
Cumulative  frequency  curves  for  BOD  concentrations  and 
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Table  5-18 

Pollutant  Removal  Efficiency  of 

Dispersive  Variable  Volume  Model 

for  Wet  Weather  Event  No.  52 


%    Reduction  for 

Item          t„ 

--1    hrs 

tp=6  hrs   tp=12  hrs 

tp=24  hrs 

Peak  BOD 

Concentration 

11 

21         41 

65 

Mean  Bod 

Concentration* 

0 

2         21 

57 

Peak  BOD 

Mass  Rate 

10 

-4**       34 

82 

Mean  BOD 

Mass  Rate 

10 

24         54 

85 

*From  flow-weighted  average. 

**Slight  increase  in  peak  mass  rate  because  the  peak 
BOD  concentration  occurs  when  outflow  from  the  sys- 
tem is  near  its  highest  discharge. 
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DISPERSIVE   VARIABLE   VOLUME 
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Figure    5-30 


20  30  40  50  60 

BOD  MASS  RATE,    1000    ibs/hr 

Frequency  Distribution  of  Output  BOD  Mass 
Rates  for  All  Wet  Weather  Events  for  Dis- 
persive Variable  Volume  Model. 
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Figure    5-30   continued 
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mass  rates  are  plotted,  respectively,  in  Figures  5-31 
and  5-32.   When  compared  to  the  input  BOD  concentration 
cumulative  frequency  curve,  it  is  clear  in  Figure  5-31 
that  substantial  incremental  reductions  in  output  BOD 
concentrations  are  achieved  for  each  increase  of  the 
system  detention  time.   These  reductions  are  even  more 
drastic  with  respect  to  BOD  mass  rates  as  shown  in  Figure 
5-32. 

Input  and  output  statistics  for  all  wet  weather 
events  are  presented  in  Table  5-19,  and  the  pollutant 
removal  efficiencies  for  peak  and  mean  BOD  concentrations 
and  mass  rates  are  summarized  i n  Tabl e  5-20,   Since  peak  . 
pollutant  loadings  are  important  factors  in  receiving 
water  quality,  the  performance  of  the  dispersive  variable 
volume  storage/treatment  system  appears  to  be  quite 
adequate  for  detention  times  greater  than  or  equal  to 
6  hours.   As  with  the  well-mixed  variable  volume  model, 
percentage  reductions  of  BOD  mass  rates  are  quite  high, 
due  to  outflow  rate  attenuation  (see  Tables  5-15  and 
5-16). 

It  is  important  to  examine  the  sensitivity  of 

the  dispersive  variable  volume  model  to  the  magnitude 

of  the  longitudinal  dispersion  coefficient,  E.   For  BOD 

inputs  from  a  single  storm  event  (see  Figure  5-5)  and 

constant  t  ,  K,  and  L,  the  BOD  concentration  response 

D 
for  different  dispersion  coefficients  is  illustrated  in 
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Table  5-20   , 

Pollutant  Removal  Efficiency  of 

Dispersive  Variable  Volume  Model  for 

All  Wet  Weather  Events 


%  RedL 

ction  for 

Item          tp=2  hrs 

tp.=  6  hrs 

*d"2 

hrs 

t„=24  hrs 

Peak  BOD 
Concentration     13 

34 

50 

70 

Mean  BOD 
Concentration*    37 

56 

72 

87 

Peak  BOD 

Mass:  Rate        12 

79 

91 

■97 

Mean  BOD 

Mass  Rate         72 

89 

95 

98  . 

*F1 ow-weighted  average 
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Figure  5-33.   The  results  underscore  the  need  for  choosing 
this  parameter  carefully  when  actual  measurements  are  not 
available.   It  is  obvious  that  both  peak  and  mean  BOD 
concentrations  from  the  storage/treatment  system  are 
reduced  sharply  by  arbitrarily  increasing  the  magnitude 
of  the  longitudinal  dispersion  coefficient.   For  example, 
when  E  is  increased  50-fold,  from  72  ft  /hr  to  3500 
ft^/hr  (6.7  to  334.8  m^/hr),  the  peak  BOD  concentration 
is  reduced  by  76  percent  and  the  mean  BOD  concentration 
by  57  percent. 

When  a  reactor  system  is  characterized  by  non- 
ideal  mixing  conditions,  such  as  prevail  in  the  disper- 
sive storage/treatment  case,  it  is  convenient  to  define 
a  dimensionl ess  parameter  (Wehner  and  Wilhelm,  1956). 


Pe  = 


^D^ 


(5.66) 


where   Pe  =  Peclet  dispersion  number. 

The  Peclet  number  is  the  inverse  of  the  dimensionless 
dispersion  number  adopted  by  Thirumurthi  (1969).   An 
ideal  reactor  with  fluid  flowing  through  a  long  rectangu 
lar  tank  (or  pipe,  freshwater  stream)  with  no  longitudi- 
nal dispersion  or  lateral  diffusion  is  characterized  by 
an  infinite  dispersion  number  (Pe  =  <»).   At  the  other 
extreme  is  the  limiting  case  of  a  fully  stirred  reactor 
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(activated  sludge  aeration  tank,  completely  mixed  deten- 
tion pond,  etc.)  with  infinite  dispersion  and  a  zero 
Peclet  number.   The  dispersive  storage/treatment  system 
is  analyzed  at  two  intermediate  stages,  Pe  =52  and 
Pe  =  104,  by  varying  the  three  parameters  that  define 
the  Peclet  number.   Figure  5-34  illustrates  the  system 
response  to  BOD  concentration  inputs  from  Wet  Weather 
Event  No.  52  for 

(1)  Pe  =  104  where  L  =1200  ft  (366  m), 
tp  =6  hours,  and  E  =  2304  ft^/hr 
(214  m^/hr); 

(2)  Pe  =104  where  L  =3600  ft  (1097  m), 
tp  =  24  hours,  and  E  =  5184  ft^/hr 
(481  m^/hr);  and 

(3)  Pe  =52,  tp  =  24  hours  where 

(i)    L  =  300  ft  (91.5  m),  E  =  72  ft^/hr 

■  2 

(60  m  /hr ) ;  and 

(11)    L  =  900  ft  (274  m),  E  =  648  ft^/hr 

(60  m^/hr);  and 
(ill)   L  =  3600  ft  (1097  m),  E  =10,368 
ft^/hr  (963  m^/hr). 
It  is  interesting  to  note  that  systems  with  identical 
Peclet  numbers  but  subjected  to  different  detention  times 
exhibit  widely  varying  responses.   Systems  with  different 
Peclet  numbers  but  subjected  to  identical  detention  times 
exhibit  responses  which  are  much  closer  together.   Finally, 
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systems  with  identical  Peclet  numbers  and  subjected  to 
identical  detention  times  exhibit  identical  responses. 

The  above  analysis  may  be  extended  by  examining 
input  and  output  statistics  for  BOD  concentration  for 
all  wet  weather  events,  presented  in  Table  5-21.   The 
same  conclusions  may  be  drawn,  so  it  appears  that  gen- 
eralizations are  possible  among  systems  of  different 
physical  dimensions  and  dispersive  characteristics. 
It  is  also  interesting  to  examine  the  system  response 
in  terms  of  another  dimensionless  parameter,  the  reaction 
rate  group  K  t^..   As  stated  previously,  biological  waste 
treatment  systems  are  designed  to  be  highly  reactive 
relative  to  the  detention  time,  and  thus  K  t^  >>  1. 
The  results  of  varying  both  parameters  that  define  the 
rate  group  are  presented  in  Table  5-22,  for  the  storage/ 
treatment  system,  for  all  wet  weather  events.   In  all 
cases  K  tj,  is  much  smaller  than  unity,  contrary  to 
secondary,  treatment  biological  waste  systems.   However, 
the  basin  mixing  effects  also  contribute  to  reductions 
in  peak  and  mean  BOD  concentrations.   Unless  very  high 
reaction  rate  coefficients,  K,  are  assumed  for  the 
pollutant  mass  it  appears  that  the  detention  time  pro- 
vided by  the  dispersive  storage/treatment  system  is  the 
dominant  factor  in  effluent  quality. 


191 


s- 

eu 

o 

, — 

«t- 

u 

OJ 

l/l 

Q. 

o 

•r- 

Vt 

+J 

Z3 

(/) 

O 

•r— 

•r— 

+J 

S- 

(O 

ra 

■M 

> 

on 

i- 

c 

o 

o 

<+- 

4->    (U 

c   > 
(U  LU 
O 
C    I. 

o  cu 

O  JC 
■!-> 
Q  (O 
O  (U 
CO  3 

4J 


E  to 


X) 
O) 


O  4J 


O  ■!-> 


192 


+J   o 


r—       .c  x: 

o  c>j 
o 


1 —      (/) 

s- 


CM 

1 — 

O 

II 

Q 

+J 

:^ 

CM 

• — 

4J 

O 

3 

11 

Q. 

Q 

C7 

+J 

•-• 

^ 

n 

c> 

u 

f- 

c: 

CM 

o 

r— 

•r— 

■ 

+-> 

o 

IJ 

II 

3 

a 

-o 

-t-> 

0) 

i^ 

q: 

&« 

■  A 

aj 

W 

c 

U3 

o 

O 

O- 

• 

irt 

o 

QJ 

II 

Ol 

o 

•4-» 

i>!i 

vo 

o 

o 

II 

o 

4-> 

i^ 

o 
o 


a 
o 


■M  E 


O 
O 


o 

+J 

£10 

c 

a) 

^ 

o 

(O 

c 

QJ 

o 

a. 

<_) 

s-    . 

o  vo 


S-  f*- 

O) 

<u 

+J 

Q.CM 

JZ 

trt  1 — 

o> 

•r- 

•1— 

•a  II 

(U 

5 

ta  uj 

1 

s 

S-  T3 

o 

O    C 

1 — 

193 


5.5  Summary 


A  rigorous  mathematical  approach  has  been  pre- 
sented to  study  the  response  of  three  storage/treatment 
systems  to  variable  forcing  functions.   All  of  the  models 
are  based  on  the  principle  of  conservation  of  mass.   Thus, 
each  system  model  is  represented  by  a  special  case  of  the 
continuity  equation,  and  its  solution  constitutes  the 
system  response.   Nevertheless,  intrinsic  differences 
exist  among  the  model  parameters.   For  example,  the  well- 
mixed  systems  are  characterized  by  a  volume  which  is 
constant  in  one  model  and  variable  in  the  other.   The 
variable  volume  models  are  characterized  by  completely 
mixed  conditions  in  one  and  nonideal  mixing  in  the  other. 
In  particular,  the  dispersive  variable  volume  model  re- 
quires a  length  dimension  and  a  longitudinal  dispersion 
coefficient.   Therefore,  the  response  of  the  various 
system  models  to  a  given  set  of  inputs  may  vary  con- 
siderably. 

To  illustrate  these  differences,  system  responses 
which  were  originally  presented  individually  are  now 
plotted  together  in  Figures  5-35  and  5-36  for  one  deten- 
tion time.   The  BOD  inputs  are  from  the  single  storm 
event.   For  the  well-mixed  cases,  an  immediate  response 
is  obtained  at  the  system  outfall  because  the  reactor 
contents  are  assumed  to  be  always  completely  mixed.  . 
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For  the  dispersive  case,  the  displacement  of  the  response 
curve  on  the  time  scale  is  a  function  of  the  system 
longitudinal  velocity  and  the  longitudinal  dispersion 
coefficient.   Thus,  a  response  at  the  system  outfall 
depends  on  the  length  of  the  reactor,  the  detention  time, 
and  the  magnitude  of  longitudinal  dispersion  in  the 
system.   For  example,  in  Figures  5-35  and  5-36  the  effects 
of  the. first  pollutant  injection  are  noticeable  at  the 
outfall  approximately  15  hours  later. 

The  BOD  concentration  curves  in  Figure  5-35  rep- 
resent actual  effluent  concentrations.   The  constant 
volume  storage/treatment  system  is  more  efficient  at 
reducing  peak  and  mean  BOD  concentrations  than  the  other 
systems.  In  any  well-mixed  constant  volume  basin  provid- 
ing sufficiently  long  detention  for  inflows,  concentra- 
tion equalization  is  accomplished  by  both  volumetric 
mixing  and  decay.   The  primary  function  of  variable 
volume  systems  is  to  provide  flow  equalization.   However, 
by  reducing  peak  and  mean  outflow  rates  the  pollutant 
peak  and  mean  mass  rates  from  these  systems  are  effec- 
tively decreased  as  shown  in  Figure  5-36.   The  same 
pattern  holds  for  the  other  residence  times  investigated, 
except  that  the  performance  of  all  three  system  models 
improves  steadily  as  the  detention  time  is  increased. 
The  high  peak  in  the  well-mixed  variable  volume  response 
curve  (Figure  5-35)  is  due  to  initial  conditions  in  the 
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basin.   These  conditions  cannot  be  accounted  for  by  the 
dispersive  variable  volume  model,  and  therefore  the  re- 
sponse for  a  high  longitudinal  dispersion  coefficient  (see 
Figure  5-33)  does  not  approximate  that  of  the  well-mixed 
variable  volume  system. 

Unquestionably,  whether  the  systems  are  well- 
mixed  or  dispersive,  the  dominant  factor  in  pollutant 
removal  efficiency  is  the  detention  time.   This  means 
that,  for  either  constant  volume  or  variable  volume 
models,  the  detention  facility  size  is  of  major  sig- 
nificance.  Based  on  modeling  studies  using  performance 
data  obtained  from  the  Humboldt  Avenue  combined  sewer 
overflow  detention  tank,  the  City  of  Milwaukee  et  al . 
(1975)  predicted  BOD  removal  efficiencies  of  about  35 
percent  for  a  tank  volume  of  1  million  gallons/sq.  mi 
of  drainage  area  (1461  cu  m/sq   km).   The  BOD  removal 
efficiency  increased  to  about  85  percent  for  a  tank 
volume  of  6  million  gallons/sq.  mi  (8766  cu  m/sq   km). 
These  results  were  projected  for  an  annual  precipitation 
total  of  28.9  inches  (734  mm).   Since  the  facility  simply 
overflowed  to  the  receiving  water  when  inflows  exceeded 
tank  capacity,  it  may  be  compared  to  the  well -mixed 
constant  volume  storage/treatment  model  applied  to  Des 
Moines,  Iowa,  for  which  an  annual  precipitation  total 
of  27.59  inches  (701  mm)  was  recorded. 
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For  an  average  detention. time  of  6  hours,  the  con- 
stant volume  model  achieved  a  mean  BOD  concentration  re- 
moval efficiency  of  42  percent  (see  Table  5-7)  for  a  tank 
volume  of  1.26  million  gallons  per  square  mile  of  drain- 
age area  (1842  cu  m/sq  km)  for  all  wet  weather  events 
during  the  year.   The  BOD  removal  efficiency  increased 
to  70  percent  for  a  tank  volume  of  5.1  mi  1 1  ion  gal  Ions 
per  square  mile  (7372  cu  m/sq  km)  which  maintained  an 
average  detention  time  of  24  hours.   The  close  agree- 
ment between  results  obtained  for  the  Milwaukee,  Wiscon- 
sin faci  1  ity  and  the  model  application  to  Des  Moines, 
Iowa  has  great  significance.   It  provides  a  certain 
degree  of  confidence  in  generalizing  the  findings  ob- 
tained throughout  this  chapter. 

Considering  the  importance  of  preventing  shock 
BOD  loadings  to  receiving  waters  as  well  as  providing 
adequate  reductions  in  mean  BOD  concentrations,  it 
appears  that  6  to  12  hours  of  detention  or  higher  should 
be  given  to  urban  runoff  flows  prior  to  discharge. 
The  requirements  will  vary  according  to  the  BOD  strength 
of  the  urban  stormwater  surface  runoff,  mixing  in  com- 
bined sewers,  and  receiving  water  conditions.   The  ef- 
fects on  receiving  water  minimum  dissolved  oxygen  levels 
produced  by  varying  the  detention  time  of  a  constant 
volume  storage/treatment  system  are  evaluated  for  the 
Des  Moines  River  in  Chapter  VII. 


CHAPTER  VI  ..    . 

STATISTICAL  APPROACH  TO  STORAGE/TREATMENT  RESPONSE 
FOR  WELL-MIXED  CONSTANT  VOLUME  MODEL 


6. 1  Introduction 

A  deterministic  approach  to  the  transient  response 
of  a  well-mixed  constant  volume  storage/treatment  system, 
based  on  the  unifying  principle  of  continuity,  has  been 
presented  in  Chapter  V.   The  statistical  approach  addressed 
herein  is  complementary  to  the  mathematical  solutions  to. 
equation  (2.8),  and  applies  only  to  a  system  governed  by 
equation  (2.8).   The  methodology . di scussed  in  this  chapter 
provides,  specifically,  a  computationally  convenient  measure 
of  equalization  of  time  variable  waste  loads.  When  the 
statistical  approach  is  not  used  in  a  complementary  role, 
some  serious  limitations  may  arise  depending  on  the  overall 
objectives. 

6.2  Equalization  of  Concentration  and  Mass  Flow  Rate 

Equalization  may  be  defined  as  a  means  of  dampening 
extreme  fluctuations  of  concentration  and  mass  flow  rate 
in  wastewater,  by. providing  a  storage/treatment  facility 
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with  a  sufficiently  long  detention  time  for  complete 
volumetric  mixing  and  adequate  decay.   Numerous  investi- 
gators have  addressed  the  need  to  provide  equalization 
facilities  for  wastewater  treatment  plants.   Effluents 
from  industrial  processes  often  contain  components  whose 
concentrations  are  highly  variable  in  time  (Wallace,  1968). 
Municipal  sewage  discharges  are  also  variable  phenomenons: 
extreme  low  flows  occur  between  2  and  6  A.M.,  and  the  maxi- 
mum flow  occurs  at  about  noon  (Metcalf  &  Eddy ,  Inc.  ,  1972) , 
Engineering  experience  has  shown  that  biological  treatment 
processes,  and  even  chemical  treatment  trains,  perform 
better  if  extreme  fluctuations  in  the  load  to  the  reactor 
system  can  be  avoided  or  damped  (Eckenfel der ,  1966).   Fail- 
ure to  account  for  the  unsteady  nature  of  influent  wfaste 
loads  can  cause  severe  deterioration  of  the  treatment  effi- 
ciency and  impairment  of  receiving  water  quality  (Novotny 
and  Stein,  1976) . 

An  examination  of  the  hydrologic  time  series  de- 
picted in  Figure  4-1  is  sufficient  to  confirm  that  extreme 
fluctuations  are  found  in  urban  stormwater  runoff.   The 
WWF  storage/treatment  system  (refer  to  Figure  3-4)  is  sub- 
jected to  variability  in  the  influent  mass  rate  and  also 
fluctuations  in  the  fluid  flow  rate.   Wallace  (1968)  and 
later  Novotny  and  Stein  (1976)  proposed  a  wide  variety  of 
analytical  procedures  to  characterize  output  variance. 
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These  techniques  revolve  around  the  definition  of  an 
equalization  function  which  is  the  ratio  of  the  variance 
of  the  output  concentration-time  curve  to  the  variance  of 
the  input  concentration-time  curve.   Novotny  and  Stein,  in 
particular,  use  power  spectral  analysis  to  describe  equali- 
zation of  concentration  for  completely  mixed  basins,  dis- 
persed flow  basins,  and  plug  flow  basins.   However,  both 
Wallace,  and  Novotny  and  Stein,  assume  a  constant  fluid  flow 
rate  through  the  reactor.  Thus,  such  techniques  are  un- 
suitable for  adaptation  to  urban  stormwater  runoff  storage/ 
treatment  systems  because  they  are  not  applicable  for  the 
case  of  varying  flow. 

In  a  recent  paper  by  DiToro  (1975),  the  author  pre- 
sents an  analytic  framework  that  addresses  the  design  of 
a  completely  mixed  constant  volume  equalization  basin  re- 
ceiving a  fluctuating  influent  mass  flow  rate  and  a  fluctu- 
ating influent  fluid  flow  rate.   Such  a  framework  is 
adopted  in  this  chapter  to  investigate  the  degree  of 
equalization  provided,  by  WWF  storage/treatment  systems 
for  the  well-mixed  constant  volume  case. 

6.3  Methodology 


It  is  quite  difficult  to  determine  the  actual 
probability  density  function  (pdf)  of  the  time-varying 
input  concentration,  C^(t),  even  for  simple  probabilistic 
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models.   The  analysis  is  complicated  even  further  if  con- 
sideration must  be  given  to  essential  hydrologic  variabili- 
ties.  These  variabilities  consist  of   (1)  heterogeneities 
of  catchment  properties,  (2)  heterogeneities  of  individual 
storm  properties,  and  (3)  random  variabilities  in  the 
climatic  time  series  of  storms  and  antecedent  moisture  con- 
ditions (Eagleson,  1972).   The  only  case  for  which  the 
pdf  is  known  is  for  constant  flow  rate  for  which  it  can  be 
shown  that  a  Gaussian  mass  rate  input  results  in  an  efflu- 
ent concentration  which  is  also  Gaussian  (DiToro,  1975). 
Thus  the  pdf  of  the  input  mass  flow  rate,  W^ (t) ,  is  ig- 
nored.  Instead,  the  inputs  are  described  by  basic  sta- 
tistical averages  and  exponentially  decaying  serial  corre- 
lation coefficients. 

The  relationships  that  follow  apply  to  a  well-mixed 
constant  volume  storage/treatment  system,  governed  by 
equation  (2.8)  or  (5.6),  expressed  as 


V  1^  =   W,(t)    -   Q(t)Co(t)    -   K  CpCt)    V 
dt  ' 


(6.1) 


where       W,(t) 


fluctuating   influent  mass   rate, 
M/T, 


Q(t)    =   fluctuating   influent   fluid   flow   rate, 
l3/T, 

C,(t)    =   time-varying  effluent   concentration, 
M/l3. 
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K  =  first-order  decay  coefficient,  1/T,  and 

V  =  time  invariant  volume  of  the  fluid  mass  in 
the  system,  L^. 


With  the  introduction  of  probabilistic  inputs  and  coeffi- 
cients the  problem  becomes  one  of  analyzing  random  differ- 
ential equations,  but  the  exact  solutions  for  such  problems 
are  seldom  available  or  quite  complex  so  that  an  approximate 
analysis  is  all  that  is  possible.   The  basic  statistical 
measures  have  already  been  presented  in  Chaper  IV.   How- 
ever, the  notation  was  of  a  general  nature,  and  it  is 
convenient  to  redefine  the  variables  again  as  follows: 


1=1 


(6.2) 


1=1 


(6.3) 


i=n 


(6.4) 


=  k1t  Z  [Q(t,)-Q] 

1  =  1  •- 


(6.5) 
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(6.6) 


(6.7) 


where   W,  =  mean  influent  mass  rate,  M/T, 

Q  =  mean  influent  fluid  flow  rate, 
l3/T, 

a  ^  =  variance  of  the  influent  mass  rate, 

^  m2/t2. 

Or}  =  variance  of  the  influent  fluid  flow  rate, 

^  L6/T2, 


a  =  standard  deviation  of  the  influent  mass 

^  rate,;  M/Ti   - 

o^  =  standard  deviation  of  the  influent  fluid 

^  flow  rate,  l3/T, 

V.  =  coefficient  of  variation  of  the  influent 


ma 


ss  rate,  dimensionless ,  and 


Vf,  =  coefficient  of  variation  of  the  influent 
^   fluid  flow  rate,  dimensionless. 


Although  the  lag  k  serial  correlation  for  the  influent  mass 
flow  rate,  ry(k),  may  be  defined  as  (DiToro,  1975) 


r,(k)= 


n-k 


n-k 

E 

i=l 


W^(t.)  -  W^ 


w^(V,)-Wi 


(6.8) 


actual  computation  herein  is  performed  by  using  equation 
(4.1)  as  recommended  by  Yevjevich  (1972)  and  Fiering  and 
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Jackson  (1971).  Assuming  exponentially  decaying  serial 
correlation  coefficients  (DiToro,  1975), 


Tw^M  =  ^^P 


kh 


(6.9) 


where   k  =  hourly  lag,  dimensionless ,. 
h  =  sampling  interval 

=  1. 1  -  t.  =  1  hour,  and 

T  =  correlation  time  constant,  for  mass 
^   rate  inputs,  hours. 

The  correlation  time  constant;  T^  is  essentially  the  time 
interval  that  separates  uncorrelated  random  influent  mass 
rates,  W^(t)  and  W^(t+T^).   The  value  of  T^  is  obtained  by 
fitting  equation  (6.9)  to  the  correlogram  obtained  from 
computation  of  equation  (6.8).   This  procedure  is  demon- 
strated in  a  later  application.   The  lag  k  serial  correla- 
tion coefficient  for  the  fluid  flovi  rate,  rg(k),  is  simi- 
larly defined  by  equations  (6.8)  and  (6.9)  using  the 
appropriate  substitutions.   Thus, 


TqO)  =  exp 


kh 


(6.10) 


where   T^  =  correlation  time  constant  for  the  influent 
^        fluid  flow  rates  ,  hours. 
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DiToro  (1975)  states  that  an  influent  fluid  flow 
rate  which  is  fluctuating  in  a  random  way  increases  both 
the  mean  effluent  concentration  and  the  effluent  concentra- 
tion variation  over  that  expected  with  no  flow  variation .  . 
It  is  reasonable  to  expect  an  increase  in  the  effluent 
concentration  variance.   The  explanation  suggested  for  the 
increase  in  the  mean  effluent  concentration  is  that  the 
average  value  of  the  ratio  of  two  random  variables, 
W,(t)/Q{t),  is  not  equal  to  the  ratio  of  the  averages. 
However,  DiToro  assumed  that  Q(t)  and  W^(t)  are  statis- 
tically independent.   The  author  of  this  work  demonstrates 
in  a  subsequent  application  that  these  variables  may  be 
significantly  correlated  when  dealing  with  urban  stormwater 
runoff.   Normalized  detention  times  are  defined  as 


n^tttA^ToT  '^•"' 


T,  =      ^D  -  (6.12 

^Q   Tq  (1  +  K  to) 


where   t  =  normalized  detention  time  for  the 

^       pollutant  mass  input,  dimensionless ,  and 

Tp)  =  normalized  detention  time  for  the  fluid 
"        flow,  also  dimensionless. 


207 


The  appearance  of  the  term  (1  +  K  tp)  in  the  denominator 
of  equations  (6.11)  and  (6.12)  is  from  a  well-known  result 
of  the  solution  to  the  conservation  of  mass  equation, 
shown  in  equation  (2.9). 

From  an  approximate  analysis  of  equation  (6.1) 
with  random  inputs  Wi(t)  and  Q(t),  and  assumption  of  the 
exponential  correlation  given  by  equations  (6.9)  and  (6.10), 
DiToro  defines  a  flow  variation  parameter 


(1  +  K  t  J  VTT" 


(6.13) 


where   m  =  equalized  variation  of  the  flow, 
dimensionless . 


Finally,  the  effluent  mean  concentration  and  the  effluent 
concentration  coefficient  of  variation  may  be  defined  in 
terms  of  the  flow  variation  parameter,  m: 


W 


h- 


Q  (1  +  K  tn)  (iV) 


,  m  <  1 


(6.14) 


and 


(1  +Tw) 


0.5 


,  m  <  1 


(6.15) 
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where   C.^  =  mean  effluent  concentration,  M/L  , 
^   and 

Vp  =  coefficient  of  variation  of  the  effluent 
concentration,  dimensionless . 


It  is  interesting  to  note  that  equation  (6.14)  is  identical 
to  the  steady-state  solution  to  equation  (2.8),  for  m  =  0. 
In  other  words,  m  is  merely  a  correction  factor  to  account 
for  the  decreased  equalization  in  the  basin  due  to  fluid 
flow  rate  fluctuations.   It  should  be  noted  that  equations 
(6.14)  and  (6.15)  are  approximate  in  that  they  apply  for 
m  <  1  with  the  approximation  deteriorating  as  m  ^  1 
(DiToro,  1975).   If  it  can  be  assumed  that  the  effluent^ 
concentration,  C2(t),  and  the  fluid  flow  rate,  Q(t),  are 
jointly  Gaussian,  the  coefficient  of  variation  of  the 
effluent  mass  flow  rate  may  be  computed  from 


1>  J-+JL+mV-2) 


2   2 


(6.16) 


where   v^q  =  coefficient  of  variation  of  the  effluent 
^   mass  flow  rate,  dimensionless. 


In  the  derivation  of  equations  (6.13)  through 
(6.15),  it  has  been  assumed  that  Q(t)  and  W-j(t)  are  sta- 
tistically independent.   The  cross-correlation  of  these 
two  variables,  at  lag  zero,  is  the  ordinary  (product- 
moment)  correlation  coefficient  and  may  be  estimated  by 
(Yevjevich-,  1972) 
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^  Z  W-(t.)Q{t.)  -4  Z!  W|(t.)  Z  Q(ti) 


WQ 


n 

0.5       .  . 

n 

r     n                      -| 

2 

i  =  l                               " 

1=1    .    • 

• 

(6.17) 


0.5 


n  r  " 

1  S  Q^(t,)--i  S  Q(t,) 
i=l        "  •-  i=l 


^2 


Equation  (6.17)  should  be  used  to  test  the  independence 
of  the  probabilistic  inputs. 


6.4  Application  to  Pes  Moines,  Iowa 

The  statistical  approach  discussed  in  the  previous 
section  is  applied  to  the  series  of  urban  runoff  inputs 
generated  by  precipitation  over  Des  Moines,  Iowa  (refer  to 
Figure  4-1 ).   The  response  of  the  we  11 -mixed  constant  volume 
storage/treatment  system  to  these  inputs  is  characterized 
primarily  by  the  mean  BOD  concentration  of  the  effluent  and 
its  coefficient  of  variation.   The  influent  BOD  mass  rate 
and  fl  uid  flow  rate  are  characterized:  by  temporal  averages 
presented  in  Table  6-1.   The;  correlation  time  constants, 
Ty  and  Tq,  are  obtained  by  fitting  exponential  curves  to 
the  correlograms  obtained  from  computation  of  serial  cor- 
relation coefficients. 
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Table  6-1 

Summary  Statistics  for  BOD  Mass  Rate  and 

Fluid  Flow  Rate  Inputs  of  the 

Urban  Runoff  Time  Series* 


Statistics 


Computed  Value 


Mean  Flow  Rate,  Q 

Standard  Deviation  of  the 
Flow  Rate,  Oq 

Coefficient  of  Variation 
of  the  Flow  Rate,  Vg 

Mean  BOD  Mass  Rate, 
Standard  Deviation  of  BOD 


Mass  Rate,  o 


W 


Coefficient  of  Variation  of 


BOO  Mass  Rate,  v 


U 


Correlation  Time  Constant  for 
Mass  Rate,  T|. 

Correlation  Time  Constant  for 
Fluid  Flow  Rate,  Tn 


1,584  cfs 
(45  cu  m/sec) 

1  ,290  cfs 
(37  cu  m/sec) 

0.81 


21  ,717    Ibs/hr 
(9,851    kg/hr) 


18,250    Ibs/hr) 
(8,278   kg/hr) 

0.84 
1.10   hours 
1 .40   hours 


♦Precipitation  year  of  record  for  Des  Moines,  Iowa,  is  196! 
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Figure  6-1  shows  the  correlogram  for  the  input  BOD 
mass  rate,  for  all  wet  weather  events.  Tolerance  limits 
at  the  95  percent  probability  level  (5  percent  probability 
of  committing  a  Type  I  error)  are  indicated.  As  stated  in 
Chapter  IV,  values  of  the  function  between  these  limits 
are  essentially  uncorrelated.  The  exponential  form  of  the 
serial  correlation  coefficient  function  is  approximately 


ry(k)  =  exp   [-  0.91  k] 


(6.18) 


where  T   is  found  to  be  1.10  hours.   Similarly,  from  the 
w 

correlogram  for  the  flow  rate,  shown  in  Figure  6-2,  its 
exponential  form  is  approximately  given  by 


rQ(k)  =  exp   [-  0.71  k] 


(6.19) 


where  Tq  is  found  to  be  1.40  hours. 

Equations  (6.13)  through  (6.16)  may  be  applied, 
from  the  information  given  in  Table  6-1 ,  to  the  storage/ 
treatment  system  for  the  residence  times  of  interest.  The 
first-order. decay  constant,  K,  is  determined  from  equations 
(5.19)  ,  (5.21)  ,  and  (5.23).   The  statistical  characteriza- 
tion of  the  storage/treatment  effluent,  and  other  pertinent 
parameters,  are  summarized  in  Table  6r2.   In  particular, 
the  response  of  an  equalization  basin  to  waste  inputs  is 
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Figure  6-1.   Lag  k  Serial  Correlation  Coefficients  for  Input 
BOD  Mass  Rates. 
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Figure, 6-2.   Lag  k  Serial  Correlation  Coefficients  for 
Input  Fluid  Flow  Rates. 
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Table  6-2 


Predicted  System  Response  to  Mass 
and  Fluid  Flow  Rate  Inputs 


Storage  Treat 

ment  Response 

Parameter 

^0=2  h 

rs   ^"^=6  hrs 

tp=12  hrs  tp=  24  hrs 

Pi  rst-Order  Decay 
Rate  Constant  *  K 

0.083 
hours" 

0.109 
^    hours"^ 

0.110  ,    0.110 
hours"     hours"' 

Dimensionl ess 
Fluid  Detention 
Time,  tq 

1.23 

2.59 

3.69       4.70 

Dimensionl ess  Mass 
Detention  Time, 

1.56 

3.30 

4.70       5.98 

Equal ized 

Variation  of  Flow, 
m 

0.47 

0.26 

0.16       0.09 

Mean  Effluent  BOD 
Concentration, 

C2 

67 
mg/1 

39 
mg/1 

27         17 
mg/1       mg/1 

Coefficient  of 
Variation  of 
Effluent  Concen- 
tration, Vp 


0.62 


0.4.6 


0.39 


0.33 


Coefficient  of 
Variation  of 
Eff 1 uent  Mass  - 

^cq' 


1.22 


1.00 


0.95 


0.91 


♦Increases . with  detention  time  due  to  increased  sedimentation 

of  suspended  solids. 
**Assumes  that  effluent  concentration  and  fluid  flow  rate  are 
jointly  Gaussian . 
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important  to  receiving  water  quality  in  two  aspects: 
reduction  of  both  mean  concentrations  and  concentration 
variabilities.   The  predicted  pollutant  removal  efficien- 
cies of  the  well-mixed  constant  volume  model  are  presented 
in  Table  6-3,  A  relative  degree  of  equal ization  may  be 
estimated  from  the  percentage  reduction  between  the  co- 
efficient of  variation  of  the  influent  BOD  mass  rate,  Vj^, 
and  the  coefficient  of  variation  of  the  effluent  BOD 
concentration,  Vq  (Novotny  and  Stein,  1976). 

From  Table  6-3,  it  appears  that  the  approximation 
provided  by  equation  (5.14)  deteriorates  for  a  value  of 
m  much  smaller  than  unity  when  the  average  detention  time 


"tpj  is  small.   By  inspection  of  equation  (6.13),  it  can  be 
seen  that  m  is  sensitive  to  the  magnitudes  of  the  average 
detention  time  and  the  normalized  detention  time  for  fluid 
flow.   The  statistical  model  predicts  an  increase  in  the 
mean  BOD  concentration  due  to  a  fluctuating  fluid  flow 
rate,  which  is  accounted  for  in  equation  (6.14)  by  the 
term  (1  -  m).  A  constant  flow  rate  would  result  in  an 
effluent  mean  BOD  concentration  of  52  mg/1.   However,  the 
statistical  model  also  predicts  a  coefficient  of  variation 
of  the  effluent  BOD  concentration  of  0.62  versus  0.84  for 
the  influent  BOD  mass  rate.   It  should  be  emphasized  that 
the  numerical  difference  between  Vp  and  v^^  is  not  an  ab- 
solute measure  of  the  degree  of  equalization.   Otherwise, 
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Table  6-3 


Pollutant  Removal  Efficiency  of  the  Well-Mixed 
Constant  Volume  Model 


Average  Detention  Time 


%  Reduction  in 

Mean  BOD 
Concentration* 


Equal izat ion 
[-V2C]100,% 


tpi  =  2  hours 


tj,  =  6  hours 


tp)  =  12  hours 


tp  =  24  hours 


•10 
36 
55 
72 


26 
45 
54 
61 


*The  statistical  average  of^  thie  influent  BOD  concentration 
is  computed  by  the  ratio  W-i/Q  and  is  equal  to  61  mg/1. 
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agreement  between  the  mean  effluent  BOD  concentrations  as 
predicted  by  the  deterministic  model  (Table  5-7)  and  the 
statistical  approach  (Table  6-2)  is  quite  good.   To  illus- 
trate, these  values  are  plotted,  over  the  same  range  of 
detention  times,  in  Figure  6-3. 

Close  examination  of  the  other  items  in  Table  6-2 
provides  more  insight  into  the  behavior  of  the  storage/ 
treatment  system.   If  the  detention  time  of  the  basin  (tp,) 
is  much  larger  than  the  correlation  time  of  the  fluid 
flow  rates  (T^),  then  the  dimensionless  detention  time  t^, 
is  large  and  the  flow  fluctuations  are  smoothed  by  the  basin 
mixing  effect.   For  example,  the  flow  variation  parameter,  : 
m,  decreases  as  Tq  increases  according  to  the  relationship 
in  equation  (6.13).   Thus,  for  an  average  detention  time 
of  24  hours  and  the  urban  runoff  time  history  for  Des 
Moines,  Iowa,  for  1968,  (1)  a  value  of  0.99  is  obtained 
for  the  term  (1  -  m  ),  so  that  (2)  the  flow  variation  has 
little  effect  on  the  computation  of  the  mean  effluent  con- 
centration in  equation  (6.14).   This  is  not  surprising  at 
all,  since  it  is  merely  a  mathematical  description  of  the 
fact  that  a  large  enough  storage/treatment  faci 1 ity  wil 1 
attenuate  the  extreme  fluctuations  in  the  influent  fluid 
flow  rates,  eliminating  interference  with  the  simultaneous 
process  of  mixing  and  first-order  decay  of  the  influent 
pollutant  mass  rates.   Of  course,  the  inverse  situation 
occurs  when  Tq  increases  relative  to  t"p. :   flow  fluctuations 
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WELL    MIXED    CONSTANT     VOLUME    MODEL 


DETERMINISTIC      MODEL 
(see    Table    5-r7  ) 

STATISTICAL    APPROACH 
(see    Table     6-3) 
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Figure  6-3.   Deterministic  Model  and  Statistical  Approach 
Predictions  for  All  Wet  Weather  Events. 
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become  a  significant  source  of  variation,  increasing  both 
the  mean  concentration  and  the  variation  of  the  effluent. 


6.5  Limitations  of  a  Strictly  Statistical  Approach 

First  of  al 1 ,  the  statistical  approach  does  not 
predict  the  peak  concentrations  which  may  be  expected  to 
enter  the  receiving  water,  nor  the  peak  mass  rates.  This 
limitation  may  not  be  critical  to  the  design  of  an  equaliza- 
tion basin,  but  is  very  important  to  study  the  effects  of 
shock  loads  on  dissolved  oxygen  levels  in  the  receiving 
body  of  water.   Of  course,  the  statistical  approach  does 
not  provide  a  transient  history  of  the  response  of  the 
storage/treatment  system  either.   Thus,  it  is  most  useful 
in  the  evaluation  of  equalization  effects  when  the  overall 
objectives  are  as  stated  in  Chapter  I. 

The  statistical  approach  is  further  limited  when 
applied  to  urban  stormwater  studies  in  another  respect: 
the  influent  mass  rates  and  the  fluid  flow  rates  may  be 
significantly  correlated.   Application  of  equation  (6.17) 
resulted  in  a  value  for  r^g  of  0.79;  consequently,  the 
assumption  of  independence  between  W^ (t)  and  Q{t)  is  not 
supported.   This  fact  does  not  appear  to  invalidate  the 
analysis  presented  thus  far  in  view  of  the  comparison  pre- 
sented in  Figure  6-3.   However,  it  seems  to  indicate  that 
probability  statements  concerning  the  effluent  concentration 
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would  be  risky.   Furthermore,  application  of  the  statis- 
tical approach  to  a  single  storm  event  proved  impossible. 
Figure  6-4  shows  a  plot  of  the  lag  k  serial  correlation 
coefficients  for  the  influent  fluid  flow  rates  of  Wet 
Weather  Event  No.  52.   Although  the  event  has  a  runoff 
duration  of  24  hours,  an  exponentially  decaying  curve  can- 
not be  fitted  to  the  correlogram.  The  widely  diverging 
tolerance  limits  indicate  the  region  of  essentially  uncor- 
related  coefficients,  not  significantly  different  from 
zero.  Thus,  a  much  longer  time  history  is  required  for 
analysis  than  provided  by  one  event.   Otherwise,  a  much 
smaller  time  interval  than  one  hour  is  required.  Observa- 
tions at  time  intervals  less  than  one  hour  are  difficult 
to  obtain  in  urban  hydrologic  studies  (McPherson,  1974). 
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Figure  6-4.   Lag  k  Serial  Correlation  Coefficients  for  Fluid 
Flow  Rate  for  Wet  Weather  Event  No.  52,  Hourly. 


CHAPTER  VII 
RECEIVING  WATER  RESPONSE  TO  WASTE  INPUT  COMBINATIONS 

7.1   Introduction 

The  choice  of  receiving  water  analyses  is  wide- 
ranging,  from  simply  examining  wastewater  concentrations 
at  the  inflow  points  to  the  use  of  detailed  hydrodynamic 
and  water  quality  simulation  models.   The  methodology 
presented    herein  is  consistent  with  the  research  objec- 
tives enumerated  in  Chapter  I  and  the  unified  approach 
described  in  Chapter  II.   The  receiving  body  of  water 
may  be  viewed  as  a  natural  storage/treatment  system 
and  may  be  represented  by  the  generalized  component 
depicted  in  Figure  2-1.   Simplified  mathematical  model- 
ing techniques  are  used  to  determine  critical  deficits 
and  resulting  minimum  dissolved  oxygen  (D.O.)  concentra- 
tions for  various  waste  input  combinations  and  stream 
conditions,  as  summarized  in  Table  7-1.   An,  abstraction 
on  the  physical  system  has  been  presented  in  Figure  3-4, 
which  is  useful  to  review  in  conjunction  with  Table  7-1. 

The  mathematical  modeling  techniques  adopted 
throughout  this  chapter  were  applied  by  the  author  to 
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Table  7-1 


Waste  Input  Combinations  and  Stream  Conditions 
(See  Figure  3-4) 


The  combined  effects  of  the  following  waste  sources  are 
modeled: 

1.  Upstream  river  flow 

2.  DWF  treatment  plant  effluent 

3.  WWF  storage/treatment  facility  effluent  from  separate 
sewer  and  combined  sewer  flow  Inputs. 

For  the  following  DWF  treatment  rates: 

1.  30%  BOD  removal,  primary  treatment 

2.  85%  BOD  removal ,  secondary  treatment 

For  a  well-mixed  constant  volume  WWF  storage/treatment 
facility  maintaining  the  following  average  detention  times 


1. 
2. 
3. 
4. 
5. 


D 


0,  no  storage/treatment 
2  hours 
6  hours 
12  hours 
24  hours 


And  the  following  stream  conditions,  measured  upstream  from 
the  point  of  urban  waste  Inflows  are  considered  for  each  time 
step: 

1 .  Ri  ver  flow 

2.  Initial  BOD  concentration 

3.  Initial  D.O.  concentration       , 

4.  River  temperature 
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Des  Moines,  Iowa  (Heaney  ejt  a1_.  ,  1 976)  to  simulate  the 
hypothetical  response  of  the  Des  Moines  River  to  the 
separate  and  combined  BOD  waste  inputs  from  (1)  upstream 
sources,  (2)  dry  weather  urban  sources,  and  (3)  wet  wea- 
ther urban  sources.   Emphasis  was  placed  on ' the  rel ative 
importance  of  these  waste  inputs  to  receiving  water 
minimum  D.O.  levels ,  as  an  example  of  a  possible  metho- 
dology for  areawide  wastewater  management  and  planning. 
The  presentation  of  results  was  of  key  importance,  and 
continuous  simulation  provided  a  means  by  which  cumula- 
tive D.O.  frequency  curves  (on  an  annual  basis)  could  be 
developed  to  demonstrate  the  deleterious  impact  of  wet 
weather  flows. 

the  emphasis  placed  herein  is,  as  summarized 
in  Table  7-1,  to  assess  the  adequacy  of  the  well-mixed 
constant  volume  storage/treatment  system  as  a  pollutant 
control  device  during  periods  of  urban  runoff.   The  WWF 
storage/treatment  system,  the  DWF  treatment  plant,  and 
upstream  river  flow  interact  to  produce  a  receiving 
water  response  below  the  urban  area,  as  a  function  of 
stream  conditions.   The  development  of  a  detailed  re- 
ceiving water  model  is  not  justified  for  the  problem 
context,  since  a  much  larger  data  base  would  have  been 
necessary  on  river  channel  geometry  and  quantity  and 
quality  measurements. 
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Some  key  assumptions,  most  of  them  typical  of  models 
for  interim  planning,  are  made  (Hydroscience ,  1971). 

1.  Temporal  steady-state  conditions  prevail, 
where  all  system  parameters  and  inputs  are 
constant  with  respect  to  time;  however,  a 
relatively  short  time  step  (1  hour)  is 
used  for  simulation. 

2.  Stream  system  parameters  (such  as  river 
flow,  velocity,  depth,  deoxygenation  rate, 
and  reaeration  rate)  are  spatially  constant 
along  the  flow  axis  throughout  each  time 
step.  '  .  , .: 

3.  All  waste  inflows  occur  at  one  point  on  the 
receiving  stream. 

4.  The  effects  of  various  natural  biological 
processes  (algal  photosynthesis  and  respira- 
tion, benthal  stabi 1 i zation)  are  incorporated 
by  D.O.  deficit  (if  none,  by  saturation) 
upstream  from  the  waste  inflow  point. 

5.  The  DWF  waste  treatment  plant  operates  at 
constant  efficiencies ,  independent  of 
hydraulic  and  organic  loadings,  for  the 
entire  period  of  simulation. 

7.2   Initial  Conditions 


InitiaT  condi tions  of  BOD  in  the  river  are  defined 
by  equation  (3.3).   Irt  subsequent  equations,  the  mixed 
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!0D  concentration  in  the  river  will  be  denoted  by  L 


Thus , 


0      m 


(7.1) 


The  assumption  that  all  waste  inflows  occur  at  one  point  is 
not  unreasonable  for  Des  Moines,  but  in  some  locations  the 
distribution  of  inflows  along  the  river  may  need  to  be 
considered-   It  is  important  to  emphasize  that  all  of  the  BOD 
contributors  in  equation  (2.4)  represent  BOD5  values.   Thus, 
the  mixed  BOD  concentration  in  the  river,  BOD^,  is  also  in 
terms  of  the  standard  BOD  test.   The  ultimate  first-stage 
(carbonaceous)  demand  is  related  to  the  BOD5  value  by 


(Lo)c 


BOD, 


1   -5K, 
1-e   I 


(7.2) 


where 


(L  )^  =  ultimate  first-stage  BOD  demand,  mg/1 , 


and 


K^        =      first-order   BOD   decay   rate    constant, 
day""". 


Thus,    the  mixed    BOD   concentration    in   the   river   in    terms   of 
the  ultimate   BOD  first-stage   demand   is   given   by 


where 


a  J 


o'c  1    -    exp      (  -    120    KO 


:■^      =      first-order   BOD  decay  constant,    hours 


-1 


(7.3) 
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The  other  initial  condition  required  is  the 
initial  oxygen  deficit,  Dq.   It  is  assumed  that  all  waste 
inflows  will  be  at  saturation.   Thus,  the  only  contribution 
to  the  initial  deficit  will  be  from  the  upstream  portion  of 
the  river.   Thus, 


D  Q 
u^u 


Q..  +  Q^  +  Q  +  Q- 


(7.4) 


where 


initial  D.O.  deficit ,  mg/1  ,  and 

D.O.  deficit  in  receiving  waters 
upstream  of  inflow  point,  mg/1 


7.3  Oxygen  Balance  of  Pol  1 uted  Streams 


Pollutant  transport  processes  in  a  stream  system  may 
be  adequately  aproximated  by  the  one-dimensional  version  of 
the  classical  convective  dispersion  equation.   This  partial 
differential  equation  is  based  on  the  principles  of  conserva- 
tion of  mass  (continuity)  and  was  presented  in  Chapter  II  as 
equation  (2.1).   The  main  sources  of  dissolved  oxygen  in  the 
stream  are  atmospheric  reaeration  and  oxygen  production  by 
photosynthesis.   The  major  si  nks  include  carbonaceous  oxygen 
demand  (CBOD) .  nitrogenous  oxygen  demand  (NBOD) ,  benthal 
demand,  and  respiration  of  aquatic  plants.   When  these  sinks  and 
sources  are  represented,  and  all  stream  system  parameters  are 
assumed  spatially  constant  along  the  flow  axis,  the  governing 
equation  is  given  by 
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f  =  E§  -   Uf  ^  ^2  (^s-^) 


-  K^L  -  K^N  +  P  -  R^  -  B  (7.5) 

where   C  =  concentration  ofD.O,  in  the  stream,  mg/1  , 

E  =  longitudinal  dispersion  coefficient, 
ft2/hr, 

U  =  freshwater  stream  velocity,  ft/hour, 

K2  =  atmospheric  reaeration  coefficient, 
hours-1 , 

C  =  dissolved   oxygen  saturation,  mg/1, 

C  -C  =  dissolved  oxygen  deficit,  mg/1  =  D, 

K,  -  deoxygenation  constant  of  carbonaceous 
'    BOD,  hours- 1  , 

L  =  remaining  carbonaceous  BOD  concentration, 
mg/1  , 

K  =  oxidation  coefficient  of  nitrogenous  BOD, 
"   hours-1, 

N  =  remaining  nitrogenous  BOD  concentration, 
mg/1  , 

P  =  oxygen  production  rate  by  algal  photo- 
synthesis, mg/1-hour, 

R  =  algal  respiration  rate,  mg/1-hour,  and 

B  =  benthal  demand  of  bottom  deposits,  mg/1-hour. 

This  equation  assumes  no  diffusion  of  pollutants  through  the 
river  boundaries  (other  than  that  which  is  included  in  the 
source-sink  term)  and  is  best  suited  to  predict  concentra- 
tions relatively  far  downstream  from  the  point  of  waste  in- 
jection.  Since  critical  D.O.  deficits  usually  occur  several  miles 
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downstream  from  the  waste  source,  equation  (7.5)  is 
particularly  well  suited  for  such  predictions.   For  fresh- 
water streams,  the  advective  flux,  is  signi  fi  cantly  larger 
than  the  mass  flux  due  to  longitudinal  dispersion 
(Hydroscience,  1971).   However,  it  is  of  interest  to  present 
a  methodology  that  can  be  applied  to  tidal  rivers,  and  thus 
longitudinal  dispersion  is  considered.   For  steady-state 
analysis,  all  system  parameters  are  assumed  constant  in 
time.   Since  it  is  desired  to  solve  for  the  D.O.  deficit 
and 


i£  =  0  .  11=    ID. 
3 1     '  9  X  "ax 


(7.6) 


equation  (7.5)  reduces  to  an  ordinary  differential  equation. 


dx2 


U  H7^  ^1  L  -^  ^N  -  h^ 


(P  -  R3  -  B)  =  0 


(7.7) 


A  basic  assumption  of  the  model  ,  stated  earlier,  is 
that  the  effects  of  the  biological  processes  (P  -  Rg  -  B) 
are  incorporated  into  the  measured  upstream  D.O.  deficit. 
The  production  and  consumption  of  oxygen  by  these  processes 
after  waste  injection  are  assumed  negligible  since  no  data 
were  available.   Although  measurements  of  total  organic 
nitrogen  and  nitrate-nitrogen  for  upstream  sources  are 
available,  modeling  of  the  nitrogenous  BOD  was  not  possible 
because 
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1.  The  amount  of  organic  and  ammonia  nitrogen 
present  in  all  wastewater  inputs  is  un- 
known ,  and 

2.  The  oxidation  coefficient  of  nitrogenous 
BOD  is  also  unknown. 

By  applying  the  above  assumptions  and  expressing  the  remain- 
ing carbonaceous  BOD,  L,  in  terms  of  the  ul timate  carbon- 
aceous BOD  and  its  reaction,  the  governing  differential 
equation  for  dissolved  oxygen  deficit  is  reduced  to  the 
expression 


E  d^D  _  ij  ^^    K,     (LJ   e 
dx2     d^    ^    °^ 


mx 


K2D 


(7.8) 


where 


U 
2E 


x  >  0 


1  + 


4K^E 


7.4  Critical  Deficit  and  D.O.  Levels 

The  solution  to  equation  (7.8)  as  a  function  of 
time  since  release  is  given  by 


D  = 


^"-oK   h 


eJt  -  e9* 


+  D^  e 


gt 


(7.9) 


where 


D.O.  deficit ,  mg/1  , 


deoxygenation  coefficent,'  hours   , 
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=  reaeration  coefficient,  hours"  ,  and 
lapsed  time,  hours. 


J        or 


2E 


Ml 
2E 


1 
1 

-v.. 

4KiE 
U2      J 

-v.. 

4K2E 
u'     J 

hours"  ,  and 


hours 


-1 


To  determine  the  time  at  which  the  critical  (maximum)  deficit 
occurs,  the  partial  derivative  of  the  D.O.  deficit  equation 
(7.9),  is  taken  with  respect  to  time  and  set  equal  to  zero 


(It=°) 


0  =- 


h   "  "^1 


JeJ^c 


ge^  c 


+  gD^e^^c      (7.10) 


Solving   for  t     the  following  expression  is  obtained 


wu 


{   3.  )    +   In 


'^l^^^c-Pp^-^Po^l' 


(7.11) 


Equation  (7.11)  may  be  simpl i f i ed  for  convenience  by  making 
some  substitutions, 


where 


elapsed  time  at  which  the  critical  deficit 
occurs ,  hours , 

=  self-purification  ratio 


Kp/  K ,  ,   and 
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ratio  of  the  initial  D.O.  deficit,  D 

0 

to  the. initial  BOD,  (L.)-,  d i men  si  on  1  ess. 


Finally,  the  critical  deficit  is  found  by  substituting  the 
value  of  t^,  given  by  equation  (7.12),  into  equation  (7.9); 


D. 


(L  )  K, 
V  o^c  1 

K2  -  K^ 


tc  _  gS^c] 


+  D^e^tc 


(7.13) 


where 


D   =   critical  (maximum)  deficit,  mg/1 


The  minimum  D.O.  level  is  calculated  as 


C  .   -  e  -  D 
mi  n    s    c 


(7.14) 


where 


min 


concentration  of  D.O.  at  maximum 
deficit ,  mg/1 ,  and 

saturation  concenration  of  D.O.  , 
mg/1  . 


The   saturation    concentration    is    determined    from   the 
regression    relationship    (American  Society  of  Civil   Engineers,  1960), 


Cg  =  14.652  -  O.41022T  +  0.0079910T^  -  0.000077774T^ 


(7.15) 


where 


water  temperature,  °C 
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7.5   Reaeration  and  Deoxygenation  Coefficients 

The  deoxygenation  coef f ici ent ,  K, ,  represents  the 
loss  of  D.O.  in  the  water  due  to  reduction  of  BOD. 
A  calibrated  value  of  0.7  day'^  is  used  for  this  simulation 
for  K-j  at  20°C  (68°F).   A  temperature  correction  and 
conversion  to  hour"   gives 


K^(T)  =  jj   K^  (20°)  1.047 


T-20 


(7.16) 


A  variety  of  formulas  exists  for  prediction  of 
the  reaeration  coefficient  K2,  almost  all  of  which  depend 
upon  velocity,  U,  and  depth,  H.   The  equation  of  Langbein 
and  Durum  (1967)  was  chosen  because  it  is  most  closely 
related  to  subsequent  procedures  used  to  obtain  U  and  H. 


where 


K^  =  2.303 


3.3 


,1  .33 


(7.17) 


reaeration  coefficient  at  20°C,  day"\ 

stream  velocity,  ft/sec,  and 
stream  depth,  ft. 


The  problem  1 ies  in  obtaining  values  of  U  and  H, 
since  the  stream  flow  varies  with  time.   In  the  absence  of 
measurements,  or  if  the  data  cannot  be  obtained  in  an 
expedient  manner  (as  in  the  ensuing  application  to  the  Des 
Moines  River),  an  approximation  can  be  made  based  on  the 
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work  of  Leopold  and  Maddock  (1953)  in  which  they  show  strong 
correlations  between  velocity  vs  flow  and  depth  vs  flow, 
namely, 


U  =  a^Q^ 


H  =  e.^Q 


where 


"l >a2'^i 


streamflow,  cfs,  and 
regression  coefficients 


The  model  presently  utilizes  coefficients  which  were 
determined  for  the  Kansas  River  System  in  Kansas  and 
Nebraska,  for  which 


(7.18) 
(7.19) 


1 

=  ■ 

1  .60 

2 

= 

0.03 

1 

■  = 

0.11 

2 

= 

0.45 

(7.20a) 
(7.20b) 
(7.20c) 
(7.20d) 


When  equations  (7.18),  (7.19)  and  (7.20)  are 
substituted  into  (7.17)  and  conversion  is  made  to  units  of 
hour   ,  the  reaeration  coefficient  is  established  as  a 
function  of  streamflow,  Q, 


2.303  [4.14350'°-^'']  1  .024^"^° 


(7.21) 


where 


stream  temperature,  °C,  and  the  last 
factor  represents  a  temperature  correction 
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7.5  Special  Probi ems 

_  o '  ■  .  "  ■ 

During  the  model  application  to  Des  Moines,  various, 
problems  were  encountered  revolving  around  the  critical 
dei'icit(D^)  and  critical  time  (t^)  equations,  equations 
(7.13)  and  (7.12),  respectively.   Due  to  the  large  number  of 
conditions  being  simulated,  situations  were  encountered  in 
which 

1.   the  deficit  load  ratio,  R^,  was  undefined 


o'  c 


(7.22) 


because  both  D^  and  (Lq)_  were  equal  to 
zero, 

2.  the  self-purification  ratio,  f,  was  equal 
to  one,  causing  equation  (7.12)  to  be 
undefined,  and 

3.  values  of  R^  were  such  that  negative  values 
of  t^  were  obtained. 

Mathematical  analysis  led  to  the  incorporation  of  certain 
modifications  and  safeguards.  Thus,  equations  (7.12)  and 
(7.13)  were  defined  only  for 

1  .   f  f      1  ,  or  g  j«  j, 

2.  (L^)^  f      0,  and 

3.  0  <  Rq  <  1/f . 

Otherwise,"   in    order,    if 
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1 •   f  =  1 .  or  g  =  j ,  then 


Dr  =  (•-)/  e  '^°"^  (7.23) 


c     o'c 


2.   f  ,  j«   1,  Lq  =  0,  then 


Dc  =  ^0  (7.24) 


3.   f  7*   1  ,  L^  =  0,  R^  >  1/f ,  then 


These  equations,  obtained  by  taking  limits,  are  not 
particul ar  to  Des  Moines  and  are  applicable  to  any  receiving 
stream  (Arbabi   ejt  a  1 .  ,.  197  4). 

7 . 7  Application  to  the  Pes  Moines  River 

A  description  of  the  receiving  stream,  urban  area 
waste  sources,  and  field  measurement  efforts  has  been  pre- 
sented in  Chapter  III,   Methodology  has  been  discussed  in 
Chapters  II,  III,  IV,  V  and  the  previous  sections  of  this 
chapter.   Before  presenting  the  results  of  the  receiving 
water  analysis,  it  is  important  to  discuss  model  validity. 
The  verification  procedure  was  preceded  by  calibration  of 
the  urban  runoff  BOD^  loading  rates  for  Des  Moines,  Iowa, 
as  computed  by  STORM.   The  dust  and  dirt  surface  loading 
factors  were  adjusted  to  obtain  an  annual  average  BODc  con- 
centration of  53  mg/1  for  urban  stormwater  runoff.   The 
above  concentration  was  the  average  value  determined  by 
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the  field  monitoring  program  in  the  separate  sewer  system. 
The  developed  mathemati cal  model  ,  as  discussed  in  the 
methodology,  simulates  the  mixing  of  stormwater  runoff  and 
sanitary  sewage  in  the  combined  sewer  system.   The  annual 
average  BOD^  concentration  of  combined  sewer  overflows  was 
computed  to  be  75  mg/1  ,  including  the  effects  of  first  flush, 
The  average  value  determined  by  the  field  monitoring  program 
in  the  combined  sewer  system  was  determined  to  be  72  mg/1 
(Davis  and  Borchardt,  1974). 

Verification  Analysis 

An  important  part  of  the  total  effort  required  to 
develop  a  mathematical  model  of  water  quality  in  a  stream  is 
devoted  to  verification  and  improvement  of  model  accuracy. 
The  verification  procedure  recommended  for  steady-state 
water  quality  models  includes 

1 .  examination  of  model  output  using  preliminary 
coefficients  on  a  diverse  set  of  data  (dif- 
ferent waste  loads  and  temperatures  under 
conditions  of  high  and  low  flow,  and  vari- 
able initial  stream  quality); 

2.  assessment  of  the  closeness  of  fit  of  observed 
field  data  to  computed  values; 

3.  adjustment  of  the  model  coefficients  until  the 
desired  accuracy  is  obtained;  and 
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4.   achievement  of  a  mathematical  abstraction 
that  reasonably  reproduces  observed  stream 
response  and  establishes  the  necessary 
validity  for  planning  purposes. 
The  carbonaceous  BOD  reaction  coefficient,  K, ,  was 
refined  during  the  verification  process  to  a  final  value  of 
0,70  day"   (at  20°C).   The  model,  of  course,  converts  to 
units  of  hour"  and  adjusts  for  temperature  through  equation 
(7.15).   The  atmospheric  reaeration  coefficient,  Kj,,  is 
calculated  internally  as  a  function  of  streamflow  and 
temperature  by  equation  (7.21);  therefore,  no  adjustment  Was 
necessary.   Since  field-measured  stream  dissolved  oxygen 
concentrations  were  given  as  daily  values  (Davis  and 
Borchardt,  1974),  storm-averaged  D.O.  concentrations  were 
determined.   Measured  daily  values  and  storm-average  computed 
values  of  D.O.  at  a  point  5.6  miles  (9.0  km)  downstream  from 
the  confluence  of  the  Raccoon  and  Des  Moines  Rivers  are 
compared  in  Figure  7-1.   Correlation  between  the  calculated 
and  observed  profiles  is  quite  good.   The  point  corresponds 
to  sampling  location  No.  6  as  shown  previously  in  Figure  3-2. 
Included  in  Figure  7-1  are  rainfall  and  average  total 
river  flow  values  for  each  wet  weather  event  (as  defined  in 
Figure  4-3).   Differences  between  measured  and  computed  D.O. 
concentrations  may  be  attributed  to  such  factors  as    (1)  the 
time  of  day  during  which  the  sample  was  taken,.  (2)  the  lag 
time  between  sampling  and  laboratory  analysis  and  the 
temperature  variations  in  the  receiving  water  during  the  day. 
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and  (3)  a  lack  of  data  on  photosynthesis ,  algal  respiration, 
and  benthic  demand.   The  time  scale  in  days  represents  the 
wet  year  beginning  on  March  8,and  ending  December  30,  1968. 
Again,  it  should  be  emphasized  that  these  D.O.  values  are 
not  the  minimum  D.O.'s  resulting  from  maximum  deficits .   The 
maximum  deficits  occur  much  further  downstream  and  water 
quality  standards  are  violated  much  more  frequently. 

A  longitudinal  dispersion  coefficient,  E,  of  50 

2 
ft  /sec  (16,722  sq  m/hour)  was  assumed  for  the  Des  Moines 

River.   Verification  was  not  possible  by  field  measurement. 
However,  its  order  of  magnitude  is  within  the  range  of  a 
moderate,  free-flowing  freshwater  stream  (Hydroscience,  1971) 
For  such  flow  regimes,  the  advective  flux  is  much  larger 
than  the  dispersive  mass  flux.   Nevertheless,  a  stretch  of 
river  may  undergo  a  transition  from  a  freshwater  stream  with 
a  high  velocity  (during  the  high  flow  season)  to  a  semi- 
tidal  river  with  low  velocity  during  the  low  flow  period 
(Hydroscience,  1971).   It  is  important  to  compile  data  on 
longitudinal  dispersion  for  large,  deep  rivers  where  the 
effect  of  dispersive  mass  flux  can  generally  not  be  neglected 
for  time-varying  waste  inputs  (Thomann,  1973). 

Results 


A  review  of  Table  7-1  is  appropriate  at  this  stage 
The  urban  waste  load  to  the  receiving  stream  may  be 
significantly  controlled  by  the  community  itself  through 
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(1)  the  treatment  rate  given  municipal  wastewaters  at  the 
DWF  plant,  and  (2)  the  size  of  the  WWF  storage/treatment 
facility,  which  determines  the  residence  time  of  urban 
stormwater  inflows.   Receiving  water  response  is  given 
for  the  combined  effects  of  urban  waste  sources  and  up- 
stream waste  sources.   Essentially  by  varying  the  degree 
of  control,  as  shown  in  Table  7-1,  10  different  combina- 
tions of  waste  inputs  are  investigated  for  the  stretch 
of  the  Des  Moines  River  extending  downstream  from  its 
confluence  with  the  Raccoon  River  (refer  to  Figure  3-1). 
in  addition,  the  response  of  the  receiving  water  is  shown 
for  the  unlikely  situation  that  no  treatment  whatsoever 
is  provided  by  the  urban , environment,  as  a  basis  for 
comparison. 

Results  are  presented  in  the  form  of  minimum 
dissolved  oxygen  cumulative  frequency  curves,  as  shown 
in  Figures  7-2  and  7-3.   The  ordinates  represent  the  per- 
centage of  critical  D.O.  concentrations  exceeding  the 
given  stream  standard  (abscissa)  during  the  wet  weather 
period.   The  critical  (minimum)  D.O.  concentration  in 
the  stream  is  computed  for  e^ery   hour  of  wet  weather 
during  the  year.   Thus,  the  ordinates  also  represent 
the  percentage  of  all  wet  weather  hours  that  exceed  a 
particular  dissolved  oxygen  level.   The  dashed  curve 
labeled  NO  TREATMENT  in  both  figures  represents  .the  re- 
ceiving water  response  for  the  unique  case  in  which 
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Figure  7-2.   Mi nimum  D.O. Frequency  Curves  for  Varied  WWF  Control 
and  DWF  Secondary  Treatment  Versus  No  Treatment 
of  Urban  Waste  Sources. 
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CONSTANT  VOLUME  STORAGE  BASIN 
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Figure  7-3.   Minimum.  D.O.  Frequency  Curves  for  Varied  WWF 
ControT!;;and  DWF  Primary  Treatment  Versus  No 
Treatment  of  Urban  Waste  Sources. 
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neither  the  municipal  wastewaters  nor  the  urban  stormwater 
flows  are  given  any  degree  of|control.   The  State  of  Iowa 
enforces  a  dissolved  oxygen  stream  standard  which  consists 
of  (1)  maintaining  a  minimum  dissolved  oxygen  concentra- 
tion of  at  least  4.0  mg/1  at  all  times,  and  (2)  maintain- 
ing a  D.O.  concentration  of  at  least  5.0  mg/1  for.  at 
least  16  hours  a  day  (State  Hygienic  Laboratory,  1970). 
From  Figures  7-2  and  7-3,  it  is  possible  to  determine 
the  percentage  of  hourly  minimum  D.O.  concentrations  that 
meet  the  required  levels  for  all  wet  weather  events  during 
the  entire  period  of  simulation,  March  8  to  December  30, 
1968.   There  was  a  total  of  449  hours  of  wet  weather 
during  this  period,  producing  a  total  runoff  of  10.28 
inches  (261  mm). 

It  can  be  concluded  from  Figures  7-2  and  7-3  that 
changing  the  DWF  treatment  rate  from  secondary  to  primary 
reduces  the  incremental  benefit  derived  from  maintaining 
higher  detention  times  in  the  WWF  storage/treatment  facil- 
ity, especially  at  the  higher  D.O.  levels.   However,  it 
is  also  clear  that  receiving  water  response  is  much  more 
sensitive  to  variations  in  the  WWF  storage/treatment 
detention  times.   The  effectiveness  of  WWF  control,  while 
providing  secondary  treatment  to  DWF,  is  summarized  in 
Table  7-2  and  compared  to  the  option  of  no  control  of 
urban  wastewaters.   It  appears  that  if  the  goal  of  main- 
taining an  absolute  minimum  D.O,  level  of  4.0  mg/1  in 


245 


Table  7-2 


Percentage  of  Critical  D.O.  Concentrations 
Exceeding  Given  Stream  Standard* 


Degree  of 
Control 


Stream  Standard 
of  4.0  mg/1 


Well -mixed 
Constant  Volume 
WWF  Storage/ 
Treatment 


Stream  Standard 
of  5.0  mg/1 


1 


•D 


0,  85% 


DWF  Treatment 


82% 


77% 


85%  DWF 
Treatment 


85 


80 


t  r,  =  6  h  r  s  , 

85%  DWF 
Treatment 


88 


82 


tp  =  12  hrs, 

85%  DWF 
Treatment 


91 


86 


tjj   =   24   hrs, 

85%    DWF 
Treatment 


95 


91 


No  Control  of 
Urban  Wastewaters 


77 


68 


'During  period  of  wet  weather 
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the  stream  is  to  be  achieved,  a  24-hour  detention  facil- 
ity will  be  needed.   For  practical  reasons,  it  is  likely 
that  an  urban  community  would  have  to  build  a  number  of 
these  facilities  rather  than  a  very  large  central  holding 
basin.   Nevertheless,  the  results  shown  in  Figures  ,7-2 
and  7-3  and  Table  7-2  indicate  the  need  for  the  higher 
detention  times  in  the  WWF  control  facilities.   It  is 
anticipated  that  further  research  into  storage/treatment 
of  urban  stormwater  will  be  needed,  but  results  indicate 
that  if  a  storage  facility  maintains  adequate  detention 
of  wet  weather  flows,  equalization  and  first-order  decay 
effects  are  sufficient  to  provide  the  degrees  of  control 
to  prevent  extended  violations  of  receiving  water  D.O. 
standards . 

It  is  interesting  to  note  that  the  mathematical 
model  predicted  stream  velocities  in  the  Des  Moines  River 
ranging  from  1.9  to  2.1  ft/sec  (0.58  to  0.64  m/sec).   At 
these  velocities,  it  can  be  shown  that  the  advective  flux 
is  more  significant  than  the  dispersive  flux  by  several 
orders  of  magnitude.   At  the  lower  levels  of  WWF  control, 
critical  distances  up  to  about  50  miles  (80  km)  were 
predicted.   By  contrast,  at  the  higher  levels  of  WWF 
control  critical  distances  up  to  about  20  miles  (32  km) 
were  predicted. 

Unquestionably,  a  detention  facility  for  urban 
stormwater  flows  has  a  favorable  impact  on  receiving  water 
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quality.   The  Humboldt  Avenue  combined  sewer  overflov/ 
detention  tank  (City  of  Milwaukee  et  al- ,  1975),  con- 
structed as  a  demonstration  project  to  intercept  over- 
flow from  a  750-acre  (231-ha)  segment  of  the  urban  area, 
provides  a  good  example.   The  detention  facility  had  a 
capacity  of  3.9  million  gallons  (14,762  cubic  meters). 
An  analysis  of  the  data  obtained  during  the  12-month 
period  between  November  1,  1971,  and  October  31,  1972, 
indicated  that  the  project  facility  prevented  approxi- 
mately 100,000  pounds  of  BOD  (45,359  kg  BOD)  out  of 
147,000  total  pounds  of  BOD  (66,578  kg  BOD)  from  being 
discharged  to  the  Mi Iwaukee  River- (City  of  Milwaukee 
ejt  a^.  ,  1975).   This  indicates  a  percentage  removal 
efficiency  of  68  percent  for  BOD.  The  Humboldt  Avenue 
detention  tank  may  be  compared  to  a  constant  volume 
storage/treatment  system  providing  approximately  20  hours 
of  detention.   Direct  effects  on  river  quality,  however, 
could  not  be  assessed  adequately  due  to  several  factors 
exerting  significant  but  undetermined  influence,  such 
as  (1)  algal  and  macrophyte  activity,  and  (2)  benthal 
oxygen  demand  from  accumulated  sediments  (City  of  Mil- 
waukee et  al . ,  1975) . 


CHAPTER  VIII 
SUMMARY  AND  CONCLUSIONS  . 

A  unified  concept  of  pollutant  mass  transport 
through  the  various  components  of  the  physical  system 
has  been  presented.   Mathematical  models  derived  from 
the  unifying  principle  of  continuity  have  been  developed 
and  successfully  applied  to  an  actual  urban  area  config- 
uration (Des  Moines,  Iowa)  and  its  receiving  water. 
Deterministic  and  statistical  techniques  have  been 
blended  to  serve  complementary  roles.   Thus,  the  inter- 
action of  urban  runoff,  control  measures  and  receiving 
water  response  has  been  demonstrated.   Since  continuous 
simulation  was  performed  for  all  wet  weather  hours 
during  the  precipitation  year  of  record,  frequency 
analyses  of  BOD  concentrations  and  mass  rates  into  and. 
from  the  WWF  storage/treatment  systems,  as  well  as  of 
minimum  D.O.  concentrations  in  the  receiving  stream 
resulting  from  these  BOD  inputs,  were  possible. 

Although  the  product  Ktp.  may  be  a  useful  indi- 
cator of  treatment  potentiaT  in  the  various  components 
of  the  urban  and  natural  environments  (e.g.,  sewer  pipe 
segments,  DWF  treatment  plants,  WWF  detention  facilities, 
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rivers,  lakes,  and  estuaries),  the  detention  time  tr,  has 
been  shown  to  be  the  most  critical  parameter  in  storage/ 
treatment  systems  subjected  to  urban  runoff  inflows. 
Mathematical  models  were  developed  for  a  well-mixed  con- 
stant volume  system,  a  well-mixed  variable  volume  system, 
and  a  dispersive  variable  volume  system.   In  all  cases, 
significant  reductions  in  peak  and  mean  BOD  loads  to  the 
receiving  body  of  water  were  predicted  at  the  higher 
levels  of  detention  (6  to  24  hours).   Results  obtained 
from  the  well-mixed  constant  volume  model  compared  closely 
to  the  performance  of  the  Humboldt  Avenue  combined  sewer 
overflow  detention  tank  in  Milwaukee,  Wisconsin.  . 

The  receiving  water  response  was  evaluated,  in 
terms  of  critical  deficits  and  resulting  minimum  D.O. 
concentrations,  for  the  case  where  the  WWF  storage/ 
treatment  system  was  represented  by  the  well-mixed  con- 
stant volume  model.   The  results  indicated  that  the  WWF 
control  device  was  effective  in  improving  minimum  D.O. 
levels  in  the  stream  for  average  detention  times  of  2,  6, 
12,  and  24  hours.   From  inspection  of  the  frequency 
analyses  of  BOD  concentrations  and  mass  rates  into  and 
from  the  other  storage/treatment  systems  modeled,  and 
the  pollutant  removal  efficiencies  predicted,  it  appears 
reasonable  to  conclude  that  the  variable  volume  devices 
would  perform  equally  well  or  better. 
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in  order, 


Several  recommendations  for  further  research  are 

(1)  Additional  modeling  needs  include  a 

more  detailed  representation  of  DWF  treat- 
ment processes,  and  extension  of  the  re- 
ceiving water  analysis  to  the  variable 
volume  WWF  storage/treatment  systems; 

(2)  Concepts  from  autocorrelation  analysis 
and  the  correlogram  may  be  useful  to  fur- 
ther characterize  the  response  of  the 
different  storage/treatment  systems;  and 

(3)  More  extensive  verification  of  the  perform- 
ance of  these  storage/treatment  systems  may 
be  possible  if  demonstration  projects  are 
financed  to  construct  prototypes  in  other 
urban  areas . 


APPENDICES 


APPENDIX  I 


WELL-MIXED  AND  DISPERSIVE  VARIABLE  VOLUME 
DIGITAL  COMPUTER  MODEL 
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